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PREFACE 


The work presented in tliis thesis is tlie out come of research is earned 
out in the field of queueing theory, at the department of Bundelkhand 
University Jhansi, while engaged in teaching as a Lecturer, in the 
department of Mathematical science & computer Application, Bundelkhand 
University, Jliansi. 

The work presented in die thesis is based on my following research 
papers. 

1- The steady state model of the machine interference 
M/M/C/K/oD(with balking, reneging and spares) 

published in XIX annual conference of the mathematical society, 
Banaras Hindu University. 

2- Response time analysis to the bulk arrival queue system in 
communication network. 

Published in proceeding of Fifth conference of the International 
academy of Physical Sciences, Bundelkliand University Jhansi. 

3- The transient behaviour of the machine interference : M/M/C/K/N 
(with balking, reneging and spares) 

Presented in Nintii annual conference of Gwalior academy of 
mathematical sciences. Institute of Technology & management 
sithouli, Gwalior. 


4- Transient behaviour of the M/M/C interdependent Queueing 
model with controllable arrival rates. 

Presented in national conference on Information Technology, 
operations Research and computing, Dau Dayal Institute of 
vocational education Dr. B.,R. Ambedkar, University, Agra. 

5- The steady state analysis of M/GVl retrial queue with multiple 
vacations 

Presented in tire tenth Annual conference of Gwalior academy of 
Mathematical sciences, Govt. Kamla Raja Girls post Graduate 
(Autonomous) college, Gwalior. 

This thesis consists of six chapters, each divided into several sections 

(progressively nmnbered 1.1, 1.2, etc,) The results in the text have 

been numbered serially, section and chapter wise, e.g. (3.2.1) means first 
result of second section of chapter third. The references cited have been 
arranged alphabetically and year wise at the end of thesis wherever more 
than one paper of the same author occur. The same year, they are 
distinguished by letters a,b,c.... etc. The references are indicated in the text 
by giving the year within parenthesis. Tliis work starts with tlie study of the 
steady state model of the machine interference with balking, reneging and 
spares. We analyze two cases of spares first is c>y and c<y and Mso discuss 
special case. 

Chapter I, deals with introduction and brief historical survey of the 
work, done in the field of Bulk Queueing model, with special references to 
the work embodied in the thesis. 
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Chapter 11, is devoted to the study of the steady state model of the 
machine interference, studied by A.L showky [133], In this chapter we have 
analyzed the steady state behaviour of the system and obtained length 
Queue of each cases. 

Chapter HI, is devoted to the transient behavior of the machine 
interference studied by A.L Shawky [133], In this chapter we have discuss 
the transient state behaviour of the system and assume die boundary 
condition then we have obtained length queue of each cases. 

Chapter IV, is devoted to transient behaviour of interdependent 
queueing model with controllable arrival rates studied by M.I. Aflab Begum 
and D.Maheswari [107]. In this chapter we have analyzed the transient state 
behaviour of the system and its laplace transform then we have obtained the 
Prob. that the system is in state, r and R, empty, conditional Probability 
state, 1, Mean Queue length and expected waiting and particular cases. 

Chapter V, is devoted to the steady state analysis of Bulk Service 
retrial queue with multiple vacations studied by B. BCrishna Kumar and S. 
Pavai Madneswari [12]. The supplementary variable techniques and 
boundary condition have been used in this chapter, then we have obtained 
the probabihty generating function of the orbit size when server is idle, 
busy and under repair. 

Chapter VI, is to devoted to Response time analysis to the bulk 
arrival queue system in communication network studied by M.P. Singh V. 
Kumar and A. Kumar [100]. In this chapter we have analyzed mean waiting 
time, mean response time and its computation of numerical residts. 


(Anis Ahmad) 
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INTRODUCTION 


Queue is an essential part of life as it has been seen 
everywhere whenever we face a queue , we have to wait till the 
clearence of the queue. There are various situations, such as van, 
post office. Railway reservation counter. Transportation system. 
Communication system etc. are some of the common queues that we 
find in our daily routeins. Queueing theory is one of the most 
important branch of applied probability theory and has vast 
application in technology and management. As far as the future 

development of technology and management may be extrapolated 

ft 

from the past and present state of affairs the need for a deeper insight 
into queueing theory will increase rapidly. It is hardly necessary to 
point out the actual queueing situations encountered in every day 
life. Production lines, the theory of scheduling and transportation, the 
designs of automatic equipment such as telephone and telegr^h exchanges 
and particularly the repidly growing field of information handling and 
data processing are new fields in which queueing situations are 
encountered.The beginning of the technological development of automatic 
telephone exchanges led to infinity number of problems which could be 
solved in a satisfectory way only by probabilistic method. 

It is Erlang (cf BROCKMEYER, HALSTROM, JENSEN [1948]), 
* A Danish mathematicion, who was the founder of queueing theory. His 
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studies , in the period of 1909-1920, are now classical work in queueing 
theoiy.Until about 1940 the development of queueing theory has mainly 
been directed by the needs encountered in the design of automatic telephone 
exchanges. After the second world war, when applications of mathematical 
models and methods in technology and organizations rose to a level hither 
to unknown, it was soon recognized that queueing theory has a lot of 
potential in most applications of wide variety of fields such as, machine 
maintenance, road traffic, aviation problems, defence operations, inventory 
management etc. It has been used to study physical phenomena like the 
semi-conductor noise problem (Bell, 1958), Panico (1969), have given 
various applications of queueing theory to business, economics and 
Science.The theory has applications to the following scientific areas. 
Catalysts in a chemical reaction, filtration process, gas molecules going 
through a whole, hospitals and the demand for medical Care, nervous 
reactions, and such psychological stress and strain as neurosis. 

In the queueing system the number of units (servers) are available 
to provide the services for the customers who arrive at their on choice of 
time. They are served by one or more servers according to a given service 
discipline, the service times being randomly distributed an that governed 
by some probabilitstic law and after being served the custmers depart 
from the system. A queueing system is specified by the following fectors: 

(1) The Input Process (Arrival Pattern of Coustomers) The 
Input process describes the way in which the customers arrive and join 
the system- Generally, the customers arrive in a more or less random ftishion 
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which is not worth making the prediction. Thus the arrival pattern can be 
^ described in terms of probabilities and consequently the probability 

destribution for interarrival times (the time between two Successive 
arrivals) or the distribution of number of customers arriving in unit time 
must be defined. Another fector of interest concerning the input process is 
the manner that arrivals come in batches or groups insteaxi of one at a 
time. In the event where more than one arrival come to the system altogether. 
The input is said to be bulk or batches type. Hie situation, where bulk 
arrival occurs the time between the successive arrivals of the groups also 
governed by the probabilistic law. 

(2) The service Mechanism (service pattern of servers) The time 
* interval from the instant of initiation of Service of a customer to the mstant 

when this service ends is called the service time. The service times Vj, 

Vj, of the successive units are assumed to be independent of each 

other and of the Input distribution, and their distribution B (v) is Called 
the Service time distribution or the holding time distribution. 

(If we can write dA(u)=a(u) du, dB(v) = b(v)dv then a (u) and b (v) are 
called the input and the service time density fiinctions respectively). 
Service may be single or in batch. The service rate may depend on the 
number of customers waiting for service. A server may work frister if the 
queue is building up or conversely he may get flustered and become less 
efficient. The situation in which service depends on the number of 
4 customers waiting is referred to as state dependent service. 
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Even if the semce rate is high, in general, customers arrive and 
depart at irreguler intervals hence the queue length wall assume 
pattern unless arrivals and service are deterministic. 

(3) Queue Disicipline The queue discipline is the order of manner in 
which customers in the queue will be served. The most common discipline 

is called 

(i) FCFS- First come first served, 

under which customer are served in order of their arrivals. 

(II) LCFS - Last come first served 

the last arrival in the system is to be provided service first of all. 

(III) SIRO- Service in random order 

in this queue discipline the arrivals are serviced randomly irrespective 

of their arrivals in the system. 

(iv) PSD - Priority service discipline 

under priority discipline, the service is of two types;- 

(1) Peeemptive Priority. In this case the customer with the highest 

priority is allowed to enter service immediately after entering into the 

system even if a customer with lower priority is already in service. That 

is, lower priority customer’s service is interrupted (preempted) to start 

service for a special customer. This service is resumed again after the 

higher priority customer is served or according to the preemptive rules. 

(2) Non-Preemptive Priority. In this case highest priority customer 

goes ahead of the queue but service is started immediately on completion 

ofthe current service. ■* 




(4) System capacity. In queueing process there is a physical linritation 
queue ofthe waiting room, so that when the queue reaches a certain leng*, 
no further customer is allowed to enter until space becomes available by 
a service completion. These ate referred to as finite queuemg situationsp.e. 


there is a finite limit to the maximum queue size . A queue with limited 
waiting room can be viewed as one with forced balking where a customer 
is forced to balk if he arrives at a time when queue size is at its limit. 


(5) The number of channel:- The service fecility can have one or more 
channels queueing processes with single channel is called single channel 
or single server queueing process, while those with more than one channel 


is called multichannel or multi server queueing process. 

Classification of Queueing Model:- Some ofthe weU known queue 

models may be classified as ; 

(i) (MM/l): (OO /FC FS):- The queueing model (M/M/1) : (OO /FCFS) 
deals with a queueing situation having poisson arrrivals (exponential 
inter-arrival times) and poisson services (exponential service tunes), smgle 
server, infinite capacity ofthe system and first come first served queue 

discipline. 

(ii) (M/M/l) : (N/FCFS) :- The queueing model (M/M/1):(N/FCFS) is 
different firom the queueing model (M/M/l):(oO /FCFS) in respect ofthe 
Capacity ofthe system. Here the capacity of the system is limited to N 
customers in the system while remaiing fects are same. 
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(iii) (MM/S) : (00 /FCFS) "nie queueing model (MM/S) : (oo /FCFS), 

as mentioned inlWsmodel(MAl/l):(cX)/FCFS), while this model instead 

of a single service channel, this model deals with, multiple servers, m 
parallel which are to limited S. For this queueing system it is assumed 
that customers arrive at a mean arrival latex and are served accordi g 
first come first serve basis at any of these servers and it is also assumed 
that only one queue is formed and each server has same mean service rate. 


(iv) (MMJS) : (N/FCFS) The queueing model(M/M/S) ; (N/FCFS) is 
different fi-om model(NVM/S) : (oO /FCFS) in respect of the capacity of 
the system. Here the capacity ofthe system is limited otN customers only 
and first come first served discipline. 


(V) (M/M/S) (M/GD) Machine servicing Model:- The machine 
servicing model is same to the model (MM/S) : (N/FCFS) except in 
queue discipline, i.e. General Discipline (GD). The system has a finite 
..lling population, of size M (say) and the probability of arrivals depends 
upon the pastbehavior.The most common application ofthe queuemgmodel 
has been to the machine servicing problem in which say S repairmen 


(called servers) are assigned the responsibility for maintaining certain 
group of M machine. In this case, the calling population is the machmes 
and an arrival corresponds to a machine breakdown. 


(vi) (M/Ek/1):(oo/FCFS):- The queueing model (M/E^/ 1 ) : (oo/FCFS) 
consists of a single service channel in which there are K identical stages 
(Phases) in series for services i.e. poisson arrivals erlangian service time 
for K-phases and a single server. 
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(vii) Power Supply Model :-The power supply model may be usedlo'"^----- 
ensure better utilization of power supply to consumers . Suppose that the 
demand of power supply from the consumers follow poisson distribution 
with parmeter and existing supply schedule of power also follows poisson 
distribuation with parameter // . 

The following terms are used in this thesis : 

(i) Busy period distribution The busy period is defined as the interval 
of time from the instant of arrival of a customer at an idle channel to the 
instant when the channel next becomes idle for the first time .For this type 
of queue, the initial busy period may start with i>l customers and hence 
may be different from other busy periods which starts with the arrival of a 
single customer. Results for busy periods of single server, single -arrival, 
poison -input queue can rather easily be generalized to the corresponding 
bulk arrival queues.For multi- server queues any one of several different 
definitions of busy period may appropriate. 

Any queueing system is described by given first the input process, 
then the service time distribution and finally number of servers. For example 
- M/ E^/1 stands for a single server queueing system with poisson arrivals 
and K-Erlang service time distribution and GI / G / S stands for the most 
general case of general independent input, general service times and S for 
server. 

Prior to 1 950, Most of the work had been done under the assumption 
that the arrival or service time distribution, either exponential or K-Erlang 
or regular. Very little work had been done for general distribution. Notable 


workers before 1950 are Erlang (M/M/S,M/D/S, M/E^l), Vaulot (M/M/ 
S), Malina (M/M/S, M/D/S), Palm (M / M/S), crommelin (M/D/S; M/G/ 
S),Khintchine (M/D/ 1, M/G/ 1), Pollaczek (M/G/1, E^/G/1) M/M/S, M/D/ 

s, M/G/S), Borel (M/D/1), and volberg (E^/G/1) . A review of the work 
during this period has been given by kendall, Saathy, N.S. Kambo, MX. 
Chaudhary and Easton . 

(n) SUPPLEMENTARY VARIABLE TECHNIQUES:- This technique 
was used in 1 942 by Kosten this method consists in enhancing the description 
of the stated of the system by including one or more supplementary variables 
so as to make the process Markovian. Thus the queueing process M/G/1 is 
made Markovian by including the variable x, the elapsed serviced time of 
the unit being serviced ,and define P^ (x,t) dx as the probability that at time 

t, the queue length is n and the elapsed service time of the unit under service 
lies between x and x+dx. Similarly the Process GI/M/1 may be dealt with 
by the inclusion of the elapsed time for the last arrival. Kendall (1953) 
[97], who called it the augmentation technique. Formulation of the method 
and mentioned that the most general process GEG/S can be studied similarly 
by using s+ 1 supplementary variables technique denoting respectively the 
elapsed time since the last arrival and the elapsed service time of each of 
the s servers. 

The name appears due to Cox (1955b) [44] who used a supplementary 
variable technique to study the queueing system M/G/1, kendall (1953) 
[97] considered this technique, which he labeled ‘augmentation but preferred 



the use of an imbedded Markov’ chain as leading to simpler calculation. 
Colony (1958, 1959, 1960b) used this technique in studies ofthe process 
of the system GI/M/1 and to study the system GI/E^/1. Jaiswal (1961c, 

1 962) applied the technique to obtain transient solution of the permeative 
resume and the head ofthe line priority queues. The same technique has 
been used by thiruvengadum (1963, 1 964) to study queueing system subj ect 
to breakdown by thimvengadum and Jaiswal ( 1 963, 1 964a, 1 964b) to study 
Machine Interference problems and by Jaiswal (1965) [89] to study the M/ 
G/ 1 queue with balking and reneging Kashyap ( 1 965a, 1 965c, 1 967) used 
the supplementary variable technique to study the transient behaviour of 
the double ended queue with limited waiting space and in ( 1 966) extending 
the use of ^pplementary variable technique has studied head of the line 
priority queues with bulk arrivals. Kashyap (1972b 1972c) used the 
supplementary variable technique to study preemptive repeat priority and 
permeative queueing discipline. 

Gaur and Kashyap ( 1 972 h) [69], using the supplementary variable 
method obtained the joint distribution ofthe queue later (1972 i) using the 
same method head ofthe line priority queueing system with general service 
has been studied to obtain an expression for the laplace transform of 
generating function of the joint probability distribution of the number of 
priority and non- priority units in the queue at time t and the number of units 
served in time t 


(III) HYPERGEOMETRIC FUNCTIONS Hypergeometric function 

is denoted by pFg(a 2 ,a 2 p, p^ ;z) or simply pFg when 

there is no confusion about the arguments of the function :F(a, b ; c ; x) or 
F [^*’’C ; x] or ^Fj [^'’’C ; x], Where the suffix 2 in ^Fj denotes the number of 
parameters in the numerator and the right sufiBx 1 in jF j denotes the number 
of parameters in the denominator of the series and the hypergeometric 
function does not change if the parameters a and P are interchanged keeping 
y fixed. 

F (a , p ;y ; x) = F (P , a ; y ; X ) and its deduction by definition 

(a)^ = a(a+l) .(a + n-1) 

(a)^, = a (a+1) „ 

(a+n) - (a)^= (a) 

(IV) CONTROLLABLE ARRIVAL RATES GeneraUy we consider a 
single server infinite capacity queuing system with the assumption that the 
arrival and service processes of the system are correlated and follows a 
bivariate poisson process. In addition to this interdependence it is also 
assumed that whenever, the queue size reached a certain prescribed munber 
R, the arrival rate reduces fi'om A^^^to and it continues with that rate as 
longer as the content in the queue is greater than some other prescribed 
integer r (r > 0 and r < R) when the content reaches R, the arrival rate 
charges back to and the same processes repeated, This type of operating 
strategy is called controllable arrival rates. In the controllable Arrival 
rates Lindly (1952) [104] analyzed a single server queueing model with 
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recurrent input and arbitrary service times with the assumption that the 
service time of the n^ customer and time interval epoch are related by a 
linear function. Conolly and Hadidi (1969) [38] have studied the single 
server queueing model with the assumption that the rate Sj^/T^is constant 
for all n , where S is the service time of the n*^ customer and T is the time 
interval between n* and (n+1)* arrival epoch, Michel and Paulson (1979) 
[114] have considered an M/M/ 1 queueing system with an assumption that 
the inter arrival times separating the arrival from that of his predecessor 
and the customer service times follows a bivariate exponential distribution. 
Lengaris (1987) [103] Borst and combe (1992) [25] have studied the busy 
period analysis of a correlated queue with exponential demand and service. 
BULK QUEUE The concept of bulk queues was implied in the work of 
Erlang. The study of queues with bulk service is originated in 1954 by 
Bailey. During the last four decades, queues with bulk service 
system with bulk service and finite capacity. The units arrive at a single 
service station and have received greater attention from researchers and 
practical scientists. 

A part from having the organized/ systematic mathematical structures, 
queues with bulk service have been used for modelling various situations 
like shiploading and unloading at a port dispatching of mail bags at central 
sorting stations, transhipment of members though elevators and at bus stop 
etc. In all these models are very interesting to note that the two major 
constitunt process play a vital role namely. 


(i) Arrival Process 

(ii) Service Process 

QUEUEING SYSTEM WITH BULK SERVICE Queues with bulk 
service are described as customers get service in batches. The size of the 
batches for service may be fixed or random variables. The service times 
of different batches are independent and identically distributed random 
variable. The maximum size of a service group is known as service 
capacity. 

Bailey (1954) [16] studied a queueing process in which customers 
arrive at random from a single queue in order of arrivals and are served in 
batches. The size of each batches has a fixed maximum. The Equilibrium 
distribution of queue length has been studied by an imbedded markov chain 
technique and its ergodicity was established when average demand is less 
than to the average supply available. Downton (1955, 1956) developed 
the waiting time distributions of these models and the limiting behaviour of 
the queue when the batch size tends to infinity Jaiswal (1960 a, 1960 b, 
1961) further studied these queueing models and considered the case when 
maximum number of units to be taken is not constant but depends on the 
number of units already present with the service as well as the capacity. 

Foster and Nyunt (1961) [58] considered a single server queueing 

system in which items arrive at the sequence of instants tj, t^ such 

that inter arrival times are inter independent randon variable with an 
exponential distribution with consideration that the items are in batches of 
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exactly k units by a single server with respect to arrival and discuss about 
the system characteristics of an E/GVl queue. Finch (1962) [57] studied 
the transient behavior of a queueing are served in batches of exactly x 
units, when x>l . The service times of successive batches are independently 
and identically distributed random variables with a common distribution, 

Foster and perera ( 1 964) [59] discussed about the queue size after a 
departure and considered the queueing problem earlier discussed by foster 
& Nyut ( 1 96 1 ) [58] under the assumption the units arrive in a poison process 
with parameter X and served in batches of fixed size k and the service time 
distribution being arbitrary. Foster (1 964) discussed a general single server 
queueing system when the units are single arrivals and batch departures as 
well as batch arrivals and single departures and compared the queue size 
with various instants that are just before arrivals just after departures and 
through time. 

Bhat (1964) [20] discuss about the bulk service queueing system, the 
arrivals arriving with a fixed time interval follows a binomials and die 
service time is a discrete variable assumed to be independently identically 
distribution. The busy period and the queue lengths are obtained in terms 
of probability generating function. Aurora ( 1 964) [9] had been studied the 
bulk service is exponential and obtained time dependent probabilities for 

V ■ 

the queue length of busy periods for at least one channel being busy and 
both channels being busy. 
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Neuts ( 1 965) [119] studied the busy period distribution of the queue 
in which the customers have been served m at a time. If there are less then 
m or more present and all at a time if there are less than m present and the 
units arrival in poisson process. 

Kashyap (1966) [94] considered the double ended queueing system 
for taxis and customer at a taxi stand under the assumption that there is 
limited waiting space both for customer and taxis and the customers are 
served in bulk. The arrivals of customers taken as general and taxis as 
poisson and used the supplementary variable technique in which the 
probability of the state is empty. The cases of erlangian and poisson arrivals 
of customer and of single stated taxies are studied as particular cases Roes 
(1966) studied the batch service Markov chain technique. Goyal (1967) 
[72] studied the bulk exponential service to obtain the steady state 
probabilities with generating function. 

Ghare (1968) [71] considered poisson arrivals and exponential 
service under bulk service pattern. This queueing problem is an extension 
to Aurora (1964) [9] for two channels, by using the generating functions 
and Laplace transformations and also obtained the transient probabilities 
and discuss about the characteristics of the system. Love (1970) [106] 
improved the steady state solution for the queueing system v^ere the arrivals 
are servied in batches by parallel servers of fixed size. He considered the 
inter arrival times as erlangian with a certain degree and the service time 
for all servers are exponentially distributed Stewart (1970) [136] considered 


M/G/ 1 batch service queueing system with balking and derived the steady 

- state probabilities and length ofthe queue. 

Murali (1972a , 1972 b) considered the steady state behaviour ofa 
discrete time, single server queueing problem in which the arrivals at two 
successive time marks are correlated the service is accomplished in batches 
and the completion of batches at two successive time marks are correlated 
through probability generating function. He obtained the number ofbatches 
waiting for the service. He also studied the transient behavior of a single 
server queueing system in which the arrivals follow apoisson distribution 
ofbatches of variable size and the output follows a general distribution of 
batches of variable size. He obtained the waiting line size when the service 
time distribution is exponential with fixed branches, phase type and 
exponential. Medhi & Borthokur (1 972) [110] analyzed the queueing system 
under three service mechanism. 

(i) The server will start the service only when a minimum number of 
customers present in the queue. 

(ii) If the number of customer in the queue lies between a and b the whole 
queue is taken into service as a bulk discipline. 

(iii) If the number of customers in the queue are more than b, the first b 
customers will be taken into service. 

Deb and serfeo (1973) [47] considered batch size discipline where 
each batch size and its time of service is controlled. To minimize the expected 
continuously discounted cost per unit time for service over an infinite time 
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horizen, they used the markov decision process i.e. dynamic programming 
approach. 

Ernest (1973) [52] considered the maximum and minimum times for 
service to all Jobs in a batch of random size with unknown number of 
servers and all the servers have the same service time distribution and also 
obtained the expected minimum batch service time for e>qK)nential, uniform 
service time distribution and poisson, binomial. Negative Binomial 
Probability mass functions for batch size. 

Villiam Makis ( 1 984) discussed the poisson arrivals and the general 
service time with a single server in which the capacity is infinite. He 
derived explicit formula for the long run expected average cost per unit 
time and properties of the average cost function and also determined the 
optimal control limit in the general case. Deb (1984) considered the 
compound poisson arrivals, to control a batch service queueing process B. 
Powell (1986, 1987) considered a queueing model in which the arrivals 
and departures are bulk on three strategies. The strategies are cancellation 
strategy holding strategy and combined vehicle holding and cancellation 
stretregy. He developed moments of file queue length and derivation of the 
queue length transform. He also obtained the distribution numerically for 
waiting times of queues with out control for vehicle cancellation strategy 
and vehicle holding strategies under the consideration of bulk arrivals and 
bulk service in the track load motor carrier industry. The mean and variance 
of the calculated distribution is strongly related to analytical expressions. 
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Inter dependent Queueing Models : - Bhat (1969) [21] Classified the 
dependent queues into five categories. They are : 

(i) Systems with correlated arrivals 

(ii) Systems with correlated services 

(iii) Systems with state dependent amvals. 

(iv) Systems with state dependent service 

(v) Systems with interdependent arrival and service process. 

Jacob and Patricia (1980) [81] considered the queueing models for 
a sequence of correlated exponential inter arrival and service times for a 
single server queue shun Cheii Niu (1981) analyzed a M/G/1 queueing 
system in which the service time required by a customer is dependent on 
the inter arrival time between his arrival and that of his predecessor using 
the concept of random variables developed by Esary. Prochan and wolkup 
(1967), Garg and khanna (1983) considered the steady state behaviour of 
a queueing system with queue length depending on additional server fecility. 

Rao ( 1 986) [1 26] analyzed a M / M / 1 queueing system based on 
bivariate poisson distribution with interdependent arrivals and services & 
the system behavior through dependence parameter and obtained explicit 
expressions. Rao (1989) [125] extended these models to inite source 
queueingmodel and queues with impatience. Gupta (1990) study a first 
come first served semi markov queue in which both the arrival and service 
mechanism are semi markov processes. The inter arrival and service times 
may depend on one another and the marginal service times is assumed to 
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be phase type the waiting time is spent in the system and virtual waiting 
time are obtained and analyzed. 

Rao (1991) considered the interdependent queueing models for 
analyzing the situations arising at communication networks, R^o (1983) 
developed interdependent queueing models for studying the performance 
measures of the resource sharing systems in computer and also obtained 
the optimal operating policies for these interdependent, queueing models 
having Erlangian service time using the dependence structure given by Rao 
(1986). Borstetal (1993) [26] had been studied the M/G/1 queue in which 
the service times of arriving customers depend on the length of the interval 
between their arrival and previous arrival. 

M/MVl Interdependent queueing model with Varying Batch Size 

Murty (1993) extended the single server interdependent queueing model to 
bulk service model. He consider a different problem which is a 
generalization of the earlier queueing models. Two variations are considered 
as follows : 

(i) The customer are served x at a time, except when less then x are in the 
system and ready for service at which all units are served and 

(ii) The batch size must be exactly x and if it is not, then the server remains 
idle until the queue size reached x, A typical situation of this model is 
transportation of men and material are the customer. For this short of situation 
in order to have optimal transhipment, it is appropriate to approximate the 
situations with the model having the interdependent arrival and service 
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process. For developing these interdependent models with bulk service 
rule we make use of the dependence structure given by Rao (1986) [126]. 

Let (t) be the probability that there are n customers in the system at 
time t. The difference differential equations of the model are - 

(t) = - (A. + p -2e) P/t) + (^ - e) P^,(t) + (p - e) P^(t) n > 1 (1.1) 

P'(t) = -(X-e)P„(t) + (n-e) tw (1.2) 

Where X and p are marginal arrival rate. Service rate and mean 
dependence rate respectively. 

Consider that the system reached the steady state, the transition equation of 
the model are:- 

- (k+ p- 2e) P^ + (X- e) P^j + (p- e) P^_^= 0 
-(X-e)P^ +(ff-e) 

hi 

Using the heuristic arguments of Gross and Harris (1974). They 
obtained the solution of these equations as 

P^ = cr " n > o, 0 < r < 1 (B) 

where r is the root which lie in (0, 1 ) of the characteristic ^nations 
[ (p -e) - (A. + p - 2e) D + (X - e)] P^ = 0 

HereDisthedifferentialoperator. 

The probability that the system is empty is 

Po = (i-r) (A) 

Where r is as ^ven in equations (A) 
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M/MVl Interdependent Queueing Model With Fixed Batch Size 

Murty ( 1 993 ) [ 1 1 7] Considered a slight variation of the model. Along 
with the other assumption, hence he assumed that the batch size must be 
exactly x, and if not, then the server waits until such time to start, in 
another case to reduce the idle time , it may be feasible to use 
interdependence relation for such type of models He assumed that the 
number of arrivals and the service of any batch are correlated and follows 
a bivariate poisson process with the dependence stmcture the difference 
differential equations of the model are 

P‘„(t) = (p - e) P'^,(t) - (X - p + 2e) P‘„(t) + (X - e) P‘„., (t),n>k (1.3) 
P'„(t) =(R-e)P'^(t)-(X-e)P'„(t) + (X-e) P',..(t),lin <k (1.4) 
P'„(t) = =-(>.-e)P',(t) + (n-e)P'„(t),n = 0 (1.5) 

V= + J_ 

12 (l-r)2 

Where r is as given in equation(B) 

The cofficient of variation of the model is 
= (l-r)[(k^-l)(l-r)+12r] 

3[(k-l) (l-r) + 2r] (1.6) 

Where r is as given in equation 

M/M*/l Interdependent Queueing model with variable service 
capacity In the bulk service queuing models Bailey (1954) [16] and 
Jaiswal (1960) [83] Assumed the model in which the units arrive at random 
from a single queue in order of arrival and are served in batches. The size 
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of each batch being either a fixed number of customer or the whole queue 
length. Which ever is smaller. Jaiswal (1961)[85] extended this model to 
this case, where at a service epoch m, (o < m <.s) if m persons are already 
present with the server then (s-m) persons or whole queue length. Which is 
smaller will be taken into service. The service rule is termed as Bailey 
bulk service rule . However in these models the arrivals and service 
processes are poissonian and Erlangian respectively with interdependent 
arrival and service process. 

In both the models the system behaviour is analyzed by obtaining the 
difference differential equations of the models and solving them through 
generating function techniques. The system characteristics like mean queue 
length. Variability of the system size, coefficient of variations are derived 
and analyzed in the light of the dependence parameter, these models are 
included in earlier models as particular cases for specific values of the 
parameters. 

Murty (1993) [117] considered the single server system with 
interdependence arrival and service process having the Bailey’s is bulk 
service with variable capacity Under the present study it is assumed that 

the server servers only at instants tj, t^ if m, (0 m ^ S) 

person are present in the waiting line at time t^then the server takes a batch 
of (s - m) persons or whole queue length which ever is smaller, vshere s is 
the service capacity. 
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Two server interdependent queue model with Bulkof Bailey’s service 

Muity ( 1 993 ) [ 1 1 7] considered two server independent bulk queueing model 
with Bailey’s service rule for analysis. The model is an extension of the 
models considered earlier or two server case, it is also considered that the 
bulk service process having two service facilities with capacity bj and b 2 
. A batch of bj units or the whole queue length whichever is smaller is 
taken from the head of the queue for service in the first channel when ever 
it is free similarly the second channel on becoming free takes b^S b, or the 
whole queue length which ever is less. 

If both the servers are idle and there is no queue the next unit to 
arrive always goes to the first service facility. Further we assumed that 
the two service facilities are independent of each other and the arrival and 
service process are independent. 

This sort of situation are more common in Marshalling yard’s with 
two engines. Elevator process with two lifts etc. This model is also extended 
to multiple server independent queueing models with Bailey’s Bulk service 
.In both the models he first develop the difference differential equation and 
solved them through generating fiinction techniques. 

Queueing system with Bulk Arrival The study of bulk have begun 
with Erlang solution of the M/Ej^/1 Since it gives implicit solution ofM*^/ 
M/1 . Explicit consideration of bulk arrival queues seems to have begun 
after several years of the work of Bailey (1954) [16] on bulk service In 
“ordinary” queueing problems it is assumed that customers arrive singly 
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at a service facility. However this assumption is violated in many real - 
world queueing situation. Letters arriving at a post office, ships arriving 
at a port in convey, people going to a theatre, restaurant, and so on are 
some of the examples of queueing situations where customers do not arrive 
one by one in bulk or groups, where the size of an arriving groups may be 
a fixed number i.e. Bulk or batch operations are common on the service 
side of some system. , 

A set of chapmen-kolmogorov (differential -difference) equations is 


dt 

^Pn(t) 

dt 


= -a + M)P. (t) + + ^ZP.K. (t)C, 

K=1 

= -^Po(0 + /'P, (t) 


n>l (1.7) 
( 1 . 8 ) 
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The Steady state Mode! of the Machine Interference 
M/M/C/K/oo (with balking, reneging and spares) 


INTRODUCTION : 

This chapter consists the queuing model M/M/C/K/ oo with balking, 
reneging and spares. C servers are available, customers arriving rate is X. 
We take Y spares so that when a machine fails, a spare is immediately 
substituted for it, if it happens that all spares are used and a breakdown 
occurs then the systems becomes short. When a machine is repaired it then 
becomes a spares (unless the system is short in which case the repaired 
machine goes immediately into service). We analyze the two cases of 
spares, first of them is C<Y and the other C>Y. 

Kleinrock [98] had worked out on the truncated poisson distribution 
for solving the problem M/M/C/K/N for machine interference but without 
assumptions of balking, reneging or spares. Gross and Harris [75] 
discussed the system : M/M/C/m/m with spares only and Medhi [112] 
treated the system M/M/C/m/m without assumption of balking, reneging 
and spares. Showky [133] considered M/M/C/K/N machine interference 
model with balking, reneging & spares with the development of technology 
of manufecturing, such as those found in con^uter and communicaticm 
system, the efficiency of machine system is critical to overall 
competitiveness. Many researchers have been worked on machine repair 
problem. 
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Gupta and Rao [76] derived a recursive method to compute the 
steady state probabilities of M/G/1 machine interference model Jain et. al. 
extended the work of Shawky [132] by including additional repairman in 
case of long queue of failed units. 

The main assumption of present chapter is an infinite source o 
infinite customers. Each with an arriving rate c servers are available, that 
the service tunes are identically exponential random variables with rate p 
also the system has finite storage room and also assume that we have Y 
spares on hand. So that when a machine foils, a spares is immediately 
substituted for it if it happens that all spares are used and a breakdown 
occurs, then the system becomes short. When a machine is repaired, it then 
becomes a spares (unless the system is short in which case the repaired 
machine goes immediately into service). Here we analyze the two cases of 
spares first of them is C>Y and C<Y and also discuss special case. 
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2.1 Case (I) 


In this case we study the maching interference system ; M/M/C/K/ oo with 
balking, reneging and spares Y > C. Then the set of birth death coefficients are 

as follows : 


X 



and 


n // 0 < n < c 

c// + (n-c)a c + l<n<F+k 


In the usual arguments of the 6- techmque the probabiliy differential difference 
equations are : 

P;(t) = - ^Po(t) + ^ = (2.1.1) 

P,:(t)= - (-^l + iw)Pu(t) + ^Pn.i(t)+ l^n<c (2.1.2) 

PJ(t) = - (>9;i + c//)P,(t) +A P^, (t) +(e// + a)P^+, (t) n=c (2.13) 

Pn(t) = - [jai+c/y+(n-c)a]P„ (t) +>92P„^, (t) +[c//+(n+l-c)a]P„^i (t) 

c + l<n<Y (2.1.4) 

Pn (t) = • ^ + 1) (n - c)«] P„ (t) + (n - Y)pl P„.i (t) + 

[c//+(n+l-c)a]P„,,(t) T+l<n<y+K 


(2.1.5) 


As t-oo, the steady - state probability - difference equations are. 
-APo + /iP, = 0 n = 0 

-(A +n//)P„+AP„., +(« + !) = 0 l<n<c 
- (/? X + c// ) Pg + A Pg.j + (c/i + a) Pc+i ~ 0 n — c 


( 2 . 1 . 6 ) 

(2.1.7) 

(2.1.8) 


-[{j3A +cju + in- c)a] P„ + y02P„., + [c// + (n + 1 - c)a] P„^, = 0 

c + l<n<Y (2.1.9) 

-[(n-F+l)y52 + c/^ + (n-c)a]P„ +(«-10^^Pn-i + + ® 

y+l<n<F+K (2.1.10) 

From equation. (2.2.. 6) 

-/lPo + /iPi = 0 n = 0 
//P, = IPo 

jo = A p 
' // ° 

Pi = p?Q Where p = Xj p 
From equation. (2.1.7) 

- (2 + np) Pn + 2Pn.i +{n + l)vWPn+) =0 1 < n < c 

(n + 1) = (2 + n//) P„ - 2P„., 

Putn = l 

2 //P 2 = (2 + /i) Pi - 2 Pq 
2P2=(p+ 1)P, -pPoWherep= 2//y 

2P2=(p+l)pPo-pPo 

2P2=pPo(p+l-l) 


2P2=P^Po 



Le = 


1 

2 \ 
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Put n = 2 

S/iPa = (.^ +2//) P 2 - AP, 
3P3=(p+2)P2-pP, 

= (P+2) j;p'Po-P-PPo 

= (p + 2)^/>.-^2P. 

= ^p’n[p+2-2] 

3P3=|;P’n 


n = 


■if-1 

3! U 3 


Put n = 3 


4/jP,=(Ji +3 /i)P,-AP2 
4 P 4 =(p +3)P3-/?P2 
1 


1 


= (p+3)-p^Po-p-p^Po 


3! 

1 


= (P +3) ^ Po -— Pq 


2! 


= (p+3)jjp^F,-j^3p^P„ 
= Z’’ Po [P + 3-3] 

4P4=^p‘P. 
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Put n = 4 


5 n P5 = (A + Afi) P4 - X P3 
5 P 5 = {/t + 4) P 4 - pPj 

= (P + 4) P„ - p p’ Po 


it 


= (P-4)±p^F,-jy 


p' P« (4) 


4! 


p 4 pJp + 4.4] 




Similarly 

P„ = -^ /7"Po Where /? = — 
n\ n 


n\ 


( oA” 




0 < n < c 


( 2 . 1 . 11 ) 
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Let Y = 1/a and p = c p/a 

From equation. (2.1.8) 

- 1 + Cp) Pc + 1 Pc.i + (cp + «) Pc+i =0 n = C 

-(y^r+^)Po + rPc.i+(^+i)Pc..=0 

(<5 + 1) Pc+, = iPr + Pc - r Pc-i 




c+\ 


_j_ 

" s+i 


^+1 / 

f Pr ^ 

U+iJ 


P + 

" U+i 


1 

y 


^+1 


Pc., = 


Py\ 1 


U+U c! 


-T' 

u> 


1 ( 1 


^“■"U+U c!l,pj 


r 


1 


.c-l 


^^+U(c-1)!Vp 


(Py\ 1 

^^ + lj c! (^p 




, Po 
v/^J 


1 


a(^ + l) (c-1)! Vp, 


1 


U + U c! 


-1 


Po + 




U+U 


1 


1 Ip 


VP. 


a(^+l)(c-l)!c 




Icp 


c! a(<5 + l)l 


- p'= Po + 


U + U c! 


1 


^,<^+1; c! 


^ P' Po - 


5+\ c! 


c+1 


^c+1 


-^1 - Pc 


ly qI 

( P r ] ± 

U + lJ c! 


1 

vp; 
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From equation. (2.1.. 9) 

-{/3X + c// + (n-c)a)Pn + y5AP„., + [c// + (n + l-c)«] P„+i =0 

c-\-\'^n<Y 

-{^Y + <5+ (n-c))P„+ Pr Pn-i+ [^ + n + l-c]P^H.i =0 
[S + n + l-c] P„+, = [/Sy +S +(n-c)]V^- fir ^n-i 


Putn= c + 1 

[^+c+i+i-c]p,«= [fir+s+c+\-c]y^,- fir p, 
( 5 + 2 ) p„„ = (fir +< 5 + 1 ) Po.1 - fir Pc 


(5+2)P.«= (y^K+^+l) 


1 


rx^‘ /'a„\ w 


c! 


f fir 


Po-ySy 


c! 




(<y+2)Pe.2 = 


c+2 


1 

c! 


c! 

^ V 
c! 


(-] 


y^rPo 


py S + \ ^ 

S+l . 


Vo+ly/ 

ifirf „ 


(5 + 1) (5 + 2) 

Put n = c + 2 

[5 + c + 2 + 1-c]P„h. 3 = [y^+5 + c + 2-c]P,,.2- fir^. 


C+1 
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(^+3) p,.3 =(Pr ^s+2) p„3 r P. 

[Pr'f 


c+1 


1 ( 

{0r+6+2)- - 
c! \fi) 


_ p - Bv — 

( 5 + 1 ) (^ + 2 ) ^ c\ 


Kf^) 




c! 


P = 

^c +3 


f -1 

UJ 


ipyf 

5+1 


j3 Y +2 


1 (x^ 


{Py) 


5+2 

,3 


- 1 


c ! v^//J (< 5 + l ) (^■^ 2 ) ( 5 + 3 ) ° 

Put n = c + 3 

(5 + c + 3 + l-c]P„. = [;3r+5 + c + 3 -c]P„ 3 - 38 r P, 

(■5 + 4) P„. = {/? r + <5 + 3) P.,, - ,» / P 0.2 

' (Pr)' P. 

(5+1) (^+2) (5 + 3) 

(Prf „ 


c+2 


{fiy + 5+'i)\ 




c! y^i) 


Pr 


c! 


'a'" 


(5+1) (5+2) 


c! 


-1 




(5+1) (5+2) 


° 5+3 




P = 

^c +4 


c! 




(Pr)' 


(5+1) (5+2) (5+3) (5+4) ° 


Putn = c+4 

(5+c+4 + l-c) P „3 =(>?;' + 5+c+4-c] P.+ 4 - 
(5+5)P„3= (;^ + <y+4)P„,-38rP, 


‘c+3 


= (y 3 ;^ + 5 + 4 ) 


1 

c! 


Pr 


c! 




VP) 


kPj 


(Pr)' 


(5 +1) (5 + 2) (5 + 3) (5 + 4) ° 


(Prf 


(5+1) (5+2) (5+3) 
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fl'l 

(Pr)' p 

'y0y + ^ + 4 J 

lY 

(<? + l)(<^+2)(<J+3) “ 

d+A 


^ 1 taY p 

(^+1) (<5+2)((y+3)(<?+4) (<?h-5) 


Similarly 


P„ = P 


c+(n-c) 


J_ 

c! 




m 




c + 1 <n<Y (2.1.12) 


From equation. (2.L10) 

-[(n-Y + l);0;i +cpi +(n-c)a]P„+(«-10>^^Pn-.+ 

[c// + (n + l-c)a]P„„ =0 Y + l<n<Y + k 

(^+n + l-c)P„,, =[(«-r+l)/?r+^ + n-c]P„- 

(n-7)>9rP„., 


Put n = F + 1 

(^+F + l + l-c) P,^2 =[{Y+l-Y+l)j3r + S+Y+\-c\Vy^, 

-{Y+i-Y)pr?y 


(S + Y + 2 -c) Pj .+2 =[2y?y + (J + Y + l-c] Pj,^i -py P^ 

//^ \y-c+i 



(Prr 

(^ + 0r-c 



1 c p pyg^+^+r+i-c 

c!^ 'L S + Y+l-c 


i^r) 


n-c 

r-c+i 


ipr 


c! ° (< 5 +r+i-c) 


^Y +2 ~ 


{fir)' 


ifir 


c\ (^ + 1)^_, (S + Y+l-c)(S+Y+2-c) ” 


ifir)' 

c\ (<? + l) 


ifir) 


' Y-c+l 
Y-c 


^fir p 

(S + Y+2-c) " 

2 {Prf 


c\ {S+l)y_^ {S+l + Y-c)(S+Y+2-c) ° 

Piitn = F+2 

(^+F+2 + 1-c)P^^3= [(7+2-7+1) j3r^S+Y+2-c]Yy^^- 

{Y+ 2 -Y)PrVy,, 

{S + Y+3-c)Yy,,= [3, Pr + S^Y + 2- c] P,,, -2 py P,,i 

= (3^y+<5 + y + 2-c)4 tM TTTTT^Po 


(^ + }' + 2 -c) 


{Prf 

(< 5 + 0 , 


2 y 8 rP. 


3 ^ + Y + 2 - C 

< 5 ’+r + 2 -c 


Pi !£rYlL 2 Byv ^Pr 

c\ {S+l)y_^^^ fi'KS-^Y+2-c 


7 y-c+i 

, F-+C +1 


Pi (>gr) 23(Brf 

(<5+0r-„, (<5+P+2-c) 
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{Prr^' 2 . 3 . {Pr)‘ 


Y+3 


c! (S+l)y_,^^(S+r+2-c}(S+r+3-c) "" 

Put n = F + 3 

(5+7 + 4 + l-c)P„ = [(r+3-7+1) ^f+iJ + Y + S-cjP 

- (r + 3-Y)/5yP,„ 


r+3 


2.3. (^rT Po 


[A fir +S + Y+3-Z] 


c\ (^ + 1),.^, (<?+F + 2-c) (J + F+3-c) 


3 fir 


p" [Pr) 


Y-cH 


ipr 


(■J + l).,*, (5 + Y + 2-C) » 


^ {/3ry-‘^' 23.(PrY 

c! (5 + l)r-c., {S + Y+2-c) 


(PrY-' 2.HPrr r 

c! (^+l)r_„, (d + Y + 2-c) “ 5 + Y + 3-C 

. _ pT i/^rY"*' 2.3.4.(;gy)^ 

c! (<?+l)y_„, (<y + y+2-c) (<5 + 7+3-c) 


4 fi r + S + Y + 3-C 
S+Y+3-C 

4 fir 


-1 


r+4 


(prY" ^3A.{pry 


cl (#+!),_. (^ + i'+l-c)(<5 + r+2-c)(5+7+3-c) (5+Y + 4-c) 
Similarly, 


n-y 


£1 P. (2.U3) 

c\ {5 + \)y_, (J+F+1-c)_^y 
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Thus P„ can be written as follows. 


— /?” Pq 0<n<c Where p= Xj ju 

n! 


]_ 

c! 


~T 


(firT ° 


c < n < Y 


(firY" (n-'iV-ifirY' p 

c! (5 + 1 ),. (5 + Y + l-c)., ' 


F+lSn<r+k 


( 2 . 1 . 14 ) 


00 

To find Pq the boundary condition, consider ^ P^, = 1 


c-1 1 

It 


+ 




n\ 


c! 


j. 

s 


(PrT 

(<s+i).- 




Pi (yg^r^ y in-Y)\ {^rY'" 1 


(A) (2.1.15) 


We now consider 


(i)Z 


iM 

( 5 + 1 ) 


n — > n + c 


h{^+ 1 )„ n! 


y iJ^rY (i)n 

h (g + l)„ n! 

= ,F, (l ;5 + l;) 5 r) 


£l f {n-Y)\ {firT' 

c\ (5 + l),_, „4(g + J^+l-c)„_, 


pT {firY- + [PyT 

c\ ( 5 + 1 ),.^ ^(5 + 7 + 1 -c)^, 


n-^n + r + l 


P° (ygy)^'^ y {n+\)n\fir [PrY 

c ! (5 + 1),_^ h {5+ Y + 1 - c)(5 + Y + 2 - c)„ 
p 1 ” («+y)n! {pr)’’ 

c\ + h (^+Y+ 2 -c)„ 
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\y-c+i 


pL Y— 

(^+1).-.., (<^+Y+2-c)„ 


{Py) 


Y~c+\ 


z 


(n + 1)! {firj 
(<5+Y+2-c)„ 

(2)„ (>gr)"n! 


(<5+0y_c+i + n! 

£l iPrT"' y ( 4 (^r)”(i)n 

(^+1W, (^+Y + 2-c)„ n! 

£l i^Z) F, (l,2;<5 + F + 2-c;Ar) 

c! 

Put in (A) we get. 

p;' “i;4+4 


xy-c+1 


BL Mil 2^1 ( 42 ;^+Y + 2-c;/0r) (2.1.16) 

(^ + 1)k-c., 

Where (1)„ = 0 where n > Y - c in the first hypergemetric function. 

(n+1)! = (2)„ 

To calculated the expected number of units in the system, a result due to Abou 
- El - Ata is used, which states that for a simple birth death process. 

L= ;iPn ^ 


dX 


P 

^0 


XK 


Y^+^,F,{l;S+Ufir) 

n=0 


dX 

V y-c+i 


+ 


,F,(l,2;5 + y+2-c;Ar) 


c! (<y+i) 

d 


(D2P. 


dX 


2^ 

lA^ 


-c+1 


/ \ 

n ” 

c4 1 1 

zZ 




\P) 





Y-i- 

La 

n=0 


; «! ft 


C-l 1 


n=0 


(n-l)! /i"J 


P y_J_ i: 

U{n-\)\ /i" 


p y — L_ o-" Where o- = X j f 2 


(ii) A Po 




B- ,F, (l;^+l;y^r) 

Cl 


APo,F, (l;^ + l;y^r) 


dX 


cl 


u> 


+ 1 P„ 




c! dX 


[,F, {1;^ + 1;,5 r)] 


AP„ ,F, {U<5 + l;/J 


E 


(1). {PrT 


S (•5 + 1 ), ■>! 


,F, {l;S + Ufi r) 


(c-1)! 


lY - ( 1 ). A" '^’Z- 

\M) n =0 


,F, ( 1 ;^ + I;y5 r) 


a 


X"J 3 ^/ 

« ” (l)n” 


(0-1)!'’“'"''” c! § (^ + l),n! 




P’ vJlUMl 


(c-1)! “ ® cl (n-1)! 


a” ^ ^ p‘ n! {0/)" 


]Fi (l;^+l;^ y) Po + Fo . X! 


(c-1)! ° ‘'cl ^(<5+l)„(n-l)! 




(c-1)! c! ^(^+1), 


,F,(i;^+i;,fi i')7:^P.+Po^ 


(c-1)! 


cl 


nj3r iPr) 


n-l 


^ {S+l){S + 2)„_^ 


,F, {1;S + UP r) 7—777 Po +Po^Z 


^ pYXi. \pr) 


n^\ 


(c-1)! “ c! (^ + l)(«5 + 2)„_, 


,F,(l;^ + l;y0 y) 


o-° p p pIy /?rn {PrY~ nl 
(c-1)! ° c!^ (<^ + l)((J-h2)„_, n! 
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,F, r) 

,F, {US + l;fi r) 
,F, (US + U/3 r) 


(c-I)! 


(c-1)! 


Po + Po 


c! 


I 




{PrT 

((5 + 1) (^+2)^_, (n- 

■1)! 

[PrT (l)n+l 


((5 + 1) (^+2)„ n! ^ 


Pr (Pr)" (2). 


(5 + l)(,5 + 2)„ n! 



.F, {us+Ufir) P.+P. iP. (2;5+2;^ r) 


umxp. 


A- 

dX 


p 


[Pr) 


r-c+1 


C! 


,F, (l,2;^ + Y + 2-c;/?r) 


C+I 


APo 2F1 (l;2;^ + r + 2-c;y9r) 






p' (^r) 


y-c+i 


;! (^ + 1) 


y-c+i 


+ 




{Pr) 


y-c+i 


d 


C\ (^ + 1)..„, 


[ jF, (l,2;5 + r + 2-c;/?r)] (2.1.17) 


(iii)(a) APq 


dX 


P° {Py) 


y-c+1 


:! (<5 + l) 


y-c+1 


APn 




A!:::! 

(<^ + l)y-c+i 


IPn 


1 ; 
i - y-c+i ^ 


ipY 


P 


a 


y-c+i 


dX c!(J+l) 


y-c+i 


2P„ 


J 

P 


^ (7 + 1) 


y-c+1 


a 


c! {5 + 1) 


y-c+i 
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^r-c+i (7+1) A 


Y+\ 


a 


7-c+l 


c! (^ + 1) 


F-c+1 


(y+l) 


>=! (' 5 +')>'-..l f‘' 

(J'+O ,< 


-P 


a 


7+l-c 


c! {S + \) 

(r+i) 


7-c+l 



y+i-c 




P‘ (/ 9 y) 


7-c+l 


c! (5 + 1),, 


K-C + l 


,F,(l,2;^ + Y + 2-c;yffr)>lPo 




(^/) 


r-c+1 


d (^+ 1 ) 


,F,( 1 . 2 ; 5 +Y + 2 -c; 5 r)(l'+l)^^ 


y-c+1 
r-c+i 


c+l 


{«>) (« /IP. 4- [ .F, (l,2;5+Y+2-c;5r) ] 

m '• 


= APo — 

® M 


,<1).(2). 5" 


n=0 


(^ + 7 + 2-c)„ n! 


^0 




« (l)„(2),n A"-‘ 

S n!(J+Y + 2-c)„ 


« (l)„(2)„ n^}/„ A" 
^ n!{«5 + Y + 2-cX 


- (l)„(2)„n(^r)'’ 

^ n!{^ + F+2-c)„ 


f (l)n ( 2 )n (ygx)” 

(n-1)! ((J + F + 2-c)^ 
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00 

p.i; 


1 




^ (a + Y + 2-c) 
1 


(l)n (2)n {PrT 
(n-1)! {S + Y + 3-c)^ 

( 2 )„.. ( 3 )„., {fiYT 


.\.2.j3r 


^ (^+r + 2-c)' • (n-l)!(<5 + r + 3-c)„_, 

1 2/3r 


{S + Y + 2-c) {S + Y + 3-c)^ 


, r jF, (2,3;<5 + 7 + 3-c;^y) 

{S + Y + 2-c) ^ ^ 


Put in(B) ( 2 . 1 . 17 ) 


(y + l) P. ,F, (l, 2 ;^+Y + 2 -c;ySr) 


+2 Pr 


iPr) 


Y-c+\ 




P„ ,F, ( 2 , 3;<5 + P + 2 -c;ySr) 


c+2 


r/zM 5 L = /I P, 




0 oo 0 


= P„ 


c-1 n 


^^,F, (2;..2;^,). 


(r+l)p- (^y) 
c! (S+l) 


Y-c+l 


,F,{l,2;^+¥+2-c;jgr) + 


Y-c+2 


2Pr {py ) 

((5+1) 


Y-c+l 


jF, (2,3;<5+r+2-c;>0r) (2.1.18) 


Y-c+2 
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22 SPECIAL CASE 

PQ'(t) = -2Po(0 + /^Pi (t) n = 0 (2.2.1) 

P' (t) = - (2 + n/i)P„ (t) + 2P„., (t)+(n + l)/iP„,,(0 

\<n<c (2.2.2) 

p; (t) = - 2 + c/i] P^ (t) + A P,., (t ) + (c/i + a) (0 

« = c (2.2.3) 

P;(t) = -[^A +c^ + (n-c)a]P„{t)+y3A P^,(()+ [c;/+(n+l-c)a] P„,(0 

c+l<n<r (2.2.4) 

P' (t) = _[(„_r + i)yg2 +c// + («-c)a]P„ (t)+ (n-Y)pX P„.,(0 + 

[c/i + (« + 1 - c)6i']P„+, (0 F + l<n<F+i^ (2.2.5) 

Let a = 0, = 1 then above equations are: 

P'(t) = -;iPo(0 + /^Pi(t) 

P ' (t) = - {A + n/i) P„ (t) + A P„., (t) + (n + l)/i P„+, (0 1 < « < c (2.2.7) 

Pc'(t) = -[^ +c/i]P,(t) + 2P,.,(0 + c/iP,^,(r) n = c (2.2.8) 

Pn (t)='- ( ^ +c/i)P„ (t)+ 2, ¥^_T^(t) + c^P„+i(t) c + l<n< Y (2.2.9) 
PMt) = -[(n-r + l);i +c/i]P„(t) +(n-Y) A P„.,(0 + c/i 

P„^,(t), P+l^n<r + i^ (2.2.10) 

as t -^00 then the steady state probability differaice equations are 

-2Po +/iP, =0 (2.2.11) 

- (2 + n/i) Pn (t) + 2P„.j + (n + l)/iP„+, = 0 (2.2.12) 

-[2+c/i]Pn +^Pc4 +c/iPc+i =0 (2.2.13) 

-XI +c/i)P„ + 2 P„4+^/^Pn.i =0 (2-2.14) 

-[(«-r+l)2 +c/i]Pn +(n-Y) 1 P„.,+ c/i P„,, =0 (2215) 
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Fiom equation (2.2. i 1) 

- AP„ + /tP, = 0 

/tP, = AP„ 
P,=-Po 
Pi = PPO 


From equation (2.2.12) 

- (a + nfj) Pn + Pn-l + (« + l) Pn+l = ^ 

(n + l) // P„., = (d + n/t) P„ -aP., 


Put n = 1 in above equation 

2 // P 2 =(>l + /^)Pi-^Po 

2 p,=(p+i)p,-pp„ 

2 Pj = (/0+l)pPo-/’Po 
2 P,=pP„(p+l-l) 



Put n = 2 in above equation 

3 // P 3 = (A + 2 //) Pj — /I Pj 
3P3=(p+ 2) ^ Pj-pP, 


l<n<c 

•1 
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= f^Po[P+2-2] 

P3=^p'Po 

Similarly, 

P„=— 0<n<c (2.2.16) 
n\ 

From equation (2.2.13) 

- ( 2 +c//)P<. +2P<..i +cfi P,+, =0 


Pc+i = (^ + Pc - ^ Pc-i 

cPcH-. =(p + ^)Pc-pPc-i where /? = ^ 


CP„,=(P+-)^P‘ 




= ^p'‘ Po [p+c-c] 
c! 



From equation (2.2.14) 

- [ /I + c//] P„ + 1 P„., + c// P„+, = 0 
C/iP„, = (2 + c//)P.-AP.., 
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l+c^n<Y 


cP„,,=(p+c)P„-pP„.i 

bove equation 


CP,.2=(P + ^)PC.1-PPC 



Put n = c+2 in above equation 

C Pj.+3 = (p+ c) Pc+2 ~ pPc+1 



Similarly 




c<n<Y (2.2.17) 
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From equation. (2.2.15) 

_[(„ - 7 + 1) A +cfi]Y,+{n-Y)X P„_, + cfi P„,, = 0 
c Pn+i = [(” - + 1)P + c] P„ - (n - P)pPn-i 
Put n = F + 1 


cP,„ =[(r+i-i'+i)p+<;]PM -(P+i-iOpPv 
cP,.,= [2p+c]P Y+\ P?y 



Putn = F+2 

c P),^3 = [( 7 + 2 - r + 1 ) p + c] Pk^2 -{ Y + 2 - Y ) p?y +2., 
cP,^3= (3/>+c)P,,2-2pPr+i 


cP,,3= (3p+c)^p' 

^ vK+l-c 

2^ip' f P. 

c! / 

cfr.,= ( 3 p+c)^/’ 
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^ c! ^ 


, r-c / \2 

P 


^1 2 ^ Po(3p+c-c) 

vc; 


Py+3 ~ , P 


,r-c 

c 1 P I 01 I 1 p 

•B* A 


c! 

1 


231^, .0 


/ xK-c / \3 

. -_p“f£^ 2.31^ 

c'. IcJ U; 


Similar/y. 

£ ~ Py+cn-y) 


= 

c. U; 


/ p 




(n-Y)\ 


c 

,n-Y 


P„ = -^ 4rr("-Y)! 1^1 Po (2.2.18) 
c! c 


Thus P can be written as follows ;- 


£ = 


n! 


P^Po 


0< n<c 


ipnjpj c<n<7 (2.2.19) 

c! \c) 


,n-Y 


^ P I 1 Y + l<n<r + k 

c! 


Y-c 


To find Po the boundary condition P„ - 1 


If; 


w =0 


c P^ 

+ ^ 


A 

n ==0 


£1 y|£ 

c! 


c! c 


r-c 


I(n-Y) 


n=Y+l 


V 

VC. 


,n-r 


p - = y 1 p. + £l " +-^ t (n-y)! F az-i^) 

“ tsar c\ tiU) c!c a;, 
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(0- ^ Z 7 


^y(pX^ 

c! ^UJ 


n! 


c y-c (^)n , , 

p ^ vc; 


c! 


n-0 


nl 


^ ,Fo(l;-;p/c) 
c! 




,/i-K 


c! c 


- Zf"-')! 


n=0 


F * (2)„ K n! 

P y P \^J 

Cl c"-^ c 


n! 


^ Y+\ 


P ■■ F 1 2-^ 


Hence 


C-I ^ n 


P-‘. = + ^ ,F. (l«Wo)+;f^ A [ 1.2 k 

n—0 


y+1 




n! c! ' c! c' 


andL = lPo — Po' 


I = 1P„ 


dk 


c-1 fS 

Z^+^ .Fo 1 1.2;-; 

U n! cl 


c! c 


c-l 1 f 2 

omp.zA - 

Un'- KP) 


f P” 

S(n-l)’- 


(2.2.20) 
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IPn 


dk 


4lFo(lr;pM 

c\ 


3 p" . -.v, p° A. 




cl M 


ipQ 


,, , , ;i‘-' p' ^ 

z?, ,F„ (i;-;Wc)<= -TTT+^Po — L 


//"c! 


c\ ^ S «'- 


C 00 

Po 7"^ iFo (i;-;/>/e) +.iPo £(1)« 
(c-1)! " 


n 


j 


n=0 


n\ 




(c-l)! 


(c-1)! 


,Fo (l;sp/c) + Po Z 


cl ^ (n-1)! 

p' (^)n-i (p/^r* 


^c! (n-1)! 

P F (I- -p/c) + P ^ Py^^hiEK. 

° (c-1)! ° cl n! 




(z'h). APq 


dk 


. y+i 


c! c 


r-c+i 2 


V^{\,2-,-\pIc) 


AP„ ,\(\,2\-,pIc) 


r+1 \ 


^3l\c! c 


y-c+i 


+ 


;iR 




° c! m 


[ ,F,(l,2;-;p/c)] 


AP„ 4 (p ’■*')= 

«i P 

7i 

= (p+i)— Po 
// ' 


(r+Op'^'Pc 
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/IPo 


Ck 


(2)„(1)„ (p/c)" 


n=0 


n\ 




S 

n=0 


C JJ- 
nl 


s 

n=0 


(2)„ (1), jp/c)" 
(n-1)! 


^ ^ {p/^) ^0 


n=l 


(2)„-i (3)„., (p/cf 
(n-1)! 


I.2(p/c)P,|; 1.2(p/c)P, ,F, 

n=0 


2 3'-- — 




-1 ^ rr A 

—n /J 

■ + 




,F,(l;-;p/c) + 


( f r i)p' 


i + % 




- F- (1 


C + 1 


cl c 
(r+i)p’'*‘ 


Fo (^<-,pIo) + 


cl C 


r-c+i 2 


Fq (1,2;-;p/c) + 


,y+2 


c!c 


y-c+2 2 


Fo (2,3;-; p/c) 


.,2; -. ; p/c) +, 


(2.2.22) 


2.3 Case (II) 

This case treats the system of Machine interference M/M/C/Ky oo with 
balking reneging and spares y < c and hence the birth - death coefficient 
are 


';L 0<n<Y 

(«-7 + l) Y<n<c 

(«-r + l) pX c<n<r + k 

0 n>0 


and 


XifJ. 

c// +(n-c)a 


0<n<c 

c + l^n<F + k 




Then, as before, the probability differential - difference equations 


are 

Pat) = -APo(t) + /iP.(t) 

P;;(t)=-(^+iw)P„(t)+>tP^,(t) + (t) l<n<p 

p; (t)={ («-r +1)1 +iVJ: 1 P„ (t)+(n'T)lP„-i (t) + (n+b/i P„+i (h 

y+l<n<c 


(23.1) 

(233) 

(2.3.3) 


pj (t) = .[ (c-r +1)^ +c//]P, (t)+(c-101P^, (t) + (c//+a) P^, (t) n=c 

Y+l<n<c (33.4) 

p^(t)=.[ («-7+l)y8t +c;^+(«-c)a]P„ (t)+(n-Y);fflP^, (t) +(c//+(n+l-c)a) P^,(t) c+l^n<Y+K 

(23.5) 
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There fore the steady state probability - difference equation are 
-APo + ^iP, =0 

.(A + n/i)P„ + ;iP„.,+(n + l)/^P.*, =0 l<n<Y (2.3.7) 

- [ (« - Y + 1) 1 + n/i] P„ +(« - P)^Pn-l + (” + l)/“Pn+l = 0 

y + l<n<c (2.3.8) 

-[{c-Y +\)PX +C;u]P.+(c-T)2P..,+(c/i + o:)P.., =0 

n = c (2.3.9) 


- [ (n - Y + 1) ;®l + + (« - c)a] P„ + (n - P„.i + 

[c// + (« + 1 - e)a] P„., = 0 c+\<n<Y+k (23.10) 

From equation (2.3.6) 

- A Pq + /i P) — 0 

//P, = APq 

p, =-p. 

Pi=pPo 

From equation (2.3.7) 

- (2. + n/j)P„ + AP„., +(« + l)^<P„., =0 \<n^Y 

(n + 1)// P„+, = (2. + n//) P„ - 2. P„., 

(n + l)P,., = (/7+n)P„-pP„., 

Putn = r 

(r + l)P„ =(p+P)Pr-pPr-i 
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Put n = 1 


2 P 2 =(p+l)^i~pPo 




2 P 2 = pPq (/?+1-1) 






Put n = 2 

3P3=(p + 2)P2-pP, 

= (p + 2)^n-f(/’P.)*2 

P ,=^^ ob + 2 - 2 ] 



Put n = 3 


4 P 4 = {p+^)P^~ P^ 2 


_2 \ 


= (/>+3)^Po-f 


v 2 ! 'y 


.3 


^Po[p+3-3] 



♦ 

m 


Similary. 


p=— 0<n<C (2.3.11) 

" n! ^ ° 

From equation (2.3.8) 

-[(n-Y + l)l +n;/]P„+(n-50AP„., +(« + l)/iP„.i =0 
(n + 1)// P„, =[(n-Y + l)A + n^u] P„ - ^ Pn-i 


Putn = F+l 


(r + 2)// P,+2 = [(r + 1 - F + iM + (F + l)jo] Pr+I - (F + 1 - F);iPy 


(y+2)p,-« = 


2^+(Y+\)i!- 


p -Ap 

^r+1 


(F+ 2) P,« = [2/> + r + 1] Pk^i - pP)- 


>Y+i 


(F+ 2) P,,, = (2 p + 7 + 1) n - P 7T ^ ^0 


Y! 

Y 

p 

Y! 

_ F 


2p+7 + l ^ 
7+1 

' 2p 

T+1 

2p^ 


y+2 


71 (7+1) (F+2) 
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Putn = r + 2 


(f+3);uP,, 3 = [(7+2-r+i);i +(r+ 2 )//]p ,,3 - 

(r+2-yUp 

\ / r+2-1 

(r + 3) p,,, = [3 /? + (y + 2)] p,« - 2 p p,„ 


(r+3)p,,3 = [3/7+(r+2); 


P 


Ip 


Y\ (7+1) (F + 2) 


P. 


2p 


p Ip 


Y! 7 + 1 


p'^ 23 p^ 


7! Y + 1 

3 p 
7+2 


P, 


r+3 


Y! (7+1) (7 + 2) (7 + 3) ® 


Similarly, 

K = Py.(„-X) = y < n < c - 1 (2.3. 12) 
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From equation (2.3.9) 

-[(c - y + l)/3/l + c/i] P, + (c - y)/lPe., + (cfi + a)Pc 

=[(c-Y + \)Pr +^]P.-(c-i')rP..i 

a 


(s + 1) p„, = [{c - y + 1)/S r + ^] p. - (c - P)r p»-. 


5 + 1 ^ ' <5 + 1 


<5 p (c-P)v p 


(c-V)/ (^F)r (c-l-Y)!p p 

(^ + 1) (S + l) F!(F + 1 ),.,.y 

p ^ (c-Y)(c-l-Y)lp"*'^ -c.pr 
Y! c(7 + 1)c.i-y P 


p^ (c-Y)!p--^ p _cr__ 
Y! (r+i)c.Y ° p(^+i) 

c l/a 
X/jj,(5+V) 


s 

S + l 


Put in (A) we get. 


(c-F + 1 ) , „ S ^ d ^ 

=^~r-T-^ Pr Pc +— P- — P-- 


S+l 


S+l s+l 


(c-Y +1) ^ ^ 
p L pyP 

c+1 ^ J 


= 0 n = c 


(A) 


From equation (2.3.10) 

-[(« -Y+ 1)J3X + cju + {n- c)a] P„ +{n- 7);^ A P^., + 

[c/i + (n + 1 - c)a] Pn+i = 0 c + l<n<F+k 
\cfj. + (« + !- c)a] P„+, = [(« - 7+ VjJ^A + c^ + (n- c)a\ P^ - (« - Pn-i 

[<5 + H + 1 - c] P„, = [(n - y + \)l3y + <? + («- c)] P„ - 

(n-y);erP„., 


Put n = c + 1 

[,J + c + l + l-c]P,«=[(c + l-y+l)/?r+<y+c + l-e]Pc«- 

(c+i-y)^y Po+i-i 

(s + 2) p,*j = [(c - y + 2)j3r + <5 + 1] P.,, - (c - y + 1)^ p„ 

= [{c-Y+2)i3r+s+\] p,-(c-y+i);gr Pc 


(c-r+i);^ p. 


(c-Y+\)Pr^. 


(c-r+2)>9/ +5 + 1 


5 + 1 


-1 


{c-Y + 2)pr 

(5+1) 


^c+2 


(c-r+l)(c-y+2)(;^)- 
(5 + 1) (5+2) 


(c-y+l)(c-y+2) p’'(c-Y)!p 


c-r 


c+2 


(^ + l)(<^ + 2) 


Y! (Y + !),_, 
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Put n = c + 2 

(<5 + 3)P„3 =[(c + 2 -r + l)j8r +<?+c + 2 -c]P „2 

-(c + 2-Y)fiy P„, 

= “ P + 3)j0 7 +2] P^.+2 " ^ + 3] ^ Pc+i 

3 {c-Y+\)(c-Y + 2) (Pr)^ P 
= [(c-Y+i)l}y +S+2]^ (5+1) (5+2) ° 


, (c-K+1),^ _ 

(c-Y^2)Pr P. 

(c-y+l)(c-P + 2)(^y)" „ 
' (^+0 


(c-Y + '^Py + ^ _ j 

<5 + 2 ^ 


(c-r+i)(c-p+2)(^rf p (c-r+3)^r 
(^71) ' (^+2) 

(c-P+l) (c-P + 2) (c-P + 3) (^r)" n 
(<5+l)((5 + 2)(5+3) 

(c-P+l)(c-P + 2) (c-y + 3)(>»/ p^Cc-Y)! p"-'" 


P.*3 = 


(5 + 1) (5 + 2) (5 + 3) 


Y! (P + 1) 


c-r 


similarly, 

P — P 

^ n ^ c+{n-c) 


/ (c-Y)! p*'’ (c-Y + !)„., (fir ) ^P + k (2.3.13) 


r!(P+l)„_ 
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The empty system probability Pq , can be found from the 

00 

boundary condition: ^Pn = 1 

n=0 


«=0 


z 


y! ^ (J'+IU 

p’'(c-Y)!p‘-^ (c-Y + l)„, {Pr)' 
y!(7+l)^_, (<5 + lU 


(/) 


^ [n-Y)\p 


n-Y 


z 


y'- iS? (J'+O,-. 

p ’’ n! p" n! 

pL ‘-f (I)n(l)„P" 

n u (l'+l)„n! 

^ ,F,(1,1;Y + I;;5) 


n — > n + 7 




p''(c-Y)'.p'-'^ (c-y-^\\.Ai}r) 


n-c 


^ r!(p+i) 




(I'+iL. 


P‘i<,-Y)\ ^ (c-y + l),(^y)" 


Y\ (r+ 1 ) 


c-r n=0 


(<5+0, 


n -> n + c 


p‘ (c-Y)! ^ (c-r+i), (Pr)" (1)„ 


y'- S 


(J'+l). n! 


p^ (c-Y)! 

F! (7 + 1) 


,F, (l,c-Y + l;^+l;>?;') 


C-Y 


(2.3.14) 
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,F, (i,i;Y+i;/^) + 



^ 1 


n^n! 



(c-Y) !_ (1 c-Y + l;<? + l;/??^) (2.3.15) 


JL = AP. 


IP. 


A 

dl 

_d_ 

dX 


t-, P"+^.F,(l>l;Y + l;p) + 


p‘ (c-Y)! 

J'! (l'+l).-v 


^V,{\,c-Y + \\S + \-,Py) 


® 1 si; 


;tP Z yl 

^ ^,n\ 


n-1 


^ ^ 1 X 


r-i n 

p« E ^ 


■ir=0 


(n-l)! 


(//)APo 




Yl 


,F,(1,1;Y + 1;p) 


;ip. 


cX 


Y ! 


,F, (1,1;F + 1;^) + AP„ 'fy J; [ 2 F. (l,l;Y+l;p)] 
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r rr^n2¥ — (l ,l ; Y + 1 ;p) 

{II)a X¥o Y\ 






(i,i;Y + i;p) 


Po 


r ^ 
P 


,F,(i.i;Y + i;p) 


.(F-l)!. 


APo 


p 

Y! 


n-l 


„ (l)n 


W.P" 

F!S(Y + ')>-*''■ 


n-l 


_P E-Y 
- ^0 V I 


^ 1. 1 ^(2)..,(2)n., P __ 

n^(r + l)(F + 2)„., (n-l)! 


= P. 


p’'"' f ( 2 ). ( 2 )nP" 
(F+l)!^ (F + 2)„nl 


y +1 


= P„-^ ,F,(2,2;Y + 2;p) 


(F+1) 


{llI)X¥o 


<31 


JSlDL. ,R (l,c-Y+l;<5+W) 

F! (F+ILf 


= jF,(1,c-Y + 1;5 + 1;M^Po 


<31 


pMc-Y)! 


F! (F + 1) 


c-r 


£l iSlYl!- A?, — [ ,F, (l,c-Y+l;<J + l;#)] 


F! (F+1) 


<31 


(///)aAPo 


<31 





p'(c-Y)! 

11 

c;i-* (c-Y)l 

_F! (F+1),., _ 

PI (F + l)e-f 



F! 


= Pn 


p' (c-Y)l 


F! (F+1) 


c-K J 


61 


(IIl)b APo ~ 2 


“ 2F1 (i’C-Y+i;^+i;^r) 
dX 


A Pfi 


dX 


s 

n=0 


(1) . (c-v+i), ^py)' 

'(■5 + 1), 


C3C 

Pc I 


, ^x 

(l)„(c-Y + l)„ /«■ 


n=0 


(^+1). 


n! 


,a(1).(c-y+i), (/irY 

(5+1), (n-1)'- 


n=0 


^\{o-Y + \)Pr ( 2)-, (»-Y+2)-, (^5') 

S (5+^) 

^^(c-Y + 1) ^ (2). (g-Y+2), (^y)' 


(5 + 1) 


z- 

n=0 


(5 + 2)„ n! 


{s 

+1) ’ 

(ID) 


IP ^ 
di 

(c- 

7! (7 + 


. 5' (c-Y)! „ 

,F, (l.c-Y + l;5 + l;^l')c 77 P» 


+ 


p’ (c-ni ^r(c-Y + l) F,(2,c-Y + l;5 + 2;;0y)P„ 

Y! (r+i)_, ^■'•1 


Hence. 


I = Pr. 


y — + 7 ^ ,F, (l,l;Y + l;p) + 
i (n-1)! (r-1) 


F (2^-Y + 2;p)+ ,V,{Uc-Y + \-,S + l;Py) 

^^P}. YI(r+l)e-Y 


P (2 Y + 2;^ + 2;y^)] (2.3.16) 

F!(F+l)c-r (‘5+1) 
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2.4 SPCEIALCASE 


Leta = Oand/?=lthenM/M/C/K/cowith spares only 
we have. 

P'(t) = aP„(t) + //P, (t) n = 0 (2.4.1) 

P: (t) = -{X+ nju) P„ (t) + ;i P„., (t) + (n + 1)A P„., (t) 

1 < n < r (2.4.2) 

P^t) = - [ (n - Y + l)/l + n//] P„ (t) + (n - Y)1 P„., (t) + 

(n + l)//P„+, (t) r+l<«<c (2.4.3) 

Pc'(t) = -[(c-E +m +c//]P, (t) + (c-Y)/lP^, (t) + 

0 fi + Pj+i (t) n = c (2.4.4) 

P; (t) = -[(n-Y + l)/l +c//]P„ (t) + (n-Y);iP„., (t) 

+ c // P„+j (t) c + l<n<Y + k (2.4.5) 

As t - 00 then the steady state probability difference equation are. 

+//P, =0 n = 0 (2.4.6) 

-(2+n/i)P„ +2P„., +(n + l)^P,,, =0 

l<n<Y (2.4.7) 

-[(ii-Y + l)/l +n//]P„ +(n-Y);iP„.i + 

(n + l)//P„^, =0 Y + 1<7i <c (2.4.8) 
-[(c-y+l)^+c/i]P, +(c-Y)4,P^, + 

c/rPj., =0 n = c (2.4.9) 

-[(n-Y + l)2.+c//]P„ +(n-Y)aP„., + 

c/iP„^,=0 c + l<n<y+k (2.4.10) 
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From equation (2.4.6) 

- APq + //P, = 0 

/^P,=^Po 

p,=%p. 

P] ~ P Pq 

From equation (2.4.7) 

- (yl + n^)P„ + X P„., + (« + 1)// P„+, = 0 1 < n < Y 

(n + 1)// P„^, = (A + n//)P„ - A P„., 

Put n = 1 

2;UP, =(1+//)P,-1P„ 

2 P, =(p + l)P,-^P„ 

= (/»+l)pp„-pp„ 

2P2 = pp„(p+l-l) 



Put n = 2 

3>ju P 3 = (A +2//)P2-AP, 

= (pP„) 

3P3=(p+2)^i>_^P,2 

= ^P,(p+2-2) 
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Similarly, 

Pn=A/’"-f’o. 0<n<F-l (2.4.11) 

n\ 


From equation (2.4.8) 

-[(« - y + 1)4 + n|j] P. + (n - 104 P,., + (« + 1)/^ P„*, = 0 


Y + \<n<c 

(n + l)//P.*, =[(n-r+l)4+n/;]P„ - (n-r)4P„, 

Put n = r+1 

(y + 2 )// p,„ = [(r + 1 - y + 1)4 + (p + i)/i)] p,,, - (y + 1 - y)4 p, 
(y + 2)// p,,2 = [2 4 + (y + 1)] p„ - 4 p. 


(y+ 2 ) [2p + (Y+\)] 

^r+i pVy 


= [2p+(y+i)] 


-E p -n—E 

(r + 1)! ° (Y + 1)Y! 


Po(l^+l) 


Y+\ 


(y+1)! 

y+\ 


p„[2p+(y+i) -(y+i)] 


Po.2p 


(y+1)! 

P 2p 


■ y+2 


Py+2=' 


(7+1)! (7 + 2) ° 


71 (7 + 1) (7+2) 
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Put n = 7 + 2 


(r+ 3 )/^p,„ = [(r+2-r+i);i+(y +2)/i]Pr*2-(r+2-y)/iP,,, 

{F + 3)//P,„ = [3^+{Y+2)M]?y^^-2XPy„ 

(Y + 3) P„ = [3^+(r + 2)]P,«-2/)P,„ 


[3p+(7+2)] 


2p- 


Y+l 


ri (7+1) (7+2) 


Po- 2 p 


(7+1) 


[3p+(7+2)' 


2p‘ 


O ^7 

p _ ^ P P 

7! (7+1) (7+2) ® 7! (7+1) 


Ip 


7! 7+1 


3p+(7+2) 


7 + 2 


-1 


Ip^ 3p 


r+3 


7! 7+1 7+2 

P" 2.3 p^ p 

7! (7+1) (7+2) (7+3) ' 

Similarly, 

_p^(n-Y)!p"-^ ^ 


P„ = P 


r+(n-y) 


P! (P+lL 


Y<n<c-l (2.4.12) 


From equation (2.4.9) 

- [{c-Y+\)l+cp] P, + (c- P) +c//P,^, = 0 
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cn p„, = [(c - r + i);i + c//] p, - (c - F) ^ Pc-i 
c P,„ = [(c - 7 + \)p + c] P„ - (c - F) />P«_, 

-[(c F + l)p c] /Y + 1')_„ “ 


(c-F)p 


F! (Y + l),.v 
cp‘'(c-l-Y)!p' 




cY! (Y + l),.,.v 


_ p' (c-Y)!p‘ *'.p ^.V 4 .n^+r-c] 


F!(Y + 1),.y 


F! (Y + l),.y 


c-K+1 


p Mo-F + l)!p-"' r p‘'(c-F+l)(c-Y)! p , 

Po« p, c(Y + l),, “ 

From equation (2.4.10) 


r!e(Y + l),.^ 


-[{n~Y + l)\A + c^] P„ + (n- Y) X P„_, + c/i P„+i - 0 
c/jPrt = [(«-F + l);i+cp]P„-(«-F) 
c P„^, = [(« - F + l)p + c] P„ - (n - F) pP„_, 
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Put n = c + 1 


cPc* 2 = [(c + l-F + l)p + c]p„,-(c + l-}0 

= [(c-F + 2)/2 + c]P„,-(c + 1-50! /JPc 


P J F!c(F + l), 


fc~Y 








cP 


ric (r+iX_, 

p'’ (c-Y + \)\ p”-^" 


c+2 


Y! 


(^+ 1 U 


P,(c-F 4 - 2 )p 


P _ pHc-Y+\)\ {c-Y + 2)p 

Y\e (F+l)^_^ 


Pc.2- 


(c-Y + 2)1 p ^ _p^ (C-Y+ 2) (c- Y+1) (c-Y)lp^ -^^^ 


F!c^(F+lX_, 


Po = 


Ylc^ (F-MU 


(c-Y+l)(c-Y + 2) A^p[ (c-Y)!p‘^-^ 


(cpY 

(c-F+1) (c-F + 2) 


F! 


(F + 1) 


c-F 


{cpf 


Put 11 = c + 2 

c = [(c + 2 - F + 1) p + c ] P,,, - (c + 2 - F) pP. 


C+I 


^ ^ • Y! c' (^+Oc-F 

(._y,2) £ £:(£z21ti)l^eP. 

c, F!c(F+l) 
p" (c-F + 2)! 


c~F 


F!c^ (r + 1) 


P„[((;-r+3)p+c-c] 


c-Y 
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cP. 


c +3 


ric^ (p+iU 

p’’ (c-r + 3)! 

nc^ (r+iL, 


_p' (c-7+3)!p'-’'*’ (c-y+3)(c-r + 2)(c-r + l)l' 

r!c> (r+i),_, <■ (cpY 

_p’' (c-r + 3)(c-7+2)(c-r+l)(c-Y)!p^„ 


Y! 


<:’ {p+>L, 


Similarly, 


p = p 

^ n c+(«-c) 


^ (c-7 + n-c)! p 

ric"'^ (7 + 1) 


c-7+n-c 


c-y 


/?"(«- 7)! f 


" 7! (7 + l)^_, UJ 


P„ = 


(c-7+l)„_, (c-Y)! p 


c-Y+ri’-c 


7! 


0"-“= (7+1) 


c-y 


(c-r+l).„(c-Y)! 

P! (P+lU 

p' (c-n: p*-’' (c-y+i):^, 




M" 


(2.4.13) 


The empty system probability Pq , can be found from the 

PC' 

boundary condition; y",Pn = 1 


n=0 




Y~\ .n c-1 


tS «! Y! U (K + l),, 

00 






Y! (Y+lL, 


P 

Vc, 






Y! (F + l)„_ 


" nip" 


n -> n + Y 


I; 


Y! ;S(Y + 1). 

Y c-i-Y p " n! 


z 




P 


F! 


^F, (l,l;Y + l;p) 


nn v P^(c-Y + l)n.c (c-YF fp^ 

F!(r-fiu 


Pi y 

Y! k (F + !),_, 


(c-F+1) 


n-c 


C) 


n n + c 


(c-F+l)„ 
p‘(c- Y)! 1^1 / 


P 

c7 


A” 


n! 


Y! (F+lX 


_y w==0 


n! 


P^c-T)! 


(1)„ (c-F + l)„ 


P 

c; 


A" 


F! (F + 1) 


c-F 


w! 


£1 (^-DL Afl..-Y^l;s£ 

Y! (I' + lLr A 


Hence 




«=o n! 
pMc-Y)! 


7! 


Y! (7+1) 


F 

2^0 


c-F 


l,c-Y + l;-;^l (2.4.15) 


(2.4.14) 
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jjr (l,l;Y + l;p) + 


OA 




dk 


n" n'^ 

yp-+p- 
Un'- J'! 


£11^ .F,fl,c-Y + l;-; 

n (r+1). , I 


Y!(Y + 1L 

V— — 

SnA"-' 

^Po Z— r-T 

n =0 ^• 


C> 


y-1 n 

p Y— ^ 

0 ^ — 


n=o'(«-l)* 

Y-l 


p. z 


^(n-l)I 


um Po 


AP. 


^;i 


^ ,f,(i,1;Y + 1;p) 


,F,(U1;Y + 1;p)^ fr^ 


fI 

Y\ dX 


,F, (l,l;7 + l;p) 


(U)a (l,l;FH-l;p) lPa 


^X 


1 X 


[Y! 


Y! 


^ i p ^ 

\ ° 1^0 (P + l)„nl 
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,Fi (l,l;Y + l;p)^Po 


1 Y X 


Y-l 


Y! /i' J 


+ 


1 ^ tl -1 


, .( 1 ). ( 1 ). —nA 

P ; p V 

f! h (>' + !)» n'- 


P ^ ^ P "V ^ 

F, (l,l;Y + l;p)^^_— Po + 77 0 + (n.l)! 


2^1 


n'' _ 1-l P D "V (^)n-l (^)n-l P _ 

,F, (l,l;Y + l;p) Po + 777^ lr(r + 2 )„ , (n-l)! 


Y+l 


P 


{Y-\)l 


,F, (1,1;Y + 1;/))P„ + ^ + 1) § (P + 2), n! 


^( 2 ). ( 2 ). 


£- — - ,F, (Ul;Y + l;p)Po + -7777T ( 2 , 2 ;Y + 2 ;p) 

y _ IV (Y +U 


{Y-\y. 

Hence. 


IP. 


_£_ 


.e-,F,(l.l;Y + l;p) 

Y\ " * 


^ Y+l 


(P-l)! 


,F, (1,1;Y + 1;/>) + 7^ 2 F, (2.2;Y + 2;p) 


(///)AP„ 


2 ? 

( 2/1 


Y!(r+iL, ' I 






p!±iDL\p\X. 

Y\{Y + \l_, ; 


F„ l,c + l-Y;-;- 


(R) 
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= A 


dX 


— A" 


= A 


C -1 


cp 


A 


A 


A- 

dX 

d 


,Fo l,c+l-Y;-;- 


dX 


f. (l)n (C + l-Y)„ 
^ n! 


n=0 




(£ 

u. 

A"-' 

(1)„ (c + l-Y)„n-^ 


n==0 


nl 


. (l)„(c+l-Y)„ 

I— 

n-0 



(«-!)'• 


l.(c + l-Y) 2 

cj n=l 



(2)n., (c+2-F)„-i y— 

(n-1)! 


(c+i-n ^ S- 

\ C J n=0 


(2)„ (c+2-iO„ 



(n-1)! 


(c + l-y) 


oj 


R 2,c+2-Y;-;' 


2-* 0 
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Put in (R) we get. 




(c-Y)! (c + l-Y) 


Y'-iY^^)c-y 


,Fo 2,c + 2.Y;-;^ Po 


Hence. 


APn 




Y!(P+l)c-r ^ 


(c-Y)!c.p'=„ ^ 

° ^ ° 




l,c+l-Y;-;H + 


p--nc-Y)!(ct l:::^ ,fJ2,c + 2-7;-;^1 

c7l(7+l),_. ° ^ 


Hence. 


I=Pn 


y.^_p“+_^ ^F,(l,UY + l;-;p) + 


y+1 


(r+i)! 


,F, (2,2;7+2;-;p) + 


7!(F + 1),_, CA 


(c 7)l(c + l Y)P__ F f2,c + 2-Y;-;— 1 
c7!(7 + l),_y V 



(2.4.16) 
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Let a = 0, P= \ and Y = 0 then we get, 




'(n + l);i 
0 


and 


0<n<c 

n>0 


fn// 



0< n< c 
n>0 


Po'(t) = - 2Po (t) + //P, (t) n = 0 (2.4.17) 

Pn (t) = - [(n + 1) >^ + n //] P„ (t) + n A P^.j (t) + 

(n + l)//P„,i(t) l<n<c (2.4.18) 


^;(t) = - [(C + I) ^ +C//] Pe (t) + C2P,., (t) + 

c//P^i(t) n = c (2.4.19) 

P„'(t) = - [(n + l) A +c//]P„ (t) + nAP„.i (t) + 

c // P„+, (t) c + l<n<F + K (2.4.20) 

There fore the steady - state probaZ)ility - difference 
equations are: 

- IPo + //Pi =0 n = 0 (2.4.21) 

- [(n + 1) /I + n //] Pn + n 2, P^.! + (n + 1) // P^+i =0 1 n < c 

- [(c -M) 2 + c //] P^. + c>lPg.i +c//P^i=0 n = c 

- [(n + 1) /I + c //] Pn + nX Pn.i + c ju Pg^.! = 0 c + l<n<F+ K 


/ 

I 

(2.4.22) 

(2.4.23) 

(2.4.24) 


From equation (4.2.21) 

;U P, = A Po 

P, = pPo 

From Equation (4.2.22) 

-[(n + l);L +n/i]P„ +nlP„., + (n + 1) //Pn+i = ^ 
(n + l)//P„,, = [(n + 1) A +n/i]P„ -nAP„., 

(n + 1) P„.,i = [(n + l)p +n] ^ P ^n-i 

Put n = 1 

2?, = [2p+ 1]P, -pPo 
= (2 p + 1) p Po - P Po 
2 P 2 = pPo[2p+ 1-1] 

2P2 = 2 p'Po 
P2= P'Po 

Put D = 2 

3 P 3 = [ 3 p+ 2 ]P 2 - 2 pP, 

= (3p+2)p'Po-2p(pPo) 
3 P 3 = 3p^Po+2p^Po-2p^Po 
P3=P'Po 

Similarly, 

Pn ^ P " Po 0 < n < c - 1 (4.2.25) 

Where p = Xf p 


l<n<c 


From Equation (4.2.23) 

- [(c + 1) 2. + c a] P. + ^ ° 

c n P,^, = [(c + 1) 1 + C //] P, - C 2 P,., 

C P,,i=[(c + l)p+c]Pe-CpP<,, 

= [(c + l)p +c]p“ Po-^^P -P"'' Kl 
= p' Pq [(c + 1)p +c -c] 

= (c + l)p p' Pq 

f c + l^ 


n = c 


p = 


V ^ J 


PP^Po 





From Equation (4.2.24) 

-[(n + l)l +c;t]P„+n2P„.,+c/tP„,=0 c + l^n<Y + K 

c^P„., =[(n + l) /I +cp]P„ -nlP„., 

c P.«=[(n + ‘)P+c]P»-"^‘’-i 

Put n = c + 1 

C P„j=[(c + 2)/)+c]P„,-(c + l)pP, 

= + 2j p + c j f — ^ ^ pp Po “ Po 




= f^jp.p»P,[(c + 2)p + c 
p ^.p^(^ + 2)p 



Pc +2 =(c + l)(c + 2) 
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Put n = c + 2 


c P ,+3 =[(c + 3)p +c]P ,+2 -p(c + 2) Pc+I 
= [(c + 3)p +c]P ,^2 - (c + 2)p Pe+, 

= [(c + 3)p+c] (c + l)(c + 2)p'('-*^] Po- 

(c + 2)p. 

/ \ 2 

= (c + 1) (0 + 2)^*= — [(c + 3)/? +C -c]Po 

V c y 




VC 


Pc +3 = (c + 1) (c + 2) (c + 3) p ' 

Similarly, 

P. = P,H.-.) = (c + l)(c + 2)(c + 3)--(.Q-n)p' 


\ n-c 

£- 

K C > 



^„=(C+1)„-C P' 


The boundary condition y^.Pn = 1 


c<n<r + K (4.2.26) 


0=0 


V=Zp”+P'E(‘'+^)n.c 

n =»0 


CJ 


c -1 ” 


(c + l)„F (1)„ 


VC 


n =0 

c -1 


n =0 


n! 


p:'=Y,p'*p‘ jF. 

n»0 ^ 


L = ;iP„ p„-' 


n -» n + c 


(4.2.27) 


PA 


= AP„ 


= P 

* n 


^2, 


i;p’+A>' .F.fl.c + l;-: 


, n =0 


.P 


<=■1 M 

;i + .F. 

n =0 


cy j 


l,c + l;s— 1 a-^ c. A“'' + 

V cy p 

78 


+ ;ip‘ 


dX 


0 )n(C + l)„ 


-1 

'.cj 


n=0 


n\ 


= p« 


c-l 


/?” + 2^0 [ 


, 11=0 


c. p" + 


P‘ s- 


(I).(C + 1)..I1 


\n 

P 

u 


n=0 


n\ 


= Pn 


c-l / 

^n.p“+ 2^0 I5C + I;-; 


n=0 


C. p"" + 


p‘i;- 


(1). (c + l). 


n=0 


(n-l)! 


X= Pn 


j)n. p" + 2 F 0 fl;c + l;-;-l c. p“ + 

_n=0 V 




. (2),,, (c + 2),,, I £ 

C / 


, 71-1 


i= p» 


(«-l)! 

^n. p" +cp" 2^0 fl;c + l;-;— 1 c. p" + 


Ln=0 


,c+l / 




^2;c + 2;-;-l (4.2.28) 

c j 


The machine availability (rate of production per Machine) is. 
M.A. = 1- L/ii: 

ITje operative efficiency (utilization) is. 


c-l 


O. E. = 1 - X(l-nM P. 


n=*0 
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The Transient Behaviour of the Machine Interference 
M/M/C/K/N (with balking, reneging and spares) 


INTRODUCTION: 

In this ch^ter we consider a queueing model M/M/C/K/N with 
balking, reneging and spares and we assume a finite source (population) of 
N customers, c servers are available, customer’s arriving rate is X. We take 
Y spares so that when a machine fails, a spares is immediately substituted 
for it, if it happens that all spares are used and a breakdown occurs then the 
system becomes short. When a machine is repaired it then becomes a 
spares (unless the system is short in which case the repaired machine goes 
immediately into service). We analyze the two cases of spares, first of them 
is C<Y and the other C>Y. 

We assume that a queueing model M/M/C/K/N for transient 
behaviour with balking, reneging and spares. The arrival and service rates 
are and resp)ectively. The important works on machine repair problem 
have been considered by many researchers. Kleinrock [ 98 ] derived 
M/M/C/K/N for machine interference but without assumption of balking, 
reneging or spares. Kness have mentioned the transient behaviour of the 
repairman problem by using singular perturbation method to scale 
equations for the number of failed machine. Mokaddis et al discussed the 
cost analysis for two units warm stand by system with two types of repair 

fecilities. Jain and Singh studied finite queueing model with random feilure 
and delayed repairs. 
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In many realistic manufacturing / prodirction situation, due to long 
backlog of feiled units. Hie units may be discouraged i.e. either balked or 
reneged. Highighi et. al.[78] analyzed a multiserver Markovian queueing 
model with balking and reneging. The M/M/l/N queue v^th general balk 
function, reneging and an additional server for longer queue was studied by 
Abou-Ei-Ata and Kote [2]. Abou - El - Ata and Hariri [4] considered 
more general model M/M/C/N queue with balking and reneging. Abou - El 
- Ata and Showky [3] discussed the analytical solution of the single server 
Markovian over flow queue with baDcing, reneging and additiomi server 
for longer queues and A.I Showky [132] considered M/M/C/K/N machine 
interference model with balking, reneging and spares. 

We develop Uansient behaviour of the machine interference : 
M/M/G/K/N(with balking, reneging and spares) by using birth death 
process and Laplace transform technique and used boundary condition. 
We discuss two cases Y^and Y<C and also discuss the special case by 
using the hypergeometric distribution of first kind and second kind. 


3.1 CASE -I 


In this case. «e study the transient solution of machine interference system 
: mm/C/K/N with balking reneging and spares Y S C. then the set of birth death 

coefficients are as follows 


We consider : 



NX 

0<n<C 

^n=- 

, NfiX 

C<»<r 

+ 

1 

y<n<y+N 


1 0’ 

n>y + K 

and 





Q<n<C 

/tn = 

1cp + («-C)a, 

C^\<n<y + K 


In the usual arguments of the 8 - technique the probability differential difference 
equations are : 

/^(0=-mPo(t) + Ml*i(t), n = 0 (3.1.1) 

p„'(0= - m +n It) Pn(t) + NX Pn-i(t) + (n+1) MPn+i(t), 1< n < C (3.1.2) 
p;(/)=-(NpX + Cp)Pc(t) + mP<>i(t) + (Cp + a)Pc+i(t), n = C (3.1.3) 

_(NpX + C\i + (n-C) a) Pn (t) + NpX Pn-i (t) + 

[Cp+(n+l-C)a]Pa+i(t), C+l<n<Y (3.1.4) 

P;(0= -[(N+Y-n) px + cp + (n-c) a] Pn(t) + (N+ Y-n +1) pXPn-i(t) + 


[Cn+(n+l-C)a] Pn+i(t), Y + 1 <n < Y + K 


(3.1.5) 


-[Cn + (Y+K-C) a ] PY+K(t) + (N-K+1) Py+K-iCt), n = Y + K 

(3.1.6) 


Talcing the laplace transform of the differential - difference equations then we 
have. 

S Po(s)-P(o) =-m Po(s) + p Pi(s), n = 0 (3.1.7) 

S Pn (s) - Pn (o) = - (NX + np) Pn (s) + NX Pn-1 (s) + (n+1) \l Pn+1 (s) 

l<n<C (3.1.8) 


S Pc(s)-Pc(o) =-(N|iX + Cp) Pc(s) + NX Pc-i(s) + (Cp + a) Pc+i(s),n = C 

(3.1.9) 

S P„(s)-Pn(o) = -[NpX + Cp + (n-C)a] Pn(s) + N|3X Pn-i(s) + 

[Cp + (n+l-C)a] Pn+,(s), C+l<n<Y (3.1.10) 

S Pn(s)-Pn(o) = -[N + Y-n)pX + Cp + (n-C)a] Pn(s) + (N+Y-n+1) 

PX P„.i (s) + [Cp + (n+l-C) a] Pn+i (s),Y+l < n<Y+K (3.1.1 1) 
S Py+k(s)-Py+k(o)= -[Cp + (Y+K-C)a3 Py+k(s) + (N-K+1) pX Py+k-i(s) 

N=Y+K (3.1.12) 

Define boundary condition 

P„(o) = S P„., (s) (3.1.13) 

From equation (3,1.7) 

S Po(s)-P(o) =-NX Po(s) + |i P,(s) 
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^ Pi(s) =(NX + S) Po(s) 

P,(s)= (m/n+ S/n) Po(s) 

P, (s) = (Np + (j)) Po(s) p = X/fi, <t> = S/p 

From equation (3.1 . 8 ) 

S Pn (s) - Pn (o) = -(NA-+ np) Pn(s) + NX, Pn-l(s) + (n+1) p PiH-l(s) 

(n+l)p Pn+l(s) = (NX+np + S) Pn(s)-NX Pn-l(s) -Pn(o) 

(n+l)p Ptrfl(s) = (NX+np+S) Pn(s)-NX Pa-l(s) - S Pn-l(s) 

(n+l)p Pn+1 (s) = (NX+ np + S) P„(s) - (NX + S) Pn-i(s) 

(n+1) Pn+i(s) = (Np+n+<|>) Pn(s) - (Np +(^) P„.i(s) 

Put, n = 1 

2 P 2 (s) = (Np +1+ Pi (s) - (Np +<|)) Po(s) 

2 P 2 (s) = (Np +<t) + 1) (Np +^) Po(s) — (Np+^) Po(s) 

= (Np+(|))(Np+<t)+l-l) Po(s) 

p,(s) = 

Put, n = 2 

3 P 3 (s) = (Np+2+(j>) P 2 (s)-(Np+(i)) P,(s) 

3 P3 (s) = (Np +,(1 + 2) - (Np -H|,) . (Np -HI,) ?,(i) 

= (Np+(p+2-2) 
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m 


3 P3 (S) = 


P3 (S) = 


Similarly 


(s) = i^P + <l>) ..p^(^s), 0<n<c 
n\ 


(3.1.14) 


Where, p = X/p, ^ = S/p 
From equation (3.1.9) 

S Pc(s)-Pc(o) =-(Np?. + Cp) Pc(s) + NA. Po,(s) + (Cp + a) Pc+i(s) 
(Cp + a) Pc^i (s) = mx + Cp + S) Pc (s) - NX, Pc-i (s) - Pc (o) 

= (NpX, + Cp + S) Pc(s)-(NX, + S) f^iis) 

(S+ 1) Pc+i (s) = (NPy + 6 + V)/) Pc(s) - (Ny + ip) Pc-i (s) 

Where, X/a = y, Cp/a = 8, s/a = vj/. S/p = <j) 

0+1 0+1 

P„,(s)= im+S+v) (Np+^y ^ (yy+y) ^f/p+0-' ^ 

S + \ Cl S + l (C-l)l 


{Nfir+S+yr) iNp+0 - _ (NA/a+S/a) Q[pjjr±^ . 

S+1 a ^ (C-1)!C ^ 


(Nj0y+S+v') {Np + <j>r . (iV;H-5') iNp+0-^ _ . 
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N^Y+S+y/"\ {Np+^l>y s, , (J^P±fL<MP±ftLSEp(^\ 

“Tm J C! ^ + l C! a “ 


NpY+S^w'\{Np+ff . iNp^(l>y S j . 


P(s) [NfiY + ¥ + ^-‘^ 

a S+\ 




NPy+W ^ 
s+\ } 


From equation (3.1.10) 

S Pn (s) - Pn (o) = -[Np?i + C|i + {n-C)a] P„(s) + NpX, P„-i(s) + 

[Cn + (n+l-G) a] Pn+i(s) [C^x + (n+l-C) a] Pn+i(s) = [Np>. + C^i + (n-C)a+S] 
P„(s)-Np?l Pn-l(s)-Pn(o) 

[C^i + (n+l-C)a] Pn+i(s) = [NP^ + C^i + (n-C)a+S] Pn(s) - (NpX+S) Pn-i(s) 

[5 + n+1 - C] P„+,(s) = [Np Y + 5 +n - C + 'P] Pn(s) - [Np y +W] Pn.i(s) 
Putn=c+1 


[(? 4- c + 1 + 1 - = [N^y + <5 + c + l- c + - (NfiY + W)Pci^) 

[5 + 2] Pc^2{s) = [N P+5 +4^ + 1] Pc.i(s) - [NP Y +'P] Pc(s) 


[Npy + S+'P+l] P 


V Po(s)-W;^ + V^]^i^^^Po(^) 


iNp+<pr 

a. 


■-{NPy + ¥\ 


NPy + <5^ + V' + 1 
J+l 


iNp^0 

Cl 


P,is) [NfiY + ¥] 


Ny9/+<J+^+l— <5— 1 

_ ^ _ 
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_ {Np + 0 {NPy+ip'f 
C\ ^+1 “ 


(,, ( 5 +ix<?+ 2 ) 


Put n = C + 2 

[5 + C+2+1-C] Pcf3(s) = (Npy +5+C+2- C+'F) Pe^2(s) - (Npy +'F) Pch-i(s) 

Y+l<n<Y+K 

[6 + 3] Po. 3 (s) = (Npy+6+2+'F) Po. 2 (s)-(NPy+T) Pc>^i(s) 

= (NPy +6+2+VP) Po(^)- 


(NPy+T) 


= Wr+y^L p (Npy + 6 + 2 +^- 6 - 2 ) 

C! (^+1X^+2) ^ 

(iVp+^X {NPr + y/f 

C! (<5 + IX^ + 2) 


Pc.3(^) 


iNp + <l>Y {NfiY + wf 

a (<j+ix^+2x<y+3) 


Po(^) 


Similarly 


P«(S)= Pc+(tM:)(s) 


^^[E±f}l.WZly^P^(s), C +1 < n<Y ( 3 . 1 . 15 ) 


From equation (3.1.1 1) 


S P„ (s) - Pn (o) = - [(N + Y-n) PX. + Cn + (n-C)a] P„(s) + (N+Y - n + 1) p^i Pn.i(s) 


+ [Cn + (n+l-C)a] Pn+i(s) 


[C[i + (n+l-C)a] Pn+, (s) = [(N + Y-n) + C^i + (n-C)a + S] P„(s)- 

(N + Y-n+l)pX, Pn-i(s)-Pn(o) 

(6+n+l-C) Pn+i (s) = [(N+Y-n) Py +5 +n-C+'P] Pn (s) - [(N+Y-n +1) py +^] Pn -i(s) 


Putn = Y +1 

(5+Y+l+l-C) Py+hi(s) =[(N+Y-Y-l)py+8+Y+l-C+'P] Py+i(s)- 

[(N+Y-Y-1 +l)py+'P] Py+i-i(s) 


(5+Y+ 2-C) Py+2(s) = [(N-1) Py +6 +Y+1-C+'P3 Py+i (s) - [(NPy +'F] Py (s) 


= [(N-1) py +5 +Y+1-C+'P] Po(s) 

Cl (c>+l)r+i-c: 


C! {d+\)Y^ 


iNp+<pf 


C! 


(^+1) 


r-c 


\N-\)fiY+S + Y + \-C + Y 
S + Y + \-C 


jNp+^f iNfiY+¥y^''^ 


Po(s) 


C! {S + l)y^ 

XN-l)/3Y+S+Y+l-C-i-i/^-S-Y-l + C ' 

<^+r+i-c 


.iNp+4f iNJ3Y+w) 


T+C-i 


a 


(S + l) 


Po(s> 


r~c 


(N -\)Py'^¥ 


J + F + l-C 


^+2 (“S) 


_ {Np-\-ff jNfiY+yy'^ QfPY+V){iN-\)fiY-¥y/\ 
C\ (^+1),.^ (<5+T + l-CX^+r+2-0 


Po(s) 
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Putn = Y+2 


(6+Y+2+1-C) Py.3(s) = [ (N + Y- Y-2) Py + 5 + Y+ 2 - C + \i/] Py.2(s) - 
[(N+Y-Y-2+1 )Py + \i/] Pyfi(s) 

(6+Y+3-C) Pyf3(s) = [ (N-2) Py + 5 + Y+ 2 - C + \i/] Pjh-2(s) - [(N - 1) Py + ¥] 


Py^l(s) 

= [(N-2)Py + 8 + Y+2-C + \i/] 

{Np + (i>f {.NPr^yr)[{N-\)fiY^¥] Po(s) - [(N-D Py + vl 

C! (<y + l)y-c is ^Y + \-CXS + Y + 2-0 


r+-i-c 


{Np + ff {Nfiy + xf/) 


Po(s) 


Po(s) 

C! (<> + l))'+i-c 


— l^fiy + yz + S + y + 2 — c 
6 +y+2—c 


= [(n-i)Py+h/] 


iNp + <l>f {NPy + O 
C\ iS + \)y^: 


r+i-c: 


[iN-2)fir + ¥] 
S+Y+2-C 


Po(s) 


P y+3(s) 


(Np+f) (Nfiy + O jNfiy + QiiN - \)Py + y^][(Y - 2)fiy + y/] 


a 


iS + l)r 


iS + Y + l-CXS + Y + 2-CXS + Y + 3~0 


Po(s) 


Py+3 


_ (Np+^fiNfiy+O 




0y 


fiy 


C!(«^+l) 




iS + Y+\-CXS+Y + 2-CXS+Y+3-0 


Po(s) 


Similarly 


Pn(s)= PY.(n-Y)(s)= 

W + l),_e iS^Y+\-0„.y 


Po(s) 
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( 3 . 1 . 16 ) 




ais+DrA^+r+i-c),^ 


Pr — Po(s) 


Where, IN+ 





fif" ' ' Pr 

Thus ^(5) can be written as follows : 


p^(s), 0<n<c-l 


n\ 


Pn(.S)= { 


(Np *0 {Nfir + V'rZ .-^j.y .<„<Y 
C! (<; + i)„ 




(Np+<f>nNj 3 r+fi^y~‘iPr) 


N + 


w 


(3.1.17) 


C!(^ + l)r-c(^ + l" + l-c)-r 


P^is), Y+l <n<Y+K 


y+k 


To find Po, the boundary condition; '^P„is)=l 


nssO 


C-1 


\n— C 


S Po (s) + 2 

«=0 w*c 


{Np + <l>) (Y/?y + r)" p.(s) + 

C! (<5 + l)_ 


iNp + <l>f{NPr + y/) 


Y-C Y+K 


iPr) 


n-Y 


N- 


YL 

fir\. 


C!(^ + Dr^ S‘■.^ iS*Y^\-c)^j 


-nw=i 


P„-1 (s) = V CA^±fil+ (^P+# y mr+wT^ + 

h= 4) w! C! *=c (<^'*"l)n-C 


CI(^ + 1),^ (5 + F + l-c)_, 


iPrT^ 

iV + ^ 



L /^J 

:W-r 


(3.1.18) 
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n wi+C 




(S + l) 


n-C 


l^ iNpy+y/)" n\ _ y (A/'^r-t-V^)"0)» = jFi (1; 5+1; Npr*- 'F) 
^ iS + l)„nl S (^ + 1)„«5 


(Np + (NjSr+rf'^ (firr' 

pr 

Cl(S + 1 )}^^: + 7 + 1 — C)„_Y 


Y+K 


(n) I 


n=:7+l 


(Np+^fiNfir+wY^ 


Pr 


aid + 1 ),^ (J + F + 1 - C)„_r 


jNp 4- (ftY iNfiy + y/) 
C\iS^\)r_c 


(firr'iN+^U 


ti iS + Y+\-C) 


lH-1 




ll+l 


C!(<? + 1) 


Y-C 


^ (d+r+i-c)„,, 


/.r ,y-C [^ + ^ 10 »F) [(^- 1 ) + -^]„ 0 »F)”( 1 )„ 

Qip+^YiNPy-^y^Y^ Py y Py_ 


Cl(^ + 1 ) 


r-c 


iS + Y + \-OU iS + Y + 2-C)„t^. 


(N.i)+i; 8 +Y+2-C; ^y) 

C!(d + l)y_c+, PiK 


Po'* (s) - ^ (iVp+^) " + 


C-1 


(ivp+#^ (jyygr+r)'^ 


«=o n! 


C! 


n=s=C 


(^+ 1 ) 


If-C 


iNp+ijifiNpr+y/) 




0 (^ + 1 ) 


T-C 


iS + Y + \-0^y 


Hence 
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y We±C+<^^i^iFi(l;S+l;NPY + ^ 

iNp+<l>f {Npy^y/^'^^\ . ^.j^+JL; 5 +Y+2-C; ^y) (3.1.19) 
C!(<J + l)r_c+i fiy 

Where {l)n = 0 when n>Y"C in the first hypergeometric function and either 

( 1 )^ = 0 or (1-N)n=0 when n>K-l in the second hypergeometric function. 

To calculated the expected number of units in the system, a result due to 
Abou-El-Ata is used. Which states that for a simple birth death process. 

T = 

~ dX ^ dp 


Thus L = XPo(s) 


dPo(s) 

dX 


=^Po(s)^ 

OA 


^tid d (d+ly^, 


(I) p„(s)f 

dJL n! dA 


ni 


ZS ra 


Po(s) 


dX 




=^Po(s) 


c— 1 


5]«(A2 + s)”"' 


B=0 


IJV 

n!p”_ 


(3.1.20) 
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(n-l)! ^ 




a=A./m ri = S/n 


X n(.s)-^\^^^^^,F,(l-,S + i;Nfir+<l') 


=X P.M [iFi (t; 8+1 ;NPy+v)] 


d (Np+(^y 
dX c\ 


Xp^(j)WF ±£1 ^[,Ffi-,S + V,NPr+¥)\ 


(II) a 


?.?.(.) Po(s) 

oX Cl 


(iv'~+-r 

a // /i 


ai c! 


a Po(^) — 


JL ^ i^x-^sy 

n‘ dx cl 


aPo(.) — — c 


// (c-1)! 


=Na Po(s) 


(Na + rjf 
(c-l)! 


?„ (5) No 


{Na+rff 

(c-l)! 
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(11) b XPo(s) [|F|(l;^+l;AW'+y)] 


oA 


cl 


f jium+y'y 

S iS+\)„nl 




dX 


„ (i)„w-+-r 

Y g a 

”=» (c^ + l)„/i! 

a 


iNp+<i>r 


'kPois) 


^ i\\n{NpX + sr-- Np 


c! or (<5 + !)„(«— 1)! 


NpTii,(s) Z 

C* n=l 


(i)„(iV;gr+tyr' 
tf(^ + lX<^ + 2)„_,(«-l)! 


(ivp+<^r 


cl 


nPyPo(5) S 


(i)„AN/3r+¥y 

^(S + lX^+2)„n! 


(Np+^y Nfir 


n(s) Z 


(2)„(Nfir+V^y 


cl s+i ^ i’S+iyd 


{Np + ((>y Nfiy 
cl S + l 


P,is) ,FX2;S+2;Nfir+W) 


Hence 


Po(s) 


dX 


Q^p+^y 

cl 


iFi(l;S + l;Nfir+¥) 
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Na , F,(l;S + \;Nfir + ¥) Po(s) + 

^ (c-1)! 


... (Np + iY ,Fi(2;^ + 2;iV^7 + V^) ^o(s) 

{5 + l)c! 


(Ill) 


dX 


iNp+(i>r p(v(N-x\+(^y,s+y^2-c-,m 

c\ (^ + 1),^,. M 


iNp +0 ^x. — XPois) ■—T]^.Nfir+lpy~‘^] 

oA 


4 - 


cl 


iS + \) 


y-c+l 


{Np + ^y i^fir±¥yiyLA. F(V\ (N-l) + — \S + y + 2-c-,yr)^o(s) 
c! (J + IW dx\' ' 


(III)a 


yp,{s) 


A 

dX 


^ ^ vy-^+l 


a a 


X?.w + 


=Np- P.(s) (y-c+l) (wi+^)’~' 

a a a 


NPy Po(s) (y-c+l) (NpY+\i/)’ 


,3^ 


(Ill)b 


n=0 


(l)„[(A^-l) + -f^ln(^r)" 

. PY 

{S + y+2-c)„n\ 


Xpa) yO).K"-’)+r^<»l'kA(^A/a)- + 

(^+y+2-c).rf aT'^ 


X Po(s) .-i-[(Ar-l)+.^l. 


115=0 


(^ -f j + 2 — c) „ n! 5./^. 
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0)„(>9r) 


" ( 


i:^(S+y + 2-c)„rit[^ ply 


¥ 


h (<5+j^ + 2-c)„(«-l)! hiS^y^2-c\r& 


'-urn 


n~\ 


{5+y+2-c\(n-\)\ 


n =0 


(P+y+2-c)„nl 


S (^ + >' + 2-c)„^,n! S(^+3' + 2 -c)„w! 


(2)„i( - 1) + / >gy] [(i\^ - 2) + y. / >gr]„ (Py)" 
(^ + ) 

dU-iy/Jir 


= S rx^ \ I 

S (d + >' + 2-cXd + >' + 3-'c)„«! 


+^io(«)X 


tS(^ + r + 2-c)„/i! 


= P(s) M^U}tr/M ^p^(2,(JV-2) + ^fr/py^^ 

(o+y+2-c) 


X Po(s) iFi(i;<^+y+2—c;-^/A) 


Hence, 


L=/>o(j) 


\c-l 




(«-!)! 


(c-1)! 


(d+l)C! (d+l)r-c+i 


2F,(1 ,(N-l>H-ypy);5+Y+2-C;py> 


(iSTp+^X M(-^^-i)+<y/>gyl 


C! 


(<y+i) 


r-c+i 


(j+r+2-c) 


2F,(2,(N-2)+x|//|Jy);5+Y4-3-C;pY) ^iF^l ; 5+Y+2-C;- m/ZX,)! (3.1 J21> 


a 


(,S+\\ 


r-c+i 
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3.2 SPECIAL CASE 


p;(0=-mPo(t) + |.iPi(t), n = 0 

P’Qy= . (NX- +n ^i) Pn(t) + NX Pn.i(t) + (n+1) pPn+i(t), H < C (3.2.2) 

P:it)= -(NpX + Cn)Pc(t) + NX Po-i (t) + (C^l + a ) Pc+i(t), n = C (32.3) 

+ C[x + (n-C) a) Pn (t) + NpX Pn-i (t) + 

[C^i +(n+l-C) a ] Pn+i(t), C+1 < n< Y (3.2.4) 

-[(N+Y-n) pX + c^i + (n-c) a] Pn(t) + (N+ Y-n +1) pXPn-i(t) + 

[Cji,+(n+l-C)a]Pn+i(t), Y+ 1 <n<Y + K (32.5) 

(^)= .[CjA + (Y+K-C) a ] Py+k ( t) + (N-K+1) pX PY+K-i(t), n = Y + K 

(3.2.6) 

Let a = 0, p = 1 and n = k then the above equations are 
/?(0 = -KXPo(t)+nP,(t) n=0 (3.2.7) 

p;(/) =-(KX+nn)Pn(t) + KXPn-i(t) + (n+l)nP„+i(t), 1^<C (32.8) 

p;(0 = -(KX+Cn) Pc (t) + KX Po-i (t) + Cn PcH-i (t) , n = c (3 .2.9) 

p;(0=<KX+c^l)P„(t)4-KXP^-l(t) + c^lP^l(^^ c+l<n<Y (32.10) 

P:( 0 =A(K+Y-n)A+cp]PSt)HK+Y-n+\);LP„^(t)+c^^^ P»+i(0 Y+l<n<Y+K (3.2.11) 
PLk (.0 = " [cm- PY+K(t) +X Py+k -1 (t) n=Y+K 


(32.12) 


From equation (3.2.7) 


p;Q)=-KXPoit)+\im^ n = 0 
Taking laplace transform 
S Po(s)-P(0) = -KX Po(s) + P his) 
^iP,(s)= Po(s)[KX+s] 


Pi(s) 


KX + s 

L M J 


Po(s) 


Pi(s>= (Kp+<1>) Po(s) Where x/p-p. S/p ^ 


From equation (3.2.8) 

K(0 =: " + »M)Pn (^) + (^) + (" + OjUP^+l (t) 

taking laplace transform 

S Pn (s)-Pn (0) = - ( KX.+np) Pn (s) + KX, Pn-l (s)+ (n+1) p PiH^l (s) 
(n+l)p Pn+l(s) = (KX + np+s) Pn(s)KX, Pn-l (s) -Pn(0) 
(n+l)p P^i(s) = (K7. + np+s) P„(s)- KX. Pn-i (s) - S P„.i (s) 
(n+l)p P^i(s) = (KX, + np+s) P„(s) -(KX,+s) Pn-i(s) 

(n+1) Pn+l (s) = (Kp +^) Pn(s)- (K-P +^) Pn-l(s) 

Put n=l 

2 Pi ( s)=(Kp + 1 + <t>) Pi (s)- (Kp +<|>) Po (s) 
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> 2 ?2 (s)=(Kp +!+<!>) (Kp +4>) Po (s) - (Kp +<!>) Po (s) 

= (Kp+(l)) Po(s)[Kp + (j)] 

= (Kp+(|))(Kp+^) Po(s) 

P2(S)=^^^ 

Put n=2 

3 pj|'s) = (Kp + 2 + <i)) P 2 (s)-(Kp+<|)) Pi(s) 

= (Kp + 2 + <})) Po (s) - (Kp +<t>)(Kp +4») Po (s) 

3 P 3 (s) -= [ Kp + <|>] Po (s) 

= p«(*) 

Similarly 

p (s) = Po(s) , 0 < n < c- 1, Where p= X/ji , <1> = S/p 

^ ^ n! 

From equation (3.2.9) 

(t)=-iKA + cf2)PAt) + K^c-x(0-^cMPcAt) n = c 

Taking laplace transform 

S Pc(s)-Pc(0) = -(KX + cp) Pc(s>+Ba PoKs) +cp Pcn(s) 

cp Pc+i (s) = ( KX. + cp +s) Pc(s)-KX- Pol (s) -Pc(0) 

cp Po-i(s) = (KX + cp+s) Pc(s)-KX Poi(s)-S Poi(s) 


(3.2.13) 
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CH Pcti (s) = ( KX + on +s) Pc(s)-(KX+s) Pm(s) 
c Pert (s) = ( Kp + c +^) Pc(s)- (Kp +*) P<rt (s) 

c Pert (s) = ( Kp + c +*) p„ (s). (Kp +(,) (^+^r‘ p„ (s) 

cl (c-i)! 

C PePl (s) ( Kp + C +(j)-c) 

d 

p^,(,)=(&t^i 

cl c 

From equation (3.2.10) 

p'.(t) = - [ ia+ cp ) P„ (t) + KX P„., (t) + cp P,H-I (t) 

Taking laplace transform 

S P„(s)- P„(0) = -(!a + cp) P„(s)+KX P„.i(s)+cp P,Hi(s) 
cp P,+,(s) = (KX + cp+s) P„(s)-K>. P„.,(s)- P„(0) 

Cp Pn+1 (s) = ( Ba + cp +s) P„ (s)- KX Pn -1 (s) - S Pn -1 (s) 
cp Pn+i (S) = ( la + cp +S) Pn (s)- (KX+s) ?„-! (s) 

Putn = c+1 

cp Pc +2 (s) = ( KA- + cp +s) Pc +1 (s)- (KA+s) Pc(s) 

cp Pc^ 2 (s) = (KA + cp+s) Po(s)-(KA+s) Po(s) 

clc cl 

c Pe. 2 (s) = (Kp +C + 4) P„(5)-(Kp + 4,) Pj(s) 

dc cl 
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=( Kp + c + 




Po(s)- 


Po(s) 


ClC 


Po (s) [ Kp+(i)+c-c] 


dc 


C P Cri-2 (s) 


c\c 


c+2 


Po(s) 


Pc+2 (S) 


_ iKp + </>) 


c+2 


Po(s) 


ClC 


P c+2 (s) 


(Kp^^r^ i 


cl 


Po(s) 


Put n = c+2 

cp Pc+3 (s) = ( KA, + ep +s) Pc+2 (s) - ( KX- + s) Pc+i(s) 




(KX + on+s) Po(s)-(KX + s) 


cl 


\c+l ( I 


- Po(s) 


C P<^ 3 (s) = (Kp +C + 




+ Po(s) [iC/7+<i + c-c] 

cl Vc) 


= (gp + C lfll p„(s) (,Kp+^) 


Cl 


P.*3(s) = ^^2:^(^iJ P„(s) 


Similarly 


Po(s) 
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P.= ?c+(n-c)= P„(s), C+l^n<Y 






From equation ( 3 . 2 . 1 1 ) 

p'^{i)= -[{K + Y -n)X + c^]F„ (t) + {K + Y-n + l)2P„_i (0 + c/iP„+, (0 

Taking laplace transform 

S Pn (s)- Pn (0) = - [(JS: + F - «)>^' + C//]P„ (s) + iK + Y-n + l)AP„_j C^) + CfiP^i (s) 

Cp Pn+i(s) = [(K+Y-n)X. + c}i+s] Pn(s)-( K+ Y-n+1) X- Pn-i(s)-Pn(0) 

C\l Pn+i(s) = [ (K+Y-n)?l + Ctl+s] Pn(s)-(K + Y-n+l)X, Pn-l(s)-S Pn-l(s) 
C\i P„+i(s) = [(K+Y-n)A. + cn+s] P„(s)-(K+ Y>n+1) X- Pn-i(s)-S Pn.i(s) 
Q (s) = [ (K+Y-n) p+ c + <!)] Pn (s) - ( K + Y-irFl ) p P^-i (s) - <1> Pn-i (s) 

Put n = Y +1 

C Py-p2(s) = [(K+Y-Y-1)p + c +(!)] Py+i<s)-(K)p Py-.i-i(s)-<|» PY+i.r(s) 
C Py+2(s) = [(K-1)p + c +<|>] Py+i(s)-(K)p Py(s)-<|> Py(s) 

C Py+2(s) = [(K-1)p + c +<|)] PY+i(s)-(KpH-<|>) Py(s) 


[(K-l)p + c +<!)] 


^ ^0" P„(s) 


= L Kp->-^ !Z1 ( iY*" P„(s) [^-1) 


C Py+2(s) 


d 




Po(s)[(K-i)p+iH 


Py+2 (s)- 


Py+2 (s)= 


d 


y-c 




(Kp+(t»)[(K-l)p+(l)] Po(s) 


fnVr''"' 


c\ 


c) \c 


(Kp+(|))[(K-l)p+<|>] Po(s) 


Put n = y + 2 

C Py +3 (s) = [ (K-2) p + c + <|)] Py +2 (s) - ( K4-Y-Y-2+1 ) p Py+i (s) - <j> Py+i (s) 
C Py +3 (s) = [ (K-2) p + c + (|>] Py +2 (s) — ( K-1 ) p Py+i (s) — <|> Py+i (s) 

= [(K-2)p + c +(|)] Py+2(s)-[(K-1)p+<|>] Py+i(s) 

(K/>+, 


[(K-2)p + c +(|)]- 


.K-c ^lYrn 


VC 


'"iY''(Kp1-it.)[(K-l)p+« Po(s) - 

Kc) 


[(K-l)p+(H 




c! 


Kcj 


Po(s) 




2 

rn 

cl 

Uj 


Uy 


c Py+3(s) 


( K p+ <1) ) [(K-1) p+<l)] Po(s) [(K-2) p+(^+c-c] 




c\ \c) \c 


( K p+ <1> ) [(K-1) pf<|)][(K-2) p+(l)] Po (s) 


Py+3 (s) 




d \c j Vc 

Hence 

Pn — P Y+(n-y) ~ 


( Kp+ ♦ ) [(K-1) p+«[(K-2) Po (s) 


cl I c 


y -c / - \ n—Y 


XKp+^y- fiY^ro 


[Kp+«..Y Po(s) 


dc^-^ U 


[Kp+(|)ln-Y Po(s) 


Pn(s) \'’-^[K^^/p]„.y Po(s), Y+l^n<Y + k 

dd \cj 


Thus Pn can be written as 


(Kp+^r 


Po(s) 


c! U 


(Kp+^f-^ fir" 


0^<c 


Po(s) C+l<n<Y 


p^-^lK + <f>fp}„.r Po (s), Y+1 < n< K 


Where p = Xyp, (})- S/p 


p„-. (3) = t . I, ^ 


_a(A^+«()' If [(*.,£] (p 


«=0 ^ ^ ^ 


CIX? n^-k-l P 


t{^ 


n =:>'n+c 


y f S±iY gf:^±iTfliL= ,F„ (1; -; (Kp + <t.)/c) 

w*0\ J ^ «=sOV ^ 


U.11 f^rn 


-►n+Y+l 



Kp + ^ 

C J w=0 


K^\ 


P yc) 


n\ 


V 


c ) P ^ 


po-' (s) = 'F" ( i--;(Kp+*) •" [^ 


2Fo(i,(K.-i)+<l>/p;-;-f) 


Hence 


p -» (s) = Y^M±^L+L jFo ( l,-;(Kp+(^) /c) — 2F0 ( ^»(K-l)+ Vprif ) 

S n\ c\ 

( 3 ^. 18 )^ 


L = 9 i Po(s)^n''(^) 

ox 


L = ?i Po(s) 


dX 




n=0 


(I) % Po(s) 


dX 


ri d 


nssO 


d 






P ^ } 



X Po(s) ~ 

nl oA 


n=0 


X Po(s)- 

n\ 


^ n{KX + sy~^^^ 
Za 

rt^O M 


£(1^^^)"-’ ^ —1^ Po(s) 


?J=0 


// («-!)! 


= |(^p+^^r'(Kp)— Po(s) 




(II) X Po(s) 


dX 


\,F,{\---,{Kp + <j>)lc 
c\ 


= X Po(s) 


dX 




► nmO 


X Po(s) 


dx 


1 ^ 




KX + s 




rfSV J 


= X Po(s) 


1 ^ (iQ+s) 


ii-i 


X Po(s) 


1 Y \ F- ) 

cits 


X Po(s) 


dh yc” 


Po(s) 


c’s” M C" 


= Po(s) 


„ {Kp+(i>y-^ 




cl i»=0 


Po(s) 


I 

zX ”'^ 


C\C^^ 




- Ko « <'H^ ' " 
PoW^S— ^ 

c!c n\ 


«-i 




0)« 


V C . 


,/I-l 


c!c tf («-l)l 


p«(s)^S- ^ 

c!c 


ni 


(2).f^ 


cic ;:5 «! 


= P<,(s)S,F(,(2;-;((Kp+4i)/c)) 

CIC 

(ni) 2 ?.wA i,(k-i)+Vp;-;^) 

5a dc c 

PtiPi^ L ± 

d/?^ InX 


ex 


[2Fo(1,(K-1)+Vp;-;^)1 
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(I„)a = 


K~+(f) 

. . 


r-c+i 


= (F-c + 1) 


zr^ . 

K — Vip 


-\Y-c 


K 


= (F-c + l) iKp+<l>y-^ Kp 


(IIDb ^^bFo ( 1,(K4)+ <|)/p;-;^)] 

OA C 




{\)MK-\)+<l>lp]„ 


n\ 


dX. 


«=s0 


{K-\)+- 


^1 f aT 

^ \Xf^j 


n=sO 


{K-D^- 


4>y- 


d ( A 


A j„ dA\pc^ 


f-T 


A.y 

dA^ 


iK-l) + - 


n=0 


iK-\) + 


<t>p 


n(Ay-'—^ + 


{per 




Kfic) S 


,v 


0)„ 


n\ 


n=0L 




f 1 A" 


rti 


^ 1 v/^y 


r 1 






-^p 


0)« 


w=0 


[(JC-D+^/pU^ (1). 


.V 


«S =0 


WI 


— A-\ j. ,F.(i;-;-fi4) 


. [(i:-i)+^//jl.^£j 0). 


■nsl 


(n-1)! 


+1^] ;i,Fo(l;-;-#^) 
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tA 






P 

cj 


,Fo(l;-;-^^) 


-[{K-\)+(j>ip]{^ MUK-2)+(i>ip--\pic)+x .FoO;-;-^) 


Hence 
L=?o(^: 


£r-l 




{Y-c+\)KpiKp+<l>y-\F,iUK-\)+^l--A) 

p c 


+ 


C!C 


Y-C+l 


, (px 

UJ 


3Fo(l,(is:-l)+-^;-;^)+ 
p c 


\Y-c+\ 


(Kp+0 - ,^{p 


c\c 


CJ 


P 


clc 


(^_i)+i 


FJ 

^Cj 


,Fo(2,(F-2)+-^;-;^) 

F c 


( 32 . 19 ) 


33 Case - II 


This case treats the system of machine interference ; 

M/M/C/K/N with balking, reneging and spares Y<C and hence the birth-death coefficients are . 




and 


NX 

N + Y- n^X 
{N + Y-n)pX 
0 


0^<Y 

Y<n<C 

C<n<Y+K 

n>Y+K 




\W 


0^<C 


C;£/ + (n-c)a, C+1^<Y+K 

Then as before the probability differential difference equation are 
-NXPo(t) + p.Pi(t) n = 0 

p;(/)= . (NX +n n) P„(t) + NX P^,(t) + (n+1) 1^<Y 

P„'CO= -[(N+Y-n) X+np.]Pn(t) + (N+Y-n+1) XP^i(t) + (n+1 ) pP„n(t), Y+l^n<C 
= -[(N+Y-c) pX+Cp] Pc(t) + (N+Y-C+1 ) X Pc-i(t) + (Cp+a) PcH-i(t) n=c 

?;(/)= -{(N+Y-n) pX+Cp+(n-C)a] Pn(t) + (N+Y-n+1) pXP,>i(t) + 

[Cp+<n+l-C)a]Pn4-i(t), C+l^Y+K 

p ;^^(0 = -[cp + (Y+K-c)a]PY+K(t) + [N-K+l]pXPY+K-i. n = Y+K 

Taking laplace transformation then the differential equations are : 

S Po(s)-P(0) =-NX Po(s)+li Pi(s), n = 0 

S Pn(s)-Pn(0) =-(NX + np) Pn(s) + NX Pn-i(s) + (n+l)p Pirt^i(s), 1 <n< Y 
S P„(s)-Pri(0) = - [(N + Y-n) X + np] Pa(s) + (N+Y-n+l)X Pn.i(s) + (n+l)p P„^i(s) 


(3.3.1) 

(3.3.2) 

(3.3.3) 

(3.3.4) 


(3.3.5) 

(3.3.6) 

(3.3.7) 
(3.3-8) 


UO 


Y+1 < n< C 


(3.3.9) 


S p^(s)_p,(0) = -[(N + Y-Qp^ + Cp] Pc(s) + (N+Y-C+l)X Poi(s) + 

[Cp + a] Pcrt(s), n = C (3.3.10) 

S (s) - p„ (0) = -[ (N+Y-n)p?i + Cp+(n-C) a] Pn (s) + (N+Y-n+ 1 ) PX P„.i (s) 

+[Cp+(n+l-C)a] P.^i(s), C+1^<Y+K 

(3.3.11) 

S Py+k(s)-Py+k( 0)= -[Cp + (Y+K-C)a] Py+k(s) + (N-K+l) pX Py+k-i(s), N=Y+K 

(3.3.12) 


From equation (3.3.7) 

S 'Po(^)-P(^) Po(s) + P Pi(s), 


(3.3.7) 


p P,(s)= (NX+s) PoXs") 

^(s)= (Np+<^)Po(s) 

From equation (3.3.8) 

S Pn(s)-Pn(0) =-(N>. + np) Pn(s) + N>, Pn -1 (s) + (n+1) n Pttfi(s), 1 <n<Y (3.3.8) 

S P„(s) - SP„.,(s) = -(NX + nn) P„(5)+NXP^,(5) + (n+1) n P^,(j) 

(n+1 ) p P,^, (s) = (NX + np+S) P„ (s) - ( NX+S) P^i (s) 

(n+1 ) P„^, ( 5 ) = (Np + n + <t>) ^ (s) - ( Np+<t>) P^, ( 5 ) 

Put n = 1 

2 P^ ( 5 ) = (Np + 1 + (t>) ^( 5 ) - ( Np+(^) Po (s) 

= (Np + <1) + 1 ) ( NfHi>) P„ (s) - (Np+ <i>) Po ( 5 ) 

=(Np + ^)Po(5) [Np+4] 
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Put n = 2 

3 P^is)^ (Np + 2 + <j>) Pjis) - ( Np+-<|)) P^(s) 

= {Np + (|i+2)^^^2li>lF.(s)-{Np + « (Np+W ?.W 


2! 


-[iVp+«^3/>o(5') 


0<n<Y-l 


(3.3.B) 


Similarly 

n! 

From equation (3.3.9) 

S ^( 5 )- Pn(0) = - (N+ Y-n) X.+n p] jP„(s)+ (N+ Y-n+l)X, P^iis) + (n+1) p, P^iis) 
SP„is)- SP„_i{s) = -[(N+ Y-rt) X,+np] F„(,s')+ ^4- Y-n+l) X P„_i( 5 ) + (n+1) ]i P^iis) 
(n+1) ji P„,,( 5 ) = [(N+ Y-n) X+nn+S] P„(5)- [(N+ Y-n+l)X+s] P„.,(5) 

(n+ 1 ) P„,(s) = l(N+ Y-n) p+n+«l)] F,(s) - [(N+ Y-n+1) p+<i»] F^,(s) 

Put n = Y+l 

(Y+2) Py+2(s)= [(N+Y-Y-l)p+Y+l+(t>] Py+i(s)-(Np+i1») Py(s) 

= [(N-l)p+Y+l+<t>] Py+i(s) -(N p+(|>) Py(s) 


[(N- 


-(N^)<^Po(.) 


. (iVp+^) 

(r+i)» 


r+i 


[(N-l)p4(|>]Po(^) 
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Py+2(s) — 


{Np+</>y (Np+MiN-l)p+S n.-x 
Yi (r+ixr+2) ^ 


Put n = Y+2 

(Y+3) Py+ 3 (s) = [(N-2) p +Y+2+(j)] Py+ 2 (s) - [(N-1) p+(i>] Py+iis) 


= [(N-2) p +Y+2+(|)] 


{Np + (Np + ^}[(.N - 1)/> + ^] 

r!(r+lXF+2) 


Po(s)- [(N-l)p+<|>] 


(Y+iy. 


P,is) 


_ {Np + <f>y (N p+MN-Dp+J ] p . s f(N-2) p+<t)] 
F.(Y + IXF+2) 


Py.3(s) = 


(Np + <t>y (Np + ^)[iN - V)p + ^][(N - 2)p + J . 
7!(r + lXF + 2Xr + 3) ® 


Similaiiy 


P = Pv 


(iVp+^)F (Np+^)„_y 


r*(n-Y) 


F! 


(Y+l)„-. 




{Np+(f>y p^-^N+^x.y 

= p(5), Y^C-1 (3.3.14) 

W + l)„_r 

From equation (3 .3 . 1 0) 

S PM- PM=-[{N+Y- + cp]P, (s)HN+Y-c+ 1);IP,_, (s) + (cp -f (s) 

[cp+a]P,,X«)= [{N + Y-n)j3X + cp + Sl^M-iN-^Y-c + WaAs)-SP,_M 
{S+l)P^M= {(iV-FF-«)^y+<y+?/K(j)-[(i\r+F-C+l)/+tjr]^^(j) 

Where , — = y , — = tf/ 

a a a 

(tf+i) ?„w= ?,(s) - 

/ . v-* + Vp-r 
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r,,, {Np + f)c-Y-\ p/„\ 


(Np+^)c-Y ^ [iNp+f){NiN-\)p+(f).....iN-c+\ + Y)p+^] 


Nn = N(N-l) (N-n+1) 


(Y+1)c.y = (Y+C-Y) (Y+1)c.y.i= C(Y+1)c.y., 




E 

fj- 


/ tl 7 -w ^ ^ 

(iV + F-c + l)— 4-— 
a a 


{Np + 4i)Y {Np + _y 


7! 


c.iY + lX_,.y 


E[(N + Y-C + \)p + <P] (Np + 


iNp + (/>)Y 


a 


Y'. (.Y + lX_r 




a 


{Np^ 

Y\ (7 + lX_, 

Put in (A) we have 

(^ + 1) E,M=\{.N + Y-n)PY+S + iir\ Pois)- 


{Np + ff iNp + <l>)c-Y 


7! 


(Y+l) 


Po{s) 


c-r 


PM- 


' {N+Y-n)pY+y/ 

d+\ 


{Np + <f>Y {Np-{-<l>)^Y j 
71 (7+lU 


(A) 
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From equation (33.1 1) 

SP„(5)- S [(Ar+y-«)y9A+c^+(«-c)a]F„(5)+(N+Y-n+l)pXP„_,(5) 

+[cp,+(n+l-c)a] PnAs) 

[c// + (o + 1 - c)a] P„^^(s)= [(N + Y-n)fiZ+cju+(n- c)a + S'] P„ (s) - 
[(i\^+7-« + l);«^ + 5]P„.,(5) 

[t? +n+l-c] P„+](-s)= [(iV^ + 7-«)>S/ + ^ + n-c + ^tr)]P„(s)-[(/V' + F-« + l)^/ + V')] 
Putn = c+1 

[^ + 2]P,,,(^)= [{N + Y-c-\)PY+S+\+yr\P^Asy [{N +Y -c)fiY+¥\Pc(.s) 
[S+2]P,,^{s)=[iN+Y-c-\)pY+S+\+w]P.M-{il^^Y-c)PY+v\Pcis) 


= [iN + Y-c-\)PY+3 + l^-w] 


{N + Y —c)Py'^¥ 
___ 


(AT^+ji)' (Sp+^s 


V. 


iY+l\-Y 


PM - 


/. (I 

= lN^Y-c)pr^y,\ ft W 


(A^ 4-7-0— + ^~| 

__ _ Ij 


{N+Y-c-\)PY+yPi 


^ + 1 


j7^^(_)_ (iVp+^)’^ {Np+<p\_, [(Ar+7-c)^r+rI(iv+F-c-i)jgr+rFo(^) 


7! 


iY+i)c-r 


(<J4-1X^4-2) 


Put n = c4-2 

XS+ 3 ) Pp,M=^ [(iV4-7-c-2)^r+‘5+c+2-c+F]?c^2(^)- [(iV+7-c-l)^r 

= \(N l Y c 2)3r ■ . 2 ■ rl lN^Y-c)PY+v;gN-^Y-c^l)Pr^^ois) 

>V 7. (Y+l)^ (5+1X54-2) 
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- + Y + 


iNp+if>y \(N+Y-c)fir+V'] Ti... 

Yl (Y+lX_y (<y+l) 


jNp + ^y {Np + <fi),_y [(iVH-r-c)/??M-v|(A^+r-^-t)y?;^y^ p r (jV + F-c-2)jgy + <?+2+vr ^ 
F! (F+l),_, (S+l) (^ + 2) 

p (r)- i^P+^c-Y \N+Y-c)py-¥\i^{N-\-Y-c-i)pr+y^N+Y-c-'tjfir^-xf>\ ji,. 

^ Y. (Y+l)^y (<5+W+2X<y+3) ° 


Similarly 


= ■Pc+<n-c)(‘S) = 



{N + Y-c)+^ 

iPry~^ 

(Np+<^^ (Np+^,_y 

L Py\ 

n-c 


(F+IU 


(-^+1V 


A(s), 


C<n<Y+K 





(Ar+F-c)+-^ 

iPry~‘ 

(Np+<^^ {Np+i^^y 

Pr. 

n-c 


/l(s),C^<Y+K (3.3.15) 


Y (Y+l)^y (S + l)„_, 

The empty system probability Pq, can be found fixim the boundary condition : y] P„ ( 5 ) = 1 


Y+K 


n-^y 


/i=o ii=r 


P)n. 


(r+r)„_y 


-+ 


(NP + <^Y (Np + </>X_y 


im 


(iV+F-c)+-^ 

Pr. 


F! (F + 1) 


c-F n^c 


(^+0„ 


»=*0 


(3.3.16) 


C-1 


¥IV- 


(D Z 


PA- 


w=f 


(F+l)„. 


«— >-w+F 


( a\ 


c-l-Y 


N + 


U (F+l)„n! 


, OA . 

PJn _ I,' 

2^ 1 


l,Y+^;F+l;p 
V P ) 
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n ->«+c 


y+jc 


m” 


m E 




{N+Y-c)+^ , 

^ ^ PrL, 


n-c 

(^+1)r^ 


{firT 

(Ar+r-c)+-^ 

0)„ 

n 


i_d + \)„r& 


= 2 ^ 


UN-^Y-c)+-^;d+\;Py 

\ fir j 


Hence 






Y\ 


iF, 


\ 


UN + Y-c)+^-,S + V,pr 

fir 




h = XP-\s) 


dX 


«! n V >0 ^ 


{Np + (l)f (iV/7 + ^)^y 


r! 


(7+1) 




c-r 


l,(Ar+7-c)+-^;5+l;>?r 


(I) 






HasO 


nl 


5 ^ 


.,(Ar-+-)' 

::^ M pi 




»i 


||1S?1 // 


S («-!)! 


iNp^<fi)c-y 

( 5 ^ + l).-r 

(33.17) 

(33.18) 


(H) 




dX 




= APo( 5 ) 


dX 


{Np + <l>y y {\)„{N + <l>lp)„p" 




TIssO 


(n)a ^o(^) jr 


a;i F! 


5/1 


(Ar 2 + ^)*' 
//’'F! 


A?o(^)F(iVl + 5 ) 

/F. 


r-i 


-iV 


= No ^ 

^ (F-l)I ^ (F-1)! 


r -1 


(n)b AFo(:y) 


5 A 


f (l)„(A^ + ^/p).P" 


= APo( 5 ) 


+ ^ /•j"') ly 7 V + — 1 

^ (F + l)„n! 5A jJ^(F + !)„«! 52 V PJ«. 


(-n (-t.y 

_^ {\UN + ‘l>lp)nP'' |?.f /it 
h (F+i)„("-0! S (r+i)„«i 


h (I'+I)..,"! ^ (F+l).-!! 


. (2).(W+^) ((Ar-l)+^).flo" . (!).(-()■ 

=y_ e. — e. +xY. ^ 

La /V- . lAFV 1 M 


nxsQ 


(r+ixr+2X^ 


(r+i)„«! 


i\)XN^il>lp)„p" 

(F+l)„«! 


^ (V ^ 


(Y + l) S iY + 2)„r2] 


to (y+i)„»! 


,F,(2,(N-l)+^/p;r+2;p)+ A ,F,(l;7 + l;-^/2) 


Hence 


^o(s)~ + J 




(7-1)! 


{Np+^y 

iNp+<f>) 

Yl 

L(7 + 1) 


^^^,FA(^-l)+<^//>;r+2;p)+;i,F,(l;7+l;-<^/^ 


aiD ^o(-y) 




{Np + <lif iNp + 4')c-r 


^o(^) 


(A^P+(^).-r 


A 

ai 


(J^ + l)c-y 

iNp+<li)\ 


jn(r+i),_, 

^F,iYN-^Y-€+iirlpr\6^\;fiY) 


F, iYN + Y-c+yr! fir\5 + \,fir) 


=Poi.s) 

^ois) 

ms) 
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\Np^<h 

(7 + 1) 


+^o(^) 


(iVp+<^)"(iV)c7+<^),_, a 


YliY + l)c-r 


dA 


miN+Y-c+ip/fir^^+Ufir) 


2Fi(l,N + Y-c + ip/fir^S+ 1; Py)Np 


c-y 


(7-1)! 


+ 


{Np + <l,y{Np + 4).-y d ^(})„{N + Y-c-^y^iPY)yPYy 


7!(7 + l), 


dAt^ 


(-!? + l)„n! 


dA 


.i\UN^Y-c+iiflfir\{pYr 


WsxO 


■msiZ 


(^ +!)„»! 

{\)„(N+Y-c^yflfiY)n d 


n=fQ 


(<J + 1)«! 


dA 


ifiry+ms)'Z 


(X)n(firr s 

^iS+l)„rAdA 


iN+Y-c+wlfir)n 


w 

-m 




.f 


J.v 
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- % (sy\\ri 


w 

P^Y. 


( 1 )„ 


^ If f ' 

I 


V 


-f i^)^iN + Y-c + wlpy)„fir .g y-^^ 
tf (^ + 1)„(«-1)! ^ (^ +!)„«! 




«==0 


(5 + !)„,»! 


- (^ + y-‘=+y'^WAy V (2>-(^'^^-‘^-'+<^^'»>->- (^y).+;i,F,(l;g+l;-v>/^) 
(<5 + l) (<^ + 2)„n! 

^ (7^ + 7 - c + v^/^ r ) ^ ^^(2, AA - 1 + 7 - c + + 2;>gr) + 2,Fi(i;<? +1;-^^ /^) 

(<5' + l) 

Hence 






7!(7+l) 


c-r 


Pr 


Po(s} 


V +Vc--}' 


(Np^i^y 

(Y-l)l 


Y-\ 


+ 


{Np+<l>y {Np + ^X_y {N + Y-c + y/l Py) 


n(7 + !),_, 


5 +\ 


:,F,C 1 ,N -\ + Y -c-^y / 1 pr,S + 2 -,py) 




7! (7+1) 


c~7 


Hence 


l=Po(^)LE 


r-i 


{Np+<t>y~'.r {Np + (l>f~' Xn-KT j., t ^ 
^ +Np ^F,(l,jy+(^//7;7+l;/>) 


/j=l 


(n-1)! 


(7-1)! 


+i y^ + <!^)'"' ^F,i2,{N-l)+<^/ p;Y^2;p)+ &^^^r A,F,(l,Y+l;-;-^/A) 


(F + 1)! 


F. 


O' + l)^, O'-!)! fh 


+ 


{Np+^)^ {Np + f),_Y {N + Y-c + ^l fir) + 


n (r+ii_y 


s+i 


JL 


Pl£±iL^lP±£^^Ffi-,5*i,-^ 

F! (F + l)^_y 


,5 +2-, fir) 

(3,3.19) 


3.4 SPECIAL CASES 


/^(0 = -NX.Po(t)+^Pi(t) n=0 

?;(/) = -(N^+nn) Pn(t) + N P„-i(t) + (n+1) |a P„+i(t), l<n<Y (3.4.2) 

p^(t) = -[(N+Y-n) X+nji] Pn(t) + (N+Y-n+1) A,Pn-i (t) + (n+1) pPn+i(t) 

Y+l<n<c (3.43) 

p^{t)— -[(N+Y-n) pX-+Cfx] Pc(t) + (N+Y-c+1) A-Pc-i(t) + (cjx+a) P(^i(t) , n— c 

(3.4.4) 

PXt) = - [(N+Y-n) PA-+CH+ (n-c) a] Pn(t) + (N+Y-n+1) pA, Pn.i(t) 

+[ Cjx+(n+l-c)a] Pn+i(t) , c+l^<Y+K (3.4.5) 
P’^^(t)= - [CH+ (Y+K-c) a] PY+k(t) + [N-K+1] Py+k-i (t), n=Y+K (3.4.6) 
Let a =0, p =1 and N=K then the above equations are :- 
i^(0=-iaPo(t)+nPi(t) n=0 (3.4.7) 

p:(t) = -{KX+n\x) Pn(t) + KX P„-i(t) + (n+1) 1^<Y (3.4.8) 

p;(^) = .[(K+Y-n) ;.+n^i] Pn(t) + (K+Y-n+1) APn-i (t) + (n+1) pPtH-i(t) 

Y+l<n<C (3.4.9) 

P;(0= -[(K+Y-n) PA.+CH] Pc(t) + (K+Y-c+1) APo-i(t) + c\i P<H-i(t) , n= c 

(3.4.10) 

P^t) = - [(K+Y-n) pA,+cti] Pn(t) + (K+Y-n+1) pA. PB-i(t) + C\i Pn+i(t) , 


122 


c-V'l^<Y+K 


(3.4.11) 


f;,(0= -cnPY+k(t) + (N-K+l) X Py+k-i (t), n=Y+K 
(3.4.12) 

From equation (3.4.7) 

/^(O =-iaPo(t) + qPi(t) n=0 

Taking laplace transform 
S Po (s) - P(o) = - K.X Po (s) + M- Pi 
\iPi(s)= Po(5)[KX+s] 

^ (s) = (Kp+ ((>) Po is) where p = , <!>= S/p 

From equation (3.4.8) 

/>;(() = - (ia+ nn ) P„(t) + K X P„-i(t) + (n+l) H P„+i(t), l<n<Y 

Taking laplace transform 

P„ (s) - Pn(0) = - (KX+ np) P„ is) + KX P^, is) + (n+ 1 ) p P„+, is) 
(n+l) n?„,(s)= (KX+ nn+s) PM- KXP.-,M- Pn(0) 
(n+l)nP„,(s)=(K>.+ nn+s) ?.(s)-(ia+s) 

(n+l) ?„,(J)= ( K(f+n+i(i) P.(s)- ( Kfrii)) P^,(s) 


m 


s 


Putn=l 


2 Pj (5) = ( Kp+ l+(j)) P, is) - ( Kp+<1)) Po (s) 

=(Kprfl+<l)) (Kp-t-(|)) Po(5)-( Kp+<|>)Po(s) 
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2 F2(5)=(Kp+<|)) Po(5) [J^P+<1>+1-1] 

Put n =2 

3P3(5)=(Kp+2+<|>) :P3(s)-(Kp+<j>):^(5) 

= ( Kp+2+<l>) P,is )- ( Kp+(t>X Kp+<1>) Pois) 


= (KP+^L [Kp+(|)+2-2] P,(.s) 


Similarly, 


F,W 0<nsY-l 

flU 


(3.4.13) 


From equation (3.4.9) 

P^(t) = -[(K+Y-n) X+ np] Pn(t) + (K+Y-n+1) X,Pn.i(t) + (n+1) p Pnf i(t) , 

Y+l<n<c 

Taking laplace transform 

S P„(sy Pn(0) = -[(K+Y-n) X+ np] P„is)+ ( K+Y-n+1) X P„.As) + (n +1) P 

(n+1) p P^M^ [(K+Y-n)>. + np +s] ]p„(s)- (K+Y-n+1) X P^M- Pn(0) 
(n+1) p F„.,(^)= [(K+Y-n) I+n p+s] F„(«)-[(K+Y-n+l) X +s] P„_,is) 


Putn=Y+l 


(Y+2) n Pr,,is)=[(K+Y-Y-\) X+ (Y+1) +S ] P,„(s)-[ (K+Y-Y-1+1) X +S] 

P ^+1-1 

(Y+2) ?„(*)= [ (K-l) p + (Y+1) + « ( Kp +.|,) P,(s) 


(Y+2) F,„is)= [(K-l) p +Y+1+ ( Kp +« P,(s) 


_ (Kp+^) 


Y+1 


a + l)! 


■[(K-l)p + Y+l+(|)-Y-l] P,(s) 


P ' ^- (Kp+(l>y (.Kp+mK-\)p+<l>) 

F4-2n } xrt /xr . tX/Tr . 


F! 


(F + lXF+2) 


Put n = Y + 2 

(Y+3) P,,,(s)= [ ( K+Y-Y-2) p+Y+2+« ?,.,(5)-[ (K+Y-Y-2+1) p + 
<HP™(s) 

= [ (K-2) p + Y +2 + 4,] ?„(s) - [(K-l) p + 4,] ?„(*) 


[ (K-2) p + Y+2 + 4. ] ^ 

TJ y, (y+txr+2) " 


-[(K-l)p + 4,]^&^?.(^) 


(.Kp+^Y (Kp+^)((K-l)p+<») ^ 


n 


(F+iXr+2) 


Po (5) [ (k-2) p + Y+2 + (|) - Y - 2 ] 
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n (7+ixr+2XF+3) " 


p^p - + iKpH)^r P 

‘^y+Cn-r) — -v -M) 


Yl (Y+l)„_ 


is) 


Y < n < c -1 


iKp-^ 


Yl 


p-'(X+^)„, _ 

^ - Po(s) 


iY + l)„- 


(3.4.14) 


From equation (3.4.10) 

pUO = - [( K+Y-n) X + c\x] PXt) + (K+Y-c+1) X Pc-i(t) +cn Pefi(t) n=c 
Taking laplace transfonn 

SPX5)-Pc(G) = -[(K+Y-n)X + cp] F,(^)+ ( K+Y-c+1) >. F,_,(^)+ cpP.^X^) 
cuF,X^)= [ (K+Y-n) X+ cp+ s ) ^(;s) - (K+Y-c+1) ^F,^(^)- Pc(0) 
cF„(^)=[ (K+Y-n) p + c + (|)) F(5) -[(K+Y-c+1) p+<l>] F_, (5) 
where p = yp , <t> = s/p , 

c^.,W= [ (K+Y<) P++ + C] 


- [(K+Y-c+1) p + 


(A) 


=> [(K+Y-c+1) p + <|>] 


jKp + S'^ iKp + ^)^r-i p(„. 

Yl (F+Dp-r-. 


(i:+7-c+i)-+-^ 


Pi 


{Kp+4^y c{Kp + ^)c-\-Y p^(^s) 


iY + l\.r 


(Kp+^y 

Y\ 


[(K+Y-c+1) p + m^P+^U-r 777^ 

(I -rijc-y 




n 


(r+iy.y 


^o(s) 


Putting in (A) we get 


C pAs) 


....... . . , . . (Kp + (l>y iKp+^X.y -5y^^_ 

[(K+Y-c) p + (j)+c] (Y+\X ~ 


(Kp+</>y iKp+<l'X_y 


F. (r+ 1 ) 


(5) 


c-r 


C ?„,W= [(K+Y-c) p + i|.+c-cl (Kp+ji)' J - W 


^c+{ (■^') 


_ \{K^Y-c)p + <y\ {Kp + <^y iKp + ^),_y 

c F iY+\X.y ” 


From equation (3.4.1 1) 

P'„ (0 = ” [( K+Y-n) X- + cp ] Pn(t) +<N+Y-n+l) X P„.i(t) + cp P„-n(t) 
SP„(5)-Pn(0) = -[(K+Y-n)X + cp] P„(5)+(K+Y-n+l)X P^i(5)+cp P^,(s) 
cp P„+i(s)= ( K+Y-n) X+ cp+s > P„is) - (K+Y-n+1) X P„_i(5r)-Pn(0) 
cp P^,(s)= (K+Y-n) X+cp + s) P„(s)-(K+Y-n+l)X P^,(5)- s P„_,(s) 
cp P„,,( 5 )= ( K+Y-n) X+ cp +s) PM- [( K+Y-n+1) X+s] P^,(^) 


c p„^,(s)= (K+Y-n) p + c + (|)] P„is )- [( K+Y-n+1) p 


c + 1 <n<Y+k 


Put n = c +1 


C [ ( K+ Y-c-1) p+c+<l)] [( K+Y-c) p+(f>] P^is) 

= [ ( K+Y-c- 1) p+c+^] * 

[( K+Y-c) P+<|,] 

(F + l),_y " r. (F+l)c_r 


[( K+Y-c) p+(|)] 

Yl (r-hi),_Y 


iK + Y-c-l)p + c+^ 


[ ( K+Y-c) p+« 

Y] (F + I)p_y C 


J, _ (Kp + ^y (Kp+^X.r [(K + Y-c)p + ^l(K+Y-c-])p+Wj p , . 

^c+ 2 W 2 -^OV*/ 


F! (F+l), 


c-Y 


Put n = c+2 


C P,^ 3 is)= [ ( K+Y-c-2) p+c+(l>] P^+2(5)- [ ( K+Y-c-2+1) p+<l>] 


[ ( K+Y-c-2) p+c+(H [ ( K+Y-c) p+w 


Ik + Y-c~\)p+ 


/>,w 


- [ ( K+Y-c-i) P++] ^ 


FI (F+1) 


is) 


c-r 
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iKp+tl>Y {Kp + <t>\_y \{k + Y-c)p+^\{K + Y-c~\)p+(I>] s 


Yl 


(Y + l) 


Pois) 


C-Y 


[ ( K+Y-c-2) p+c+(}>-c] 


?„.(*) = K K+Y-c) p+«[( K+Y-c-1) P+IH 


\(K + Y-c-2)p+^] 


^-Pois) 


Similarly 


p = p 

* n •* c4-(«-c) 


XKp^■^Y (Kp+^\.r [(.K + Y-c)p+,p\^ 

Y\ (F+l),_, c"-' 


j,_ (Kp + ^y iKp+^X_y [(K + Y-c)p+^^_, ^ 

” n (r+i),_, 



f 

kJ- 

p- 

1 

+ 

0 

1 

-f 

I 


(Kp+^y 

K PJ 

c-Y 

L P\ 

n—c 


n (y+iu 

(Kp+^y XK +<f>lpy_y [{K + Y- c)+^/pl_ 


r. 


(Y+1)c.y 


p"-^ Pois) c<n<Y+K (3.4.15) 


Hence, when a=0, P =1, and N=K then M/M/C/K/K with spares only then 




n\ 


Pois) 


" ' n (F+i),., 


L (J^p+<Z>/ iKp+<^X_r [(K-^Y-c)p+ 
Y\ (r + i),_, 


Pois) 


to «'• tr F! 


p A: + 


PJ„. 


(F + l)„- 


- + 


(Kpt^y 'f \K*Y-c*4,lp\„ p„-Y^ 


F! (F + 1), 


c-Y 


F-1 


/i»0 


{Kp+^r ^iKp±^) 


Y c -1 


n\ 


n 


m-Y 


k P)n-Y ^n-r 

(F + 1). 


w 


(Kp+^y 


iK + ^X_y 


r+K 


FI (F + 1), 


C-K «*=C 


'K + Y-C + 


.tt-C 


-y (Kp+<^y ^ (Kp+^y . iK+^/py„ipy ^ 

h nl F! S (F + l)„n! 


iK + Y-c + <f,/p)„p''il\ 


FI (F + 1), 


'<?-F 


/IKSO 


cV 


ISO 


0<n<Y-l 

Y<n<c-1 

c^<Y+K 

(3.4.16) 

(3.4.17) 


(5^^+«e^,F,(lJC+*/p;Y+l;p) + 
^ «! n 


2Fo(1,K+Y-c+Vp; - ; (p/c)) 

n (r + i)^_Y 


(3.4.18) 


L= AP,(j)^ 

OA 


j. (^+S%(f:p+^) 2Fi(lJC+(,/p;Y+l;p) 


w =0 


m 


+ J&+Y-c+(l)/p; - ; p/c)] 


(3.4.19) 


(I) 

5/1 


«! 5yl S "i 


^.wl ^ r-- ' - 

S n! 


= wl 


(II) A?,{5) jF,(I Jt+,ti/p; Y+1; p ) 

OA 1 1 — 


..i?,(42Fi(i.k+Vp: y+1; p ) ^ . + iKe il t A ^f,(i jk+Vp; y+ 1; p ) 


rma -, 5 (^p+(^) _ , 5 // _ // i:_ (Ap+» 

^ ^ '^51 Yl ^dA n 71 p (7-1)1 ^ 


ai)b 4^ 2F,(1 JK:+4 ,/p; y+ 1; p) I p- 

a;i aA ^ (7+i)„/j! 
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= X, Y n (—y-^ i. H-X, V 0)« P” fitLy 

h, (j + O^"! n n ii (V + l)„n'. ^ 

“ (^)w(-^ ^ ! P^ n “ n) p" A 

“§ (1^ + 1).(«-1)! ''""^S (F^'^l’" 

n=sl !•* 


iK + <t>tp)p^ (2)„(K-^+^lp)nP'' 


(r+i) ^ 


(F + 2)„«! 


+?. ,F,(l;Y+l;-«j>/X) 


(F + l) 


2F,(2^-l+<i»/p; Y+2; p) + X ,F,(1 ;Y+1 ; -^X) 


Hence 


mis)-^ — - ^ - - -2Fi(lJC-H)/p; Y+1; p) = 
/! 


Kp. (f^tC'‘ .2F,(l,K+Vp; Y+1 ; p) + <^^2F,(2JfC-l+.|)/p; Y+2; p) + 


(Kp+^y 


■X ,F,(1;Y+1 ;-<!>/>-) 


^ 2Fi(l.K+<|>/p; Y+1; p)+ ^^^t^2F,(2JK:-l+4»/p; Y+2; p) + 


x iEP^} . ,F,(I;Y+1;-<|>/X) 


(III) >W> is)^ 2 Fo (1 ^+Yhc+,|>/p;-; p/c) 

OA XI VJr 


1 ^ (A /7 + ^) 


(r + l)._K 


■2Fo(1,K+Y-<;+Vp;-; p/c) A 


SA FI 
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z — 2Fo( 1 JC+Y-c+fp;-; p/c) 


F! (F + 1) 


(in)a A — 2 Fo( 1 JfC+Y-c+(f)/p;-; p/c) 

oA 




i^)„H>ipf(p 


n\ \c 


a)„ (-<i> 


n\ \ c 


= A. ,Fo(l;-;-<l»/c) 


Hence 




(Y-1)! 


+ t<^p,F,(2,K-l+ <i>/p; Y+2; p) + X ,F,(1;Y+1; - ^ + 

2Fo(1 ,K+Y-c+ (t»/p; - ; p/c) 

(F + l),_r 


iKp+f l_ {Kp + ^X^^ ,Fo(l;-;-«t»/cll 


(3A20) 


F! iY + \X.r 

The machine availability (rate of production per machine) is M.A. — 1 - L / K 
The operative efficiency (utilization) is O.E, = 




I* 



INTRODUCTION : 


■f 


Queiieing models provide the basic fiamework for efficient design 
aiwi analysis of many practical situations. Along with several other 
assumptions it is customary to consider that the arrival and service 
processes are independent A queueing model in which arrivals and 
services are correlated is known as interdependent queueing model. 

In interdependent queueing model, Lindley (1952) [104] analyzed a 
single server queueing model with recurrent input and arbitrary service 
times with diat assumptions the service times of the n* customer and time 
interval between arrival epoch are related by a linear function. 

The present chapter discusses a queue model of the type M/M/C 
where, queueing system with the assunqjtion that the arrival and service 
processes of the system are correlated and follow a bivariate poisson 
distribution. In addition to this interdependence whenever the queue size 
reaches a certain prescribed limit R, the arrival rate reduces from A<o to Xi 
and it continues with the reduced rate X.i as long as the number of customer 
in the queue is greater than some other prescribed integer r (0<r^) and 
when the system reaches r, the arrival rate chaises back to and the same 
process is repeated. This operating strategy is called controllable arrival 
rates with prescribed integer R and r. 


1^4 


In general waiting line of a queueing model increases due to slow 
serving rate or fast arriving rate. There are several research papers in which 
the service rates are controlled to decrease the queue length. Queueing 
models with controllable service rate have been discussed by Gray W., 
Wang P Scott, M and C.M. Wang P. Many years ago, ConoUy and Hadidi 
(1969) [38] have studied the single server queueing model with that 
assumption the rate / T„ is constant for all n, where Sn is the n*^ customer 
service time and in the time interval between n* and (n+1)* arrival 
epoch. Mitchel and Paulson (1979) [114] have considered and MM/1 
queueing system with an assumption that the interarrival times separating 
the arrival fiom that of his predecessor and the customer service times 
follows a bivariate exponential distribution. Their analysis is based on 
simulation of the waiting time process. ConoUy and Choo (1979) [37] have 
provided means for exact calculation of the waiting time for a generalized 
correlated queueing with exponential demand and service. Lenganis 
(1987) [103] , Borst and Combe (1992) [25] have studied the busy period 
analysis of a correlated queue with exponential demand and service M.I. 
AftabBegum & D. Maheswari (2002) [107] have studied the M/M/C 
interdependent queueing model with controUable arrival rates. Here we 
consider the arrival and service processes of the sii^e server infimte 
capacity queueing system are correlated and follows a bivariate poisson 
process having the joint probability mass fimction of the form 


p(Xi=Xi, X2= X2;t) = e 


mm{ 2 ) 




\itin - 




£ (ef)-' [(/l,- - - J)!(X, - J)\ 


J=0 


with parameters ?ii (I =0,1), \x as mean arrival rate, mean service rate and 
mean dependence rate (covariance between the arrival and service 
processes) respectively. The operating strategy of the system is controllable 
arrival rates with prescribed integers R and r. 




Transient State equation 


/{.'(O) - - (io-e)P„(0) + (/i-e) P.(O) (4.2.1) 

P„'(0) = -(2-o+«/^-(« + Oe)P„(0) + (^-e)P„-.(0) + (n + l)(//-e)P„,,(0) 

1 < n < c-1 (4.2.2) 

Pn'(O) = - {X, + cju-{c+\)e)V„{0) + {2.o-e)P„_i(0) + c (//-e) P„,,(0) 

c < n < r-1 (4.2.3) 

p;(0) = - iA,+cju-{c + l)e) P,(0) + (A, -e) P,.,(0) + c (//-e) P,,i(0) 

+ c(//-e)P,,,(l) (4.2.4) 

W) = - (2„+c/i-(c+])e)P„(0) + (io-^)Pn-.(O) + c(/i-e)P„,,(0) 

r + 1 < n < R-2 (4.2.5) 

P,J ,(0) = - (/lo + t.';i-(c + l)e) Pr.,(0) -I- (A,, -£?) Pr.2(0) (4.2.6) 

P 4 i(l) = - (^i+t^/^-(c+l)e)P,+,(l) + c(//-e)P,+ 2 (l) [4.2.7) 

P„'(I) = - (A,+c//-(c + l)e)P„(l) + c(//-e)P„,,(l) + (A 2 -e)P„.,(i) r+2 < n < R-1 ( 42 .g) 
P' ( 1 ) = . ( 1 , +c/j-ic + l)e) Pj,(l) + (A, -e) P„.,(l) +iA^ -e) Pr.,(0)+c (/^-e) Pr^i(1) (4 2 9 ) 

P;(l) = - (2,+c/^-(c+l)e)P„(l) + (/l,-e)P„.,(l) +c(//-e)P„^i(l) n < R+1 (4.2.10) 

Define boundaiy Condition P„(0) = SPn-i (0) (4.2.11) 

From equation (4.2. 1) 

P„'(0) = - (2„-e)Po(0) + (//-e) Pj(O) 

TakingLaplaceTransform- 
S Po(0) - P(0) = - (A^-e) Po(0) + (/i-e)Pi(O) 


(//-e) Pi(0) - (/lo-e+5) Po(0) 


A(0) 


A(0) 




V n-e n-e) 
^X^-e , 5 ^ 


Po(0) 


Po(0) 


V//-e ^x-e) 

From equation (4.2.2) 

i^(0) = - (2o+«//-(«+l)e)P„(0) + (>^o"^)Pn-i(0) 
+ (n+l)(//-e)P^i(0) 1 < n < c-1 


Taking Laplace Tarnsform - 

S P«(0) - P„(0) = - (Xo +n^-{n + l)e) P«(0) + (Xq -b) P,,-i(0) + 

(n + 1) (// - e) P «+i (0) 

(Ao +n/i-(« + l)e + 5) P«(0) = (2o-e) P«-i(0) + (n+l) (;i-e) P„+i(0) 

+ sP„-i(0) 

(n + 1) (//-e) P«+i(0) = [2o+«yU-(« + l)e+5] Pn(0) - (^Q-e) P;h-i(0)- sP„-i( 0) 
(w+l)Pn+i(0) = - (Ao-e)/(//-e) Pn-i(0) 

(n+l)P„+i(0) = [2o+«//-(«+l)e+«s]A/^-e)P„(0) - (2o-e+^)/(/^-^) P"-»(0) 

put n = \ 

2P2{0)=^[Xo+p-2e+s]/(jJ-e)P\iO)-iXo-e+s)lip-e)PoiO) 


2P2(0) = 

2 P 2 { 0 ) = 


Xn-e + p-e+s 


p-e 




Xp-e+s 


L M-e J 


Po(0)- 


I J I J 


*.11 
-+- -+ 1-1 

e p~e 


p-e p-e 


Po(0) 
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put n-2 
3^3(0): 

3P3(0): 


Zq +2p-3e + s 
M~e 


P2{0)- 


ZQ~-e + 2(p--e) + s 
p~e 


Zq -e + s 


Z.-e s 
+- 




iM-e fi-e 


P,iO)- 


A„-e ^ s 


H-e n-e) 


/l„-e ■ s 


fx—e //-c 


Po(0) 


3^3(0)= 2 


P,(0)-- 


^0 ^ J ^ 

p-e p-e 




^o(O) 


^ S 

0 +2 + 2 


H-e n-e 


V.\fi-e fi-e 


PoiO) 


Similarly P„(0) ■ 


n! 


Z„-e . s 


ln-e ^i-e\ 


Fo(0), 0<«<c-l (4.2.12) 

A-ej 

From equation ( 42 . 3 ) 

i^;( 0 ) = -(2o + c// - (C + l)e) ?„ ( 0 ) + ( 2 o - e) P„_, ( 0 ) + c(/2 - e) ( 0 ), c < « < r - 1 
Taking Laplace Transfrom 

sJn ( 0 ) - P„ ( 0 ) = -( 2 o +cp-ic + l)e)Pn ( 0 ) + ( 2 o - e)Pn-i ( 0 ) + dp - e) P«+i ( 0 ) 

c(/i - e) P„« (0) = (2o + c// - (c + l)e + s)P„ (0) - (Aq - e + 5) Pn-» (0) 


C P n+\ (0) — 

cP„+i(0) = 
put n = c 
cPc+i(0) = 

cPc.i(0) = 


^0-^ , c(/2-e) , 5 

' + ^ r • ■ ■■ 


ju-e ju-e /u-e 


P"(0) 


Aa c 


- + ■ 


I 


- + C + - 


p-e) 


Pdoy 


'Ao-e 


\P-e p-ej 

\ 

s 


P»-i(0) 


^Q-g 

V M-e 

'hzl 

^p-e 


- + C + - 


- + C + - 


^0 ^ 1 
'/^-e //-e 


//-e 
" ^ 


Pc(0) 


p-e p-e 


An -C 5 

_o + 


p.-j(0) 


p-e p-e 


Pc-l(O) 


\_ 

c\ 

1 


An -e 


V_ 


- + - 


yp-e p-e) 


Ao-e, 5 


cPc*,(0) 


C?c.l(0): 


1 f A„-e 


(c-1)! 


Po(0) 

V-' 


-+■ 


cl I p- e p-e 


Po(0) 


p-e p-e 

Aq — c 


Po(0) 


-+c- 


-c 


L^-c yU-e 


^0 ^ I -f_ 

//-e //-e 


?»(0) 


P..i(0) = - 


An-e s 

0 . + - 


\C+1 


cc\\p-e p-e 


PoiO) 


139 


put n = c+l 


cPc+2(0) = 


f 


X.-e s 

_o — -+c-^ 


V 


p-e 


r 


Pc+l(0) 


Xo ^ 


Xn-e , s 

—0 + C + 


1 I An -e 


{p-e p-e 

V c+1 


Pc(0) 




//-e j c c!\^ //-e 

p.(0) 


/>.(») 


o 

1 

1 [ 

1 

O 

V 

s ^ 

h 

1 

H 

1 

c\\ 

^p-e 

p-e) 


1 I /l„-e 


A 


c+l 


cc \\ ju-e 


Po(0)l 


An-e S 

_Q c 


p-e 


p-e 


Xo-e 




c+2 


ccl \ fi-e ju-e 


Po(0) 


Pc+2(0) 


1 (x,-e 


\ C+ 2 


c^cV 


p-e) 


Similarly Pn(0) = Pc+(n-e)(0) - 


Po(0) 


An-e 


Po(0) 




c < « ^ r (4.2.13) 


\p-e p-e) c 

From equation {A2A) 

f '(0) = -( Aq + c /2 - (c + \)e)P, (0) + ( Ao - c) (0) + c (/2 - e) P„, (0) + c (/2 - c) (1) 

Tafeng- Laplace Transform 

SPri0)-P^{0) = -(Ag +c/i- (c + l)e)i’r(0)+(2o -e)Pr-i{0)+c{p-e)PrH{0) +c{p-e)Pr*iil) 
c{p - ej)Pr*^ (0) + (1)] = [a 0 + c /2 - (c + l)e + ^ JPr(O) - (Ag - C + S)Pr-l(O) 

An +c/2-(c + l)e + 5 

c(/ 2 -e) 


P„l(0)+Pr.,(l) = 


Ag-e + c(/2-e) + .5 

I c(/2-e) 


Ag -e ^ .y 


[/2-e /2-c 


w r ^0-^ . ^ - 




Po(0) 




pi-e fi-e 


\Xn-e . s 


ir+l ^ 


e'-'c! (c(/ 2 -e) c{p-e) 
1 


An ■” C . ■S 


/ 2 -e /i 


5 r' Pn(0) 
-ej c'"’"' c! 


[/ 2 -e / 2 -e 

P,*i(0) = 


c'-‘*'c\ 

nr+l 


i2Z£+_i-.,j 

/ 2 -e / 2 -eJ 


Pr*. (1) 
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From equation (4.2^) 

P;(0) = -(Ao + c//-(c+l)e)P„(0) + ao-e)P„_,(0) + c(//-e)P„,,(0) r + \<n<R-2 
Taking Laplace Transform 

SPn (0) - P„ (0) = -(/lo + c// - (c + \)e)Pn (0) 4- ( Ao - e)P«-i (0) + c(// - e) (0) 

c{p - e)P«+i (0) = [Zq +cp-(c + l)e + 5 ]Pn( 0 ) - (Zq -e + s)P n-i (0) 

putn=r+l 

c(ju-e)Pr+2(0) =[\ +c//-(c+l)e+>s']Pr+i(0)-(.^ -e+5)Pr(0) 

, s T' ^.(0) 


— +CJU— (c + l)e+.S'j 


ju~e id-e J 






Z,-e^ s Xm 

\^jU-e }j,-eJ d~^c\ 


(^-e 



[p-e 

i 

L 

Mc\^ 


\ +C/i-(c+l>+5-c(//-e)] 


-Pr+i(l)[;io +Cfj.-ic+\)e+^ 


■ + ■ 


C(//-e)Pr+2(0) 

" p-e p-e j 
- Pr+] (1 ) [Aq - e + c{p - e) + 5] 


Z(\ 6 


,r+l 


Po(0) 


[(io -e + c(p-e) + s- c{p - e)] 


cPr.2(0) 


p-e p-e ) 




PM 


Pr+t(l) 


In -e 


•+■ 


Ip-e Z-ej 


Zn~e s 

-J— + C+— 


p-e 


P..2(0) = 
Pr-.2(0) = 


Zc^-e 


\ 


r+2 




p-e p-e) 


An -e 


\ 


r+2 


■ + ■ 


p-e p-e) 


PM 

M'^cl 

PM 


Pr+x(S) 


-Pr+l(0) 


p-e 
\-e 


c{p-e) c{p-e) 


+ 1 


1 + 


An e 




c{p-e) c{p-e) 


Pr.2(0) = 




. r+2 


- + ■ 


jU-e l^-e) 


PoiO) 


PrAD 


C Cl 


l + B 


l-B 


(l-B) 


T n ~ ^ S 

where B = — h— r 

c(ji-e) c(//-e) 




Pr.2(0) = 


A-n C 


. -|- . 


\^jU-e fi-e 


ML 


-PrH-l(l) 


1-P 


Pm? n = r + 2 

c(/i — e) Pr +3 (0) = \Aq + cju — (c + V)e + s\P r+z(0) — (Aq —e + s) Pr+i (0) 
Pr+3(0) = 


Aq 4-c//-(c + l)e + 5'1 

fA^-e^ s ^ 

P (O') — ( 

ciM -e) J 

\fi-e ^-ej 



r 


Aq — 6 -h S 

c(//-e) 


A 


Ao-e 


fj.-e fi-e 


r+l 


n(0) 

Mc\ 


Pr.x{\) 


o 

1 

^ 1 

n(o) 

~ AQ-e + c{p-e) + s ^ 

[p-e 

p-ej 


c{p -e 


+ Pr-.X{\) 


A.-e + s A(,-e^c{^-e) + s 


c{fi-e) 


, An —e-\-s 

1^0 


A 


c{id-e) \ c{fi-e )\ 



1 


Aq - e + c(//-e) + 5 ^ 

1 

1 

p-ej 

r+2-c 

c c! 

c{p-e) 


+ 


Pr+l(l) 


Aq—s + s _ j ^ 
c(^-e) c(/j-e J 




Aq ^ 


s r Po(o) 


ju-e M-e 

P.*i(1)[b-(1 + 5)"] 




:k_£+i+^ — 1 


//-e 


p-e 


+ 



Pr.3(0) = 


A^-e s ' 


■ 4 * ' 


\ ju-e ju-e 


+ Pr.^ (1) [b-\-B~ -2b] 


f 1 ^'■+^ 

An -e s 
— + 


[p-e fi-e) 

^C\ ^ 


c c\ 


Pr.^iO) = 

Pr.3(0) = 

Similarly 


'1Z^ + ^ 


\/a-e n-e) 


2{\ & 


'i^o — c> S 
J. 


y.3 


1^/i-e ^-ej 


^-?.,(i) 

cl 

Po(0) y, 


'(l-B)(\ + B^ +B) 


l-B 




r+3~c t 

C Cl 


l-B^ 

l-B 


Pr,(0) = Pr^n-ri0) 


A^fx c 






where B = 


A^rx ^ 


y/u-e fj,-e 
s 


Po(0) 


- + • 


c{p-e) c(//-e) 


c”~‘cl 
r + 1 <n<P-l 


Pr.l(l) 


/ 1 _ 5»-r ^ 

l-B 


(4.2.14) 


From equation {A. 2. 6) 

P R-\ ( 0 ) = — (Aq +Cju — (c + l)e) P^_, ( 0 ) + {Aq — e)P ^_2 ( 0 ) 

Taking Laplace Transform - 

S Pr-\ ( 0 ) - P;j_, ( 0 ) =-(Aq+cp-(c + l)e) Pr-\ ( 0 ) + (Aq - e) Pr-i ( 0 ) 
— {Aq + c// - (c + l)e + 5 ) Pr-i ( 0 ) + {Aq - e + 5 ) P/f -2 ( 0 ) = 0 


-(2o +c^-(c+l)e+5) 


+(2o-e+-s) 


Arx-e s 
— — + — 


Ap-e ^ s 

KM-e M-ej 

s«-2 


_P ,m 


'"1-5''“’"'''^ 


V 


\-B 


yn-e n-e) 
f 


“ H- cfj , — (c + l}e + s) 


A) £ 


p ,m 

s«-I 




1-5 


= 0 


yp-e p-ej 


Po(0) 


+{A) -e + s) 


^-g | -y 


x/e-2 


ji-e p-e 


Po(0) 


Pr^lO) 


(Xq +cju-(c + l)e + s) 


1-5 


- (/Lq —e + s) 


1-5 


= 0 


X. -e s 
_o + 




^ju-e jj-ej 


n(Q) 


[-Xq -c(^-e) + e-s + ciju -e) 


+ P.+,(1) 


{(^ - e) + c(ju - e) + 5} 


1-5 


-(Aq -e + 5) 


1-5 


R-2-r 


1-5 


= 0 


' Aq - e s ' 
_u + 


//-e 
1 


+5 .ni 




1 + 


c(M-e 
Xq — s + s 


1-5' 


1-5 


1-5' 


1-5 


= 0 


5rH-l(l) 

5.+i(l) 

5.+i(l) 


1 + 

V 5 

^ 1 + 5 
5 


^ Z' DR-2~r^ 




1-5 


(l-g«-.-r) 

1-5 


Aft — e 5 
_2 + 


. «-i 


fi-e fi-e) 


PpiSi) 

c'*-'-" d 




^ 1 ^ 


fAo-e , 5 ) 

Po(0) 

1-5 J 


1 1-5 J 


I 

! ^ 

r - 

! 

c c\ 


V -D yv 

(1 + 5)(1 - 5 ) - 5(1 - 5^^-^-^ ) 

5(1-5) 


Aft -e 5 

J2 + 




J KP-e 5-e. 


PoiO) 


4-5-5 


5r.l(l) 

Pr.l (1) = 

Pr.l(l) = 
5r.l(l) = 

5r.l(l) = 

5r+l(l) = 

5r.,(l) = 


R^r 


i5 + 




1-J5 


/?-r 


5(1-5) 


( Xq -e 


\ 


5(1-B) 

f An - e S 

\M-e M-^J 


Aq ^ 


A 




J V 


ju-^ J 


PoiO) 


/?-i 


PoiO) 


s 


(2 ^ n 


Ao-e 


+ ■ 


\yL-e jLi~e 
(X. 


J 


• + • 


\^i-e fj.-e 


c c: 

fo(0) g(l-g) 

e'-'-'c! (l-B'-') 

fo(0) c(;<-e) g(l-g) 

B„(0) 1 B{l-S) 


An -e 


+ ■ 


c''"V!5 (1-5''-'') 
5o(0) 5'(l-5) 


//-ej c'‘-V!5'(l-5'‘-') 

5o(0) 5' (1-5) 
G^“"d (5'’ - 5^) 


Aq € 


V 


\fi-e fi-e) 


Xq e ^ s 

+ 




^fi-e fi-e) 


c c\ 
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Pr+l(l) = 


/In - e ^ 


c{fi - e) c{fj. -e)J cl 
(1 - B) [b'' - B^y 


B’^ (l-5)[5'^ 


c\ 


Q 

= B^^’’ {\-B)~[b'' -B^^ y Fo(0) 
cl 


PrH (1) = B''"(1-B)^^P„(0) 
where A = \b'^ - ]~ 


cl 


(4.2.15) 


Substituting Equation (4.2.12) in Equation (4.2.1 1) then we get. 

'>o(0) 

c! 


n(0) = 


{ X^-e s 


\ /u-e ju-ej c" cl 




V 


1-B 


n(0) = 


Xq e 


Y 


c{fi-e) c{/j~-e) 


d d 


l-B 


= — AP, (0) 

d cl 


l-B' 


\-B 




C ^R^rC _jgRylp^^Q>^^Y-B^-'') 


cl 


cl 


= B'B"-'PM)— - B^*' -£''r'^’o(0)(l - B"-') 


cl cl 


c 

c! 


B'^*'—{B' -BV-PoWO-'B”"') 


= ^P„(0> 


cl 


r ryn-r 


B^B 




(B" -B^) 





cl 


d 


= AP,iO)—B''(B"-B^) 
c\ 

Similarly 

P{Si) = B’'{B\-B'')—APM r + l<«</?-l 
" c\ 

where = (4.2.16) 

B= + ^ — - 

c{n-e) c{iJ,-e) 


From equation (4.2.7) 

K+\ (1) = {\ + cfi-{c + l)e)P^+i (1) + c{/J. - e)p^+2 0) 

Taking Laplace Transform - 

5P.+1 (1) - Pr^i (0) = +cp-{c + l)e)P r+i (1) + c{p - e)Pr +2 (1) 

c{p — ^)p r +2 0) ” +cp — {c + V)e + 5)Pr+l (1) 


r 


P.+2(l) = 


Ij -e + c{p-e) + s 




\Pr.l{l) 


V c{p-e) J 

Pr+l(X) = (1 + + D')P r+l (1) 

where D = ■ D’ = —— - 

c{p-e) c{p-e) 

From equation (4.2.8) 

p „ (1) = -(ylj + c// - (c + l)e)P„(l) + c{p - e)P„+i (1) + (>^ - ^)Pn-\ (1) 

r+2<n<R—\ 

Taking Laplace Transform- 

S'P«(1)-P„(0) = -iXi+cp~ic+l}e)P„(l)+cip-e)Pn^il) + 

(;l, -e)P«-i(l) 

SPn (1) - SPn-i (1) = -(A, + cp-ic + l)e)Pn{\)+ c{p - e)P «+i (1) 4- 

(A^-e)P„-.(l) 


146 


p„+l(l)= 


A, - e + c(// - e) + s ^ 




A^—e+s 

c(^-e)^ 


P-(l) 


c(//-e) 

F„.i(l) = (l + D + D’)F„(l)-(D + D')F„-i(l) 

Fut n = r + 2 

(1) = (1 + £> + D')Pr.2 (!)-(£) + D')Pr.l (1) = (1 + D + D')(l + D+D’) Pr^^ (1) - (D + D')Pr^l(l) 

Pr^l) = [{l + D + Df -(D + D')]P r+1 (1) 

Pr+3 (1) = [l + i)' + D'^ +2D + 2DD' + 20 -D- £>']p.+i (1) 

P.^.3 (1) = [l + + D'^ + 2DD' + D+ D']Pr+i (1) 

= [{l + D + D') + (D+Dfp r+1 (1) 

p„3(l) = lz(P±PZp„,(,) 

l-iD + D') 

Putn = r+3 

p^, (1) = (1 + D + D')Pr ^3 (1) - (D + i)')Pr+2 (1) 

= (l + D + D')[(l + D+ Dy - (D + D')]Pr+i (1) - (D + 0)(l + D+ D')Pr^i (1) 

= [(1 + D + D’f -{D + D')(l + D + D')-{D + D’){1 + D + D')]Pr+i (1) 

= [(1 + Z) + D'f - 2(1 + D + D')(D + D')IPr+i (1) 

= (1 + D + D')[(l + D + D'f - 2{D + DOjPr+i (1) 

[1 -{D + D')](l + D + D'){{\ + D + Dy -2{D + 0)] 


[l-(D + 0)] 


Pr+l(l) 


^ [1 - (D + D')^][\ + (D + D ')^] 
[1-(D + D')] 

Pr^4(l) = ~ P^(\) 

\-{D+D') 

Similarly 

P„(l) = Pr+(«-.)(l) = Pr+l(l)| 

-e 


f\ 2 


Where D ■■ 


■a 


\-{D + D’)”-’ 
\-{D + D') 
s 


c{fi-e) c{fi-e) 


P r + l (1) 


r + l<«<P (4.2.17) 
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= (l + iC + J^')Pr+l(l) 


\-K 


R-r 


\-K 






V 


1-^ 


- (1 + (1 - B)— ^Po(O) 

c! 


:(i+is:+is:')p.+,(i) 


\-K 


R-r 


l-K 


■K{l + K)Pr^x(\) 




V 


\-K 


\-K 


[(1 + ii: + is: ' )(i - /i: - is:(i + ii:)(i - ) - (i - i^^ )] 


i-is: 


^ ^fr_^R-r.X ^^2 _ - K + -1 + i^"] 

1-is: ^ 


P«.2(l) = Pr+l(l) 
P/J^2(l) = P.+l(l)i 


\-K 

__ 


Put n=R+2 

p«.3a)=a+^«^2(i)-j^/w(i) 

~X-KP^\r. 


=a+iq 
=a+. 

=(l+Kj^ 

=a+4 


i-is: 

l-K 


Pr^(l)-^(l+K)PR(l)-KPm(l)-^iu(oi 


PniQ-iqa+^'^a: 


-^2_^+2--L 


l-K 

k^-k!^^ 


l-K 


M-K(^+K)Pr^m 


Prrfa)-ii3:i+^'Hi(i: 



-EPmiG) 


|+J20r(5^ -Pp^A^Of 
d 


\+KB^(l-Bf^AP, (G) 


=(i+ii:)| 


2 xrp+2~r 


K^-K 


l-K 


P„,(l)-ii:(l + K)P,^i(l) 


V-K 


R-r^ 


1-J5: 


+ K'Pr+,(l) 


l-K 


R-l-r 


l-K 


(1 + K)(K" - - K(1 + K)(l - K^-' ) + K^ 


t'l — jl^R-t-r ^ 

l-K 


l-K 
+ K(1-K) 


+ KP,«(1) 


KK'j 

l-Kj 


Z:lO) [j^:^ - + K^ - + K«-^"’ ^K^ + K*-'-"^ + K^ - K^'^*’ + K - K^ - k] 

1 — iC 


■■Pr.x(li 




l-K 
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P«.3(l) = Pr..(l) 

P«.3(1) = P...(1)| 
Similarly 


\-K 

^R+3-R _ ^R+3-r 


p„(l) = P«+('>-«)0) = P '•+i(i) 


jrn-R_^«-r 


\-K 


:Pr.l(l) 


K" 


P„(1) = P,.,(1) 


\-K ] 

IN/I 


(P)+D‘y 


(D + /)'r''-(D + D')- 


l-(I> + £)’) 


n>R 


(4.2.18) 


Performance measures of the Model 


Using i/ze system size distribution we can analyze the system behaviour 
of the mode/ by deriving the system characteristics under equilibrium. 
l)The probability that the system is in state (/ = 0,1) is 


^(0 = Sp «( 0.0 = W 


«=o 


Therefore 


C-1 


w=0 


/>(0) = £i>.(0) + X^"(0)+ 

■wssr+l 


M=0 


TB” = + 5 ^"’ + — _ + 5 '- =f--^ 

^ l-B 






Wasr+l 


^(Mlpo(0) + — Po(0)| 
p(0)=2^ »! " - a ’ 


C D'*+l 


1-p 


«as0 


+ C 

c! 


APo{0)B^ 


j^r+l 

i-p 


P^(P-r) 
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w =0 


nl 




Fo(0) 


-B 


r4-l 


+ 


d 


[ b ’' -B^ YPo {0)B'^' 


r nR 


B’'-B 

l-B 


c ! 


[b'' -B^YPo{0)B^^’'iR-r) 


^(^/^„(0) + f 4 Vo (0) 

L "■ 

n=0 




B' 


l-B 


U'-y 


Po(0) 


B 


r+\ 


l-B 


>r4-l 


^MO)^ 

d l-B d 


^ (B'' -B^y^PoiO)B^*’’{R-r) 


F (0) = I 




n=0 


n\ 


C! 


l-B 


A{R-t)B 


R+r 


Po (0) 


Also 


P(l) = ^P n(l) + ^P «(1) 

n=r+l rp=R+\ 

'l-iD+aT" 




/7=r+l 
00 


+ 


n=R+l 

Let D+D' = <D 

R 


L 1-( D + D ') 

( D + i 7)-^-( D+zyr 

l - p + zy ) 


{D+lTf 


?(!)= £ P „,0) 

n=rr+l L ^ ^ 


W=/?4-l 


1 -cD 


?r.l(l) 


1-^ i 1- 


IJCSp+l 


R _ i^-r \ 

~ d > 


/j 


— — [(1 - 0) + (1 - ) + — - + (1 - ) 


y ( i >-^ - < D "'^ 


1 -<D 




■ — — oo| 


P-i(l) 


(l-d> 


[(1+1 + + i )-(< I >+€>^+ — -+ d >^"^) 
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(4.2.19) 




■<D 


-r\ 


1-0 


= Pr+l(l) 


<I)^+' (1 + (J) + + oo) ] 


( 1 - 0 ) 


(R-r)- 


1-0 


'^O'^ -O 


-r\ 


1-0 


o 


«+l 


1-0 


= Pr^l(l) 


(ig-r)(l-0)-0 + 0^~''^'+0-0 


/?-r+l 


P(l) = P..i(l) 


R-r 

1-0 


(i-^r 


= Pr+l(l) 


R-r 


1-{D + D') 


(4.2.20) 


{ii)The probabiliy [Po(0)] that the systemis empty can be calculated from 
the normalizing condition- 

fp„ (0) + =P(0) + P(1) = 1 

n=0 ri=r+l 

Substituting forPifr) and'Pil) in the above equation and simplifying 

=l:^^p.(o)+- 

^ ^ V ^ 

+ P...(l)((P-r) /(I -(£) + £>’))) 




■A(R-r)B’’" iPo(O) 


l-B 


n=0 

f 


c\ 




l-B 


A{R-r)B 


R+r 


\ 


\Po{0) + 


B^^’'(1-B) 




APo (0) {{R - r) /(I - iD + D'))} 


gSp „( 0 )+‘’ 

«=0 


cl 






l-B 


A(R - r)B 




iPo(0) + 


r 


Xq e 




Brjl - B ) R-r 


\ciju-e) c{p-e)) 


d (r-B^) (l-(D + D')) 


7;r^»(0) 


I 



n /oC 


to «! 




B''{\-B) R-r 
c! {B- -B’')\-{D^D') 


{BCy C 1 B' 


■ + 


n=0 


1-5 


I -(£> + £)') 


to «'• 


1-(D+D') 
(B-{D + D’)\ 


^{Bcy cM B^ 


( 4 . 2 . 21 ) 


c! [1-B 

(///)TheprobabiU|f that the numberof units in the systemlies betweenr and R is obtainedas 

P(r < « < 5 ) = ( 0 ) + ( 1 ) 

"l-CIH-D')""! 


w— »• 


wsr+l 


l-(i)+D') J 

= 5 ''— ^ 5 o( 0 )[( 5 '' - 5 ") +— +( 5 "'"' 

c! * 


-(D+D') 


[{(1 - (Z ) + £>')} + {(1 - (D + Z)')^ )+ {l - (Z) + Z/)^ )+ + {l - (Z) + zy )] 

= B''—APo(0)[(B^+B^"' + + 5"’')-(5'* +5^ + +5^^)] 

d 


Pr.^(}) 

l-iD + D') 


B^—APoiO) 

cl 


R-r-lf^iD + DT' 

{R-l-r)B^ 


1-5 




1-(D + Z>')L 


R-r -I- '^{D + D'y-'' 


«3=r+l 


= B^—APo(Q'i\ 
c\ 


' B^ -B 


1-5 


-(5-l-r)5'' 


5r.l(l) 


1 -(Z) + Z>')L 


/?-r 


5-r-l-5](Z)+ Dy 




c! 










c! J[ 1-5 




!-(£» + £>■) 






1-5 
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\-(D+D') l-(£) + £)')t? 


/‘oCOF 


.?.(0)(=%V 




(R-1- rXl -(D+D')-(R-r- 1)(1 - B) 


(l-B)^(D + D'r 

n^l 


(1-5) 


-^5« (5-1-r^ 


l-(£)+D’) 


(1 -(£) + £>■) -(1-5) 


1-(D + D ) 


(l-(i) + Z)')) 


-IzlL^y^D+D'r 

l-(D + D)ti 


1— 0Tg )_-g(D ^ z)')” 

(1-5) 1 1-(D + D') 1-(D+Z))S 


(4.2.22)) 


(D + D) 1-(D+Z))S 
(iv) The conditional probability that the system is in state - 0 when the system size lies 
between g and R is 


R-l 




Pio/r<n<R) = ^^ 


P(r ^n< R) 

-B‘^)^APoiO) 


Pom 




B' 


L - . g _^|;(Z)H-i))4l 

[(1-5) 1 1-(D + Z>) l-iD + D)^, jj 




cl 




Po(0) 


fr^') 

■ 

c 

B' 

U'-j 



(1-R) 


[ l-iD+D) 1-(D + D')^ ■ 


(5'-_5«)-i 

^(.B~-B‘) 


Xl-B) 

\<,R-r)(B-(D + D) l-B Vm , nr 

l-(D + D) 1-(D + D)tr 


4- + B^-^)-B^(l + l + + 1)] 

1/ ■ „J(P-r)(P-(D + D')) 1-B DV 

-B'^iR-r) 




r Di? 


B''-B 

1-B 


X. 


1-5) 




!-(£>+£>') 


l-iD-hDytt 
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(v) The conditional probability that the system is in state - 1 and ^ven that the 

system size is between r and R is obtained as. 


«=r+l 


Mean queue length > 

We calculate the expected number of customers in the system let Ls denote the over 
average number of customers in tire system then we have 
Ls = Ls +Ls, (4.2.25) 


where 




n=0 

C-1 




«=0 


d 


w=r+l 


R-\ 




/i=r+l 


/j=r+l 


Po(0) 


/ 2 =C «=0 

S'=c5‘^+(c+1)B‘='"'+(c+2)S‘^*2+ 00 

55 = cB‘=^’+(c + 1)B‘^^'+ — 

(1 - B)s = cB‘^ + B^^^ + B^*^ + 00 

= cB^+B‘^*\l + B + B^+ — oo) 


B 


c+l 


s = 


s = 


(l-5)s = cB^-f 

1-B 

cB^ B^*' 

1-B^ (l-B)^ 

(l-B)cB‘= + B^^* 

(1-5)^ 

B^ (c-(c-l)B) 

(1-5)^ 

R-l R-r-2 

(11) AB’ + + 

n=0 


n+r+1 


S = (r + 1)B"' + (r + 2)B'*^ + oo 

Bs= (/- + l)g''^^ + 00 

(1 - B)S = (r + 1)5''^’ +B'*^+ 00 

= rB^*' + B'*' + B'*^ + oo 

= rB'*'+B'*'(l + B + B^ + oo) 

^ 1 ^ 


= rB'*' + B'*' 




1-B 
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^ 1- 




,r+l 


l-B (l-Bf H-B) 


(1 - BY 

B^{R-rB) 

a-Bf 


B'*\(T + l)-rB) 

a -BY 


where r + 1 = R 


fi-i 


AB^ Y, nS" = AB' 


/ 7 :=r+l 


B\R-rB) 

(I -by 


= AB 


H.r (R-rB) 

(I -by 


^+1 r 1 

{IU)AB^*' J « = AB^*' [(r + 1) + (r + 2) + — +{r+(R-\- r)J 

= AB’^^'' [r(l + 1 + 1 — ^ - +1) + (1 + 2 + 3 + — +(i? - 1 - r))] 


n-r+\ 


= AB'^^'' 
= AB'^^^ 
= AB'^^^ 
= AB'^^^ 


r(R-r-l) + 
(R-r- 1 ) 
(R-r-l)l 


R-l-r 


r + - 


2 

R-r 


[2 + (i?-l-r-l) 


2 J 

2 r + R-r ^ 

2 J 

^(R-r-l)(R + r)^ 


V 




= AB''^''{R-¥r){R-r-\)l 2 
Lso — 


\)B ^ AB'"\R -r B )_ _ ^ _ y 


(l-B)" (l-B)" 

■Po(0)‘^% + ^M(a!o"/"*)P«(®) 


( 4 . 2 . 27 ) 


«=0 


Ls, = £ <(l) + 2<0) 

wssr+1 n^R 

' x-jP^D'T' 

I \-(D+D') J 
1 


n\ 


= E 

n^r-^ 

= ?r«(l) 


Pr..(l) + X«-(^ + ^T 


n»R 

^ n- YniD+DT''A(I^ + Dr" -iD + DT^n.iD+D^^ 

■ wsix-4.1 W**/? 


(£> + £)0~^ -(£> + !)')' 
4-(Z) + £>') 


Pr+l(l) 


(i-(D+i>'))j_„t;i .til 

Simplifying as in previous case 


= Pr^(l) 


'[i-(d+i 7)1 


-(^+r)(R-) — +^llM-£l)"*-(B+l7)'^]^rt(D+£l)" 

2 »r+l /eR 


(/)|}<£>+Z?)'^ =(r+l)(iD+i7)+(r+2)(D+J[?)" +(r+^(Z)+Z?)' + Mr+R--l-r)(p+i::^f-^ 


157 



s = (r + \)iD+D') + (r + 2)iD+D'f+ +(r+(R-l-r)(D + D')^-^-' 

(D+D')s=(r + l) (D^D'f+ +(r + .R-2-r)(D + Z)')'^'^'' +(r + R-r-l)(I>+Z)')‘^'^ 

Cl-(D + D')s = ir + l)(D + D') + (3+ D'f + (D + JD')"'"' + (^ + D'f" -RiD+ DY"' 

= riD+D’) + (D+D') + (D+D'f+ — (D + - R (D + 


, (D+£>')-(D + £>'f -^ 
1-(D + D') 


= r(D + £>')+— — ^ - Yy. ’ R(D + D')^" 


r(D + D') (D + D')-(D + D')"’''^' R(D + D') 

S ' "4" 


R-r 


5 = 


1-(D + D') [1-(D + D')]' 1-(D + D') 

^ r (D + D') ( ^ - D + D') + (D + ) - (D + - R (D + [(1 - (D + D^)] 

[1-(D + D')f 

r (D + D') - r (D + D')^ + (D + D') - (D + -R(P + D')'*"'' +R(D + 


[1-(D+D')f 


(//) f;n(D + D')’ 


n=R 


s = R (D + DY+ iR + 1) (D + + (R + 2) (D + + 

(D + DOs = R(D + D')''^*+(R + l)(D + D')''^^+-- 
(l-D + D')s = R(D + DO*+ (D + D')^'"' +-^ — oo 


1-(D + D') 
^ R(D + D0^ ^ (D + DO 


R+l 


[1-(D + D')] [1-(D + D')] 

R(D + D')'' [1-(D + D')] +(D + D') 


R+l 


[1-(D + D')r 


R(D + D')" -R(D+JD^’ +(D + D') 

[1-(D + D')r 

_^ri( 

FpT5o]L2 

r{D + D’)^-iD+D') + (D+D'f-'*' +R(D+D'y'' - R(D + D'f-"*' +j(Z)+D')''' -(£)+/)')-'•]* 


kR +1 






' R(D + D'Y - R(D + D'f*' +(D+ 

1 [l-(D + D')f j 

^{R+r)(R-r-l)+\-riD+D’)+r(D+D'Y-(D+D')+(D+ + ffY'' 
- R(D + D'Y-'*' +R-R(D+ D'Y" - R(D + £>0 + R(D + + (£) + Z)‘ ) - (D + D'Y"*' }Il - (/> + £>')? ] 
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1 -(£) + £)') 



^{R + r)(,R-r-\) + ^-R{D + D')-r{D + D') + ri,D^D'f }/[l - (£> + } 

P J ^ ihrrnVr + 

1 -(£>+£>')] 2 ^ I [\-iD + D')f J 


= P,,, (1) ! \LiR + rXR-r-\)HR-r(D + D')) /(I - (Z) + D')) 

1-(D + D')JL2 


~ ^so -^si 

= - (i? - rB) (aB'^‘^I(\ - B)) 

--AB'^*''(R + r)(R-r~\) ]Po(0){c7cl)+i]«(ao7«')^o(0) 

2 n=0 

+ ^ (1) [-(i? + r)iR - r - 1) + (i? - r(D + D')) /(I - (D + D')) 

(l-iD + D^)) \2 


(4.2.28) 


(4.2.29) 


(yil) Ejqjected waiting time: using the Little formula the e5q)ected 
waiting time of the customer in the system is calculated as 


W, = LJX 


(4.2.30) 


Where A is the actual mean arrival of the system and it is given by 
A = 2oP(0) + iiP(l) 

Substituting the values of P (0) & P (X)Jrom equation (4.2. 1 8) & (4.2. 1 9) 


«=o 




tm 


Particular Cases 

(1) Single Server Poission queueingmodel with controllable arrival rates 

When C = 1 the M ! M ! C'xvX erdependent queueing mod el with conbollable 

arrival rates reduce to M ! M H int erdepend queueing mod el with controllable arrival rates. 

Thus when C = \We have the following 

For \<n<c-\,Pn (0) is not possible. 

For \<n<r 
P„(0) = ao'’Po(0) 

For r <n < R — \ 

Pn (0) = [ao” - (0)[5'mce when C = 1, 5 = ) 

WhereA = [oTo' - ]* 

Equationin(4.2.\7) gives 

l-(a, +a2)"' 


SimilarlyPn (1) = Pr+i (1) 


P„(l) = Pr+l(l) 


(«! +«2) 

Where Pr+i (1) = (1 - ao)APo (0) 

a, 


l-(a, +a 2 ) 

J(a, + 02 )”'' - (a, +aiy' 


1 - (a, + or2 ) 


r + 1 < n < R (fFhereC = 1 D = or,, D* = eta) 
n>R + l 


P(0) = 


l-ar 


- A(Ji - r)ao 


R+r 


Po (0) From equation 


P(l) = [(^f-r)/(l-(a,+a,))JP r +1 (1) From equation 

From Equation (4.2.29) gives 

L{S) = -{R- ra^ )A ao'"' /(I - )) 

JL0 


(ii)M\M\C Model 

When e = 0 andX^ = Aj = 2 transient behaviour of the M ! M ! C 
interdependent and queueing model with controllable arrival rates 
reduces to the usual M f M ! C Model [ly 
Hence, When e~o, s = 0, Aq = A, = A(r = i?) 
equation{A2.A€) gives 


(4.2.18) 

(4.2.19) 


— ... 

fAY 

Yn 

p„(0) = 




/>o(0) 


1 < n <e-l 


For C^n<r 

P„(0) = (A//t)"(l/c"'"c!)Po(0) From (4.2.16) 
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P„ (0)(r <n< R-l) and Pn (l)(r + l<n<R) are not possible 
From equation (4.2. 17) gives 


Where 


P-i(l) 


n > i? + l 


Equation (4.2.20) gives 


AiR-r)B’'*''{B-{D + D') )/(l-(D + £>') ) 


Equation (4.2.29) gives 
, \b‘(c-{c-\)B) 


AB'^^''{R + r){R-r 


'^nictQ lni)PQ{()) + Pr+\{\)l{l-{.D + D''f}* 

W= 0 ' 

[ 2(5 + rXi! - r - 1 ) + (Ji - r(i> + D’ ))/(!- (D + Z) 


tS «! (1-P) 

Where p = Xjjjc 

Thus the resuslts of usual MIMIC Model are obtained as a special case. 


Cficvpter - 5 


The Steady State Analysis of M/G^/1 Retrial Queue 
With Multiple Vacations 


INTRODUCTION 

Retrial queueing system are characterized in a such a way tliat 
arriving customers who find the server busy join the retrial group and try 
again after random intervals for tlieir request. In tliis chapter we consider 
customer’s arrival one by one at a moment but service to all customer’s 
is provided in bulk. A large number of applications related to 
transportation process such as by road transport, rail transport, airplane 
etc. involves bulk type of queue process. 

Research on retrial queues with non-exponential retrial times is 
very limited. Tlie first ever work reported in 1977, by Kapyrin [93] the 
author analyzed and derived an exact analytic solution for the M/G/1 
retrial queue with general retrial time distribution. However, this 
methods was found to be incorrect and the results are erroneous as stated 
by Falin in 1986. Anotlier account of tlie subject is die treatment of a 
G/M/S/0 queueing loss system with retrial (Pour Babai, 1987), where 
returning customers are assumed to be a new form, dependent input 
process, input process, whose parameter are approximated fi'om those of 
the overflow process from the service station. Retrial queues are widely 
used as matliematical models of several computer system, packet 
switcliing, network, sheared bus local area, network operating under the 
carrier sense and it is assumed that at every arrival epoch, a batch of 

■ i . 

primary customers arrives with a given probability. 


If the server is busy or down due to on vacation at the arrival 
epoch, then all these customers join the group whereas if the server is 
free then one of the arriving customers begins his service and the others 
leave the service area and enter a group of blocked customer called orbit 
according with FCFS discipline. 

Single server queues with vacation have been studied extensively 
in the past. These models arise naturally in telephone switching system, 
computer communication system etc. most of the analysis for retrial 
queues concern the exhaustive service schedule Artalezo [6] and gets the 
service policy Langaris [103]. Chaudhary and Templeton have provided 
a comprehensive review on bulk queues and their application Langaris 
and Moutzoukis [101] have considered a retrial queue with structured 
batch arrivals, preemptive resume priorities for a single vacation model. 
Yang and Templeton [148], Falin [53], Falin and Templeton [54] and 
Artalezo have presented a comprehensive survey and related work on 
retrial queue, queueing system with batch arrival are common in a 
number of real situations. In computer communication system, messages 
that are to be transmitted might have random number of packets. 
Comparable work on optimal control policies for batch arrival case is 
seldom found in the literature. This motivates us to develop a realistic 
model for queuing system with batch arrivals. Chaudhary and 
Templeton [33] have provided a coni|)rehensive review on bulk queues 
and their application. Langaris and Moutzoukis [101] have considered a 
retrial queue with stmctured batch arrivals. Preemptive resume priorities 
for a single vacation model. 


It is assume that a retrial tin^ starts when the server completes a 
services and becomes free so that customer waits for 
completion of repair period. If the server is idle, an arriving customer 
(new or retrial) must starts the server, which takes negligible time. If the 
server starts successfiilly then the customer get services immediately, 
otherwise the repair to be done then the server commences service 
immediately and the customer must leave for the group and make a 
retrial at a later time. 

The service of customers are independent random variable 
with common distribution B (x), density function b (x), Laplace - 
Stieltjes transform B*(q) and first two moments Bi and B2. Similarly the 
successive repair times are independent and identically distributed with 
probability distribution function H (x), density function h(x), Leplace - 
Stieltjes transform H* (q) and first two moments D, and D2. 

Further, the vacation time has distribution function v (x), density 
function V (x) Laplace - Stieltjes transform V* (q) and first two 
moment Vi and V2. 

The state of the system at time t can be described by the Markov 
process {N (t); t^ 0} = {(C (t), X (t), Xo (t), Xj (t), X2 (t), X3 (t)); t^O} 
where C (t) denotes the server state 0,1,2, or 3 according as the server 
being fi^e, busy, down or on vacation respectively) and X (t) 
corresponds to the number of customers in the orbit at time t. If C (t) = 0 
£md X (t) > 0 then Xo (t) represents the elapsed retrial time. If C (t) = 1, 
then xi (t) corresponds to the elapsed service time of the customer being 
served. If C (t) = 2 and X (t) > 0, then X2 (t) represents the elapsed 
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repair time, if C (t) = 3 and X (t) > 0 then X 3 (t) represent the elapsed 
vacation time at time t. The function r (x), m (x), h (x) and n (x) are tiie 
conditional completion rates (at time x) for repeated attempts, for 
service, for repair and for vacation, respectively. 

i.e. 

r (x) = a (x) [ 1 -A (x)]'* m (x) = b (x) [ 1 -B (x)]'‘ 


steady -state distribution : 

we study the steady state distribution for the system under consideration. For the process =0}, 

we define the probability. 

P„ix,t)clx^ P{c(t) = 0,x(i) = ri,x = ^,(t)<x+dx]foT t = 0,x = 0&n = 'i 
Q^ix,t)dx= P{c{t) = l,x(0 = n,x = ^,(/) < X + £&}for / = 0,x = 0 & n = 1 
R„(x,t)dx = P{c(t) = 2,x(t) = n,x = ^^{t) < x + ^}for / = 0,x = 0& « = 0 
V„{x,t)dx=P{c{t) = 3,x(0 = n,x = ^,{t) < x+rfr} for t = 0,x = 0& « =0 

By supplementary variable technique we obtain the following system of equations that govern the 
dynamics of the system behavior. 


dP„{x,t)^ SPniX’j )^ -(A + r(x))P„(x,0,n= 1,2,3..- 


dt 


dx 


dt dx 

-ix + juix))Q„(x,t)+^^^^ CiQ„^,(x,t), n= 1,2,3... 

dt dx 

aR,(x,f)^ dRi(x,t) _^ _(^;i + jj(x))R,(x,t) + riixyZ^kRXx,t) 
dt dx 

£M£!£l+ -{/l + 7(x))iZ,(x,0+ Jjix) C*if^*(x,/),n = 23,4... (5.2.5) 


( 5 . 2 . 1 ) 

(5222) 

(5.2.3) 

(5.2.4) 


-(X + V(x))V,ixj)+V{x)^CtVXx,t) 
dt dx M 


.(A+v(*)F.fc')nx)c.»'„,(j,».n= 

dt dx 


The boundary conditions are :- 

« ® “f 

Pn(0,t)= p Jfi„(x,/)/r(x)<ir+ ^R„ix,t)r}{x)dx+ jF,(x,0F(x)<ir, n-12»3-.. 


(5.2.6) 


(5.2.7) 


(5.2.8) 
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Qo ( 0 ,t) = a IPi{x,t)rix)dx 

0 

00 


CO 


Qn (0,0 = a (x,t)rix)dx + oA Jc* {xj)dx, n = 1 ,2,3 
0 0 


( 5 . 2 . 9 ) 


(5.2.10) 


Ri (0,0= a ^P^{x,t)r{x)dx 

0 

CO 00 

Rn(0,0= a \P„{xJ)r{x)dx+ aA^C^P„^^^^{x,t)dx, n=23,4... 

0 0 

00 

Fo(0,/)= j Q^{xj)iLi{x)dx+ 1° V^{x,t)v{x)dx 

0 


(5.2.11) 


(5.2.12) 


(5.2.13) 


V„{0,t)= \Q„{x,t)fi{x)ctx n= 1,2,3... (5.2.14) 

0 

We assume that condition a2c, {B^ ^-qV^) + a^\ + c^AD^)+ c,(l- j 4*(2.)) < 1 is fulfilled so that we can 
set the limiting densities. 


P„{x)= Lim,^„ P„(a:,0 forx = 0 & n = 1 
Q„ (x) = Q„ (x,/) for X = 0 & n = 0 

R„ (x) = R„ (x, /) for x = 0 & n = 1 

V„ (x) = Lim,^„ F„ (x, 0 for X = 0 & n = 0 

letting /-w! in equations (5.2.1 ) - (5,2. 1 4) we get 
dP (x) 

—5^= -(A + r(x))P„(x) n= 1,2,3.... 
dx 

~ ^^^ = - ( A + /r(x))go (x) + 

ax i-i 


(5.2,15) 


(5.2.16) 
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dQM) 


dR,{x) 


■--{A. + /i( x))2„ (x) + Mix)Cf; {x) n = 1 ^,3 - • . 


■ (2 + r]i.x))R^ (x) + T]{x)^ i?, (j:) 


+ Tj(x))R„ix) + 7jix)Cf:R„^g(x) n = 2,3 4... 
ax 


(5.2.17) 


(5.2.18) 


(5J2.19) 


dV,{x), 


= - (A + v(x))Fo (x) + v(x) Cjf F; (x) 


(5.220) 


dV„{x).. 


■(;i+v(x))F„(x)+v(x)CjfF„,.jc(x) n- 1,2,3. 


The steady state boundary conditions are : 


(5.221) 


OO «3 w 

Pn(0)= F jg„(x)/^(x)<3!x:+ Ji?„(x)77(x)c6c+ Jf„(x)F(x)c 6[:, n=l,2,3-.- (5.222) 


Qo (0) = a JF,(x)^(x)^^r 


(5.223) 


00 w 

Qn (0) = a Jf„^i {x)r{x)dx + oA^Ci, ix)dx , n = 1 ,2,3 . 


(5.2.24) 


Ri (0)= a J/J (x)r(x)<5&(: 


(5.2.25) 


00 ov 

Rn (0) = a jF„(x)r(x)a!xr+ aX {x)dx, n=2,3,4 


(52.26) 


00 

Vo(0)= ]Q^{x)^ix)dx+ |Fo(x)F(x)c& 


(5227) 


q\Q„{x)^^{x)dx, n= 1,2,3. 


(5228) 


And the normalizing conditions is 
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(5.2.29) 


^ ]p„(x)£fe+ \Qn{x)dx+'^"\R„{x)dx+'^"\v„{x)dx: 

n~l 0 ^*=0 0 «=1 0 W3?0 0 

We define the probability generating fiinctions 


P(X,Z)= ^P„(X)7^ 


n=\ 


Q(x,z)= £awz" 


«t=0 


R(x,z) =f^R„ (x) z" V(x,z) = X V„{x) z" 


/J=sO 


The following theorem discusses the steady state distribution of the system. 
Theorem 1 If the system is under steady state that is 

oOc, (5, + 9 1^ ) + a (1 + c,/l£), ) + c, (1 - / (2)) < 1 
then the steady-state distribution ofN(t) = = 0}is obtained as 


P(x,Z) 


^ Fq (OXP + qV' (2))[l - V' (2(1 - c(Z)Z-^^ ))]e-^[l - ^(x)] 


(5.2.30) 


D{Z) 


Q(x,Z) = 


R(x,Z) 


_ V, (0) tz(P + qv' (2))(l - (2(1 - c{Z)Z ))\c{Z)Z A' (2X1 - c(Z)Z )jg-^ [1 - Bjx)]} 

DiZ) 

(5.2.31) 

Fq (0)^(P + qV\Z)il - F-(2(l - c(Z)Z-^^)))|c(Z)Z-^^ A'(A)(l - c(Z)Z-^'*^ )>-^[l - fl-(;c)]} 

D(Z) 

(5.2.32) 




K mp + qy'mUz)Z'"’' + (-tXl -c(Z)Z-^*^ yjoB' (M\-c(Z)Z-^'^ )]+ oZtf • (A(1 -c( 2)2-^'^ )) - zje-^[l - F(;c)] 

EKZ) 


(5.2.33) 

Where D(Z) = 

^Z)Z-“ +/(2Xl-c(2)Z'«)}|^*(2(l-c(Z)2-^))jCF+^F(^l-c(Z)Z«)^+aZW'(2(I-c(Z)Z-^)) }-Z 
and the probability Vo(0) is to be determined fiom the normalizing condition. 
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Proof : - 


Multiplying both sides of equation (5.2.1 5) by z" and summing with respect to n -1 ,2,3 - ■ • • we obtain 


= -f^iX + r(x))P„(x)^ 

dX 


^-iA + r{x))P{x,z) 
dx 


(5.2.34) 

Now multiplying equations (5.2. 1 6) & (5.2. 1 7) by z" and summing over n= 0,1,2. oo yields. 

(5.2.16) 


-(A + Mix) Qo(x)+ Mix)ZC,Q„{x) 
dx 


i=l 


— ^ ^ — (2 , + //(jt)) Qn(x) + )U(-x) Ck QrH-k(x), n 1 ,2,3 . 
dx 


(5.2.17) 


+ y ^gg Mz" = -( 2 .+//(x)) 

t^^ dx ^ dx 


+ 


//(^) 


+'^QAx)z" 

/jssO ^*=0 

i]Qa(x)z'' +c*e,.*(x) 2 '’ 


= - (2 + ;U(X)) Q(x, z) + 


dx 


MMf,{CAQiix)+Q2(x)+----+Q,-i(x)+QAx)+Q.My 

^ffssO 

= . (X + /j(x)) Q(x,z) + 
dx 

//(x)[c* Ij2l (x) + fi2(x) + + !2*-l(x) + l0i(x) + 0*+i(x) + .....'=O }z j 

= - (1 + Mx)) Q(x,z)+ Mx)tc/|;Q0*^(x)Z'’ 

nwO \ 


dx 


- (2- + fj{x)) Q(x, z) + //(x) ^ g*+, (x)Z""*^ X! ^ 

ffjBsO 
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= - a+M^)) Q(x,z)+ /z(x)C(Z) 

n=0 

= - (A + /d{x)) Q(x, z) + ^{x)C{Z) Z'^*g(x, z) 

M5zI 1= . (Z + Mx)) Q(X,Z) +Mx)C(Z) Z-^'Q(x,z) 
8x 


=-[z + //(x){l - C(Z)Z"'* }1q(x,z) 

dx 


(5.2.35) 

Now multiplying equations (5.2.1 8) & (5.2.19) by Z" and summing over n = 0,1 ,2 yields 

(52.18) 


-(Z + 77(x))i?,(x) + Tj{x)f^C,R,(x) 
dx 


/a I 


- (Z + rjix)) Rn(x) + rj{x)C, (x) , n = 2 , 3 , 4 . 
dx 


( 5 . 2 . 19 ) 


= -{X + rj{x)) 

nssO '>=0 


|;/?,(x)Z'’+Xi?„(x)Z'’ 


WasO 


ff«0 


+ 


= - (Z + 7(x)) R(x, z) + 7(x) 
dx 


J C* («, (x) + ^2 (x) + ... + Rj, (x) (x))Z” 


w=0 


_(z + 7(jc)) R(x,z)+ 7(x)[c*(/?,(x) + /?2(x) + ...+ if*(x)+i?*,,(x) + ....cx>)Z'’] 

& 


• (Z + 7(x)) R(x, z) + ij(x)^Ct 




z" 


jt»i 


= - (A + 7(x)) R(x, z) + 7 (x)£ (x)Z'’ * £ CjZ* 

7a=0 A-l 

= -(A + 7(x)) R(x,z)+ 7(x)C(Z)Z"“/?(x,z) 

, [z + 7(x){l - C(Z)Z-^* })/?(x,z) 
dx 


(5.236) 


171 



Again now multiplying equation (5.2.20) & (5.2.21) by Z" and summing over n - 0,1 ^.....yields 


(A + F(jc)) V(x,z) + V{x) CiZ)Z-^^ R(x,z) 


PnO) = P ix)pix)dx + (x)jj(x)dx + jF„ (x)V (x)dx 


±p.io)z-='L ^|aw/<w*z"+ZJ[^.W’7(^)*F+X 

w=sl w=0 0 _ «.0 « 

00 00 00 00 

P(0,Z) = P \Q{x,z)p{x)dx+ \R{x,z)r]{x)dx+\V{x, 2 )V{x)dx- p\Q^{x)pix)dx- 


I Vq {x)V (x)dx (5 .238) 


From (5.2.23) & (5.2.24) 


Ji^(x)r(x)ife 


|?„^,(x)r(x)dE)c+ al jQP„^t^,(x)i& 


n+k-^l 


Q„ (0)Z" = X (xyixydxz" + a 


Q(0,z)= - \Pix,z)r{x)dx +—Z-^' |C(z)P(jc,z)iic.c(z) 


a 

z 


^P{x,z)r{x')dx + Xc(_z)z~^’‘ JP(j:,z)<ic 
.0 “ - 

From (5.2.25) & (5.2.26) 

00 

Ri(0)= Cf J/!(x)K^)ci)(: 

0 

00 00 

Rn(0)= a JP„(x)r(:>f)i&+ aX^C^P„^i,{x)dx 

0 0 

|;i?„(0)Z" = aY^\p„{x)rr{x)dx+ a^-^'f^C,z"\P„,,i_x)dxz 


/»=1 0 


n-\ 0 


(5.2.39) 


(5.2.25) 


(5.2.26) 


00 

R(0,z)= a \P{x,z)r{x)dx+ aM:{z)Z'^^ \Pix,z)dx (5.2.40) 

0 0 

From equation (5.2.27) & (5.2.28) 

00 00 

Vo(0)= ^Qo{x)^i{x)(h+ ^Vf^{x)V{x)dx 

0 0 

CO 

Vn(0) = q ^Q„ {x)fi{x)dx , n = 1 ^,3 • • • • 

0 

|]F„(0)z” j0„(x)ZXx)c&+ jFo(jr)F(x)a!r+ Ja(^)M^)^ 

w=0 /i=I 0 0 0 

^F'„(0)z" =q\Q„ix)f^{x)dx- q^Q^{x)fi{x)dx+”^Vaix)V{x)dx + J’(2oWM^)‘*c 


(5.2.27) 


(5.2.28) 


n=0 


to « "• 

V(OZ:)= q\Q{x,z)nix)dx + \v^{x)Vi,x)dx-^ {\-q)]Qoi.x)M{x)dx 


m 


(5.2.41) 


00 w 

V(O^) = q ^Q{x,z)^i{x)dx+ jVo(x)V (x)dx + P ^Qo{x)p{x)dx 


Solving the partial differential equations (5.2.34)- ( 5.2.37) 

r 

-^x~^riu)dv 


P(x,z)=P(0^)e 


(Xx,z) = Q(0,Z) e 




R(x,z) = R(0,z) e 






V(x,z) = V(0,Z) e" 

From equation (5.2.1 6) -(5.2.20) we get 

j( 

-(A + M(x)a(x)+ M(x)tc,Q,(x) 

dx /=! 


(5.2.42) 

(5.2.43) 

(5.2.44) 

(5.2.45) 

(5.2.16) 


The boundary conditions when x=0 so gives Qo(x) - Qo(0) So the differential equation (5J2.16) on 
solution gives 


Qo(x) = Qo(0) e 

_(^ + F(:c)Fo(x)+ F(x)XQF,(x) 

dx /«i 

-Ax~{\-C\)Vy(u)du 

Vo(x) = Vo(0)e 
From Equation (5.2.46) 

1— C ) f jf4{d)dy /A\ / \ Jf 

^ QMMix)dx 


(5.2.46) 
(5.2.20) 

(5.2.47) 


|a(x)//(x)<i)c= I e 


j0o(x)/i(^)^ = Qo(0) J e ^ Mix)dx 
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Qo(0)= \Qo(x)Mix)dx 

B*{A) ^ 

00 

Where B*W= 

0 

Similarly 

vm’’ y^]KMVU)d^ 




dx 


)J*PdwW« 


dx 


Where K*(/l)= |F(x)e 
0 

Using (5.2.48) and (5.2.49) in (5,2.27) we obtain 

00 00 

Vo(0)= ^Q^{x)fiix)dx+ jVoix)V(x)dx 

0 0 

Vo(0)= Qo(0) fi* (A) + Vo(0) V*(A) 

Qo(0) 5*a)=Vo(0)[l-F*(/l)] 

Now substituting (5.2.42) in (5.2.39) - ( 5.2.40) and on simplification 

oo 

Q(0;Z)= -[jp(x,z)r(x)d)c + AC(z)Z~^^ jPix,z)dx 

^0 0 

X 

= -[JP(0,z)e r(x)dbc + 2c(z)Z-^*f^^^^ » dx] 


= -F(0,Z) 

Z ■ 




00 r{u)4u 

rix^dx + Zc{z')Z~^'‘ je “ iibrl 

0 


X ' * 

CO 00 -fr(ii)dw 

jXe * <i!c= “ dx 


(5.2.48) 


(5.2.49) 


(5.2.27) 


(5.2.50) 


(5.2,39) 


(A) 
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-rr(«)t/K 

e ” I Xe~^dx 


GO 


j -J r{u)du 

— e ° 
dx 

\ J 




dx 


-jr(u)du 

e " e 


-Ax 


QO -jr(u)du 

je~^ e ° r(x)dx 


= 1- je “ r(x)dx 

0 

Putting in (A) 


Q(0,Z)=-P(0,Z) 

Z 


X 

00 -Ax-^r{u)du 

Je " r{x)dx 

+C(z)Z-^* 

( \ \ 
00 r{u)du 

1 - Je ° r(x)dx 

0 


0 

V J 


a 

z 


P(0,Z) Ia*(A) + c(z)Z-^HI-A*(A)] 

= - P(0, Z) [c(z)Z"“ + ^ * (ZXl - c(z)Z-'* )] 

Z 

Q(0;z) = - P(0, Z) lc(z)Z-'* + A * (Z)(l - c(z)Z‘'* )] 
z 

Now substituting (5.2.42) in (5.2,39) ~~ (5,2.40) and on simplification 

00 CO 

R(0;i)=ajP{x,z)r(x)dx+aMz)Z~^'lPix,z)dx 


=a 


« -^-\r(u)du “ 

Jp(0,z> » r(x)£&+aMz)Z-'* JP(0,Z> » dx 


= aP(0,Z) 


« -yLt-|r(B>dB “ -A*-Jr(ii)ii» 

Je “ r(zr)<fc+;}c(z)Z""*Je » dx 


(5.2.51) 
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At-|r(a)t/K 


r(x)dx + c(z)Z 


P(0,Z) [a * (A) + c(z)Z-'‘ (1 - ^ * a))] 
P(0,Z) [c{z)Z-'* + A * (ZXl - c(z)Z-'* )1 


(5.2.52) 


mathematical manipiilatioii, we get 


03 » 

V(OZ;)= q lQ{x,z)fi{x)dx+ jFo(x)P(x)ii)c+ g 


Q(x;z) = Q(0,Z) e 


lQoix)juix)dx 


^d(P)=--— \Vo(x)Vix)dx 

0(0^) = - P(0, Z) 1C(Z)Z-'* +A* (ZXl - C(Z)Z“'‘ )] 

z 

Putting in (5.2.41), we get 

V(0,Z) = q jQ(x, z)pix)dx + J F(,(x)F {x)dx+ p JS© 


p(x)d&c + Fo(0)F * (Z) + PQM^ * (^) 


:SEp(0,Z) [c(Z)Z-^'+A*im-C(Z)Z-^^)] jj 


fi(x)e “ dx 


Fo (OF * a) + p[f<, (0X1 - W)] 

^P(0,Z) [C(Z)Z-^* +Z*(ZX1-C(Z)Z-^*)1 5*(Z(l-C(Z)Z-^*)) + 

z 

VamV*iX)(l-P)+PV^(0) 

^P(0,Z) [c(Z)Z-'* +Z*(ZXl-aZ)2‘'*)l 5*(Z(l-C(Z)Z-'*)) + 

Z 

Fo(0)[F+9F*(^)] 


(5.2.53) 


Similarly substituting (5.2.43)X5.2.49) in (5.2.38) we get 

00 00 00 

?(0^)= p\Qix,z)fii.x)dx+ \Rix,Z)rj(x)dbc+ \V{x,Z)V{,x)dx-P^Qo{x)fzix)dx - 

0 0 0 0 

ljy,(x)V(x)dx (5.2.38) 


Jf 

« -Zi-(I-c(Z)Z'“ )[;'(“)* ®- -Z*-(l-c(Z)Z-")Ji7(‘'W“ 

= pjg(0,Z)e " /i(jc)<ix:+ ji?(0,Z)e ^ ?7(^)iir + 


]v(0,Z)e 


-Zx-(l-c(Z)Z-”)fK(«)</« 


”)|k(«). 


F(x)d&c- F2o(0)5*a)- F„(0)F*(^) 


P(0,Z) = pq(oz:)b* (Z(i - C(Z)Z-^* )) + r(oz:) h* (z(i - c(z)z-^* )) 
+ V(0Z^)V*(1(1-C(Z)Z'“))- 

P(OZ:) = PQ(OZ;)B* (2(1 - C(Z)Z'“ ))+ R(0 j:) H* (2(1 - C(Z)Z'“)) 

+V(OZ:)V* (2(1 - C(Z)Z”^‘ )) -VXO) (P + ^ F * (2)) 


(5Z.54) 
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Again using (5.2.5 1)- (52.53) in (5.2.54) and solving P(0^) after some algebraic manipulation, we 
get. 

Q(0,Z)=-P(0,Z)[c(Z)Z-'*+^*(i)(l-aZ)''')] (5.2.51) 

R(0,Z)= aP(0,Z) [c(Z)Z-^* +^*(^X1-C(Z)Z-^*)1 

V(0,Z) = ^ P(0, Z) [c(Z)Z-'* + (1 - C(Z)Z-^* )A * m] B * (Z(l - C(Z)Z-^* )) + 

z 

Vo(0) (P+^K*(Z)) (5.2.53) 
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[(i,,-z(z)j-i)r)*A-i] i(r)*A^+d) z(o)°A=z 





a[(„.z(z)D - iXr) / + n-Z(Z)j] .[(( ,^-Z(Z) J - 0^^ 





Proof 







Response Time Analysis to the Bulk Arrival Queue System 
in Communication Network 


INTRODUCTION: 

Queueing theory has been recognized as one of the most powerful 
mathematical tools for making qualitative analysis of communication 
network. In the early 1960s it was realized that queueing theory would 
prove to be an effective tools for studying the throughput, response time, 
and other measures of performance of the network. Since that time, the 
literature of queueing theory and network applications have virtually 
exploded with analytical models for different networks. Such as for the 
telephone system the data transmission rates were not sufficient to meet 
the demand, consequently, planning for special data network be^ 
earlier in Europe then it did in United States. Since there was not yet a 
great demand for data services, planning was directed towards to meet 
the future requirement. The planners of some European nation were m 
many ways more ambitious then those announced by AT & T in United 
States, whose digital network were ejqiected to be un-switched. The stop 
and wait, and continuous error detection and re-transmission schemes 
are widely used in such network. The researchers primary focus was the 
development of an appropriate model for analyzing the delays 
encountered in establishing a virtual circuit through a switch based 
LAN. ■ 

The method developed were applied to many related problems i,e. 
multiple link, telephone calls, receiver transmitter interaction between 


switches and shared peripherals in networics. One of such approaches 
developed by Wilhelm, Neil C. [145] to analyze a disk subsystem in 
which all of the disks had the same load. Whitt, W. studied a similar 
method applied to markovian queueing networks as a way of extending 
the approximations developed by Whitt, W. Kelly, F.P. considered an 
approximation technique for a general class of circuit switched tele 
traffic models. 

In 1986, Fredericks, A.A. [60] developed an approximation 
methods for analyzing the performance of a virtual circuit switch based 
LAN. The basic system consists of N input and N output parts, each 
with it band width divided into L equal segments for the virtual circuit 
firom a given input to output parts, it is necessary to obtain a time, slot in 
each part. The author suggested a simple analytical approximation 
methods for the delay encountered in setting up such a circuit Kumar, 
A. Singh, M.P. and Kumar V. [100] also si^ested a n^thods to the 
M/M^ /1 bulk service queue model for the analysis of mean response 
time of the communication network. Graph theory plays an important 
role in modeling and analyzing several problem concerns vdth various 
aspects of networking field. Also graph tteoretic modeling find 
extensive applications in queueing theory for the networking problems. 
In this chapter we consider a bridge network, consists of four systems 
namely Si, S2, S3 and S4, which are inter connected by communication 
links. In this network we shall have to decided the path sets while a path 
in a multi-graph G consists of an altemafing sequence of vertices and 
edges of the form vq, Cj, vi, 62 , V 2 , ..».en.i, v^i, e^, Vn wi^re each 


edge e; contains the vertices V;,] and Vj (which appear on the sides of e,- in 
the sequence). The number of edges (n) is called length of the path. 
When there is no ambiguity. We denote a path by its sequence of 

vertices (vq, Vi Vn). The path is said to be closed if Vq = Vn 

otherwise we say that the path is from Vo to Vn or between Vq and Vn or 
connects Vq to Vn. 

A simple path is a path in which all vertices are distinct (A path in 
which all edges are distinct will be called a trail). A cycle is a closed 
path in which all vertices are distinct except Vo = Vn. A cycle of length k 
is called a k-cycle. Here we assume that a network having “n” 

heterogeneous system, N = {Si, S 2 , Sn}, which are inter 

connected by communication links. One of the “n” system is called a 
source system and one of the remaining “n-l” system is labeled as 
terminal system. At the source system, the jobs are arriving in the form 
of data packets and follows the poison process. The queueing model 
used here in the bulk arrivals (M^/M/1) type. If the source system is 
busy, them these jobs are store in buffer and processed latter. The job 
propagates over the link from source towards terminals system till they 
reach to terminal system. Once a job reaches to terminals system, the 
terminals system sends an acknowledgment of correct processing of the 
job to the source system the aspects of the model is the service 
mechanism. Each arrivit^ data packet of jobs has variable number of 
modules, which arrive in bulk and is then served by the system of the 
network. In this model, we have obtained the mean response time. For 
mean response time analysis we have listed the path sets of the 




networks. The tasks follow any one path fiom processing for source to 
terminals system depending upon the allocation of path by the source 
system and it is assumed that arrivals of tasks occurs to as an ordinary 
poison process with mean X, and then they are served on the basis of 
first come first serve discipline. Each job of data packets has variable 
number of modules and their arrivals X at a time except when less then 
X is in the data packet. The amount of time the required for the service 
of any job is an exponentially distributed random variable. 


In this chapter, a model for mean response time is analyzed, 
which is based on graph and queueing theories for making quantitative 
analysis in networking. Then the all the path sets of network have been 
listed and the M’^/M/l bulk arrival queueing model is then applied for 
determining the mean response time of each path set. 


6.1 THE PROPOSED METHODS 


Let us consider a bridge, network, consist of four system namely. 
Si, S2, S3 and S4, which are inter connected by communicatioiL 


Terminal 


Source 





The system Si is defined as source system, while the system S3 is 
labeled as terminal system. The job enter from the source system and get 
released firom the terminal system after getting processed. The path sets 
of the network are listed below, along with then number of system in a 


5 3 number of system (N) = 2 

S2 - S3 number of system (N) = 3 

5 4 - S 3 number of system (N) = 3 


6.3 MEAN WAITING TIME 


The Kol mogrov equation are given below 


As t -> 00 then the steady state differential difference equations are 


Solving the equation with the help of generating function 



00 00 

C(Z)=2 CnZ"=2 CnZ"(IzI<l) 

n-0 


Then using from equation (6.3.3) 


(6.3.6) 


00 00 CO rt 

0 = -X2. P„z".^S P..,z"' + x2: X Pn-K U z“ (6.3.7) 

«=0 «=1 z n=0 «=1 S 




n 

Here we see that is the probability function for the sum of 

4=1 

oo n n 

the steady state system size and batch size ^ ^ Pn-k Ck zf’ =21 

Jk«l ' ^ 

± P^uZ^‘*=C(z)P(z) 

«=1 

Then equation (6.3.7) we get 

0 = -A.P(z)-p[P(z)-Po]+ ^ [P(z)-Po] + A.C(z)P(z) 


P(z) [x+m-M-^-C{z)]=\x?q - ^ Po 


P(z) 


P(z) 


P(z) 


P(z) 


{X + iS)z - fi- X^C(z) 


= pPo(l--) 

z 


//.Po(z-l) 

(X + ju)z-//- X.zC (z) 

/f.Po(z-l) 

[//(z-i)-;i.z(C(z)-ix 

-//•Po(z-l) 

/i(l-z)-A.z(l-C(z))J 


(Iz Ul) (6.3.8) 


and using the condition P (1) = 1 we get 
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1 = Liniz-^i P(z) = Lirri; 




Since numerator and denominator are both zero. So we use L’ Hospital 


rule 



(|a.z| <1) 


(63.12) 


C(z) 


z(l - a ) 


\ — az 

From equation (63.8) 

//.Po(l-z) 


P(z) = 


P(z) 




//(I -/?)(! -2) 


A.z 

f 

I 

1 

1 

/^(l - z) - xJ. 

z(l-a)" 


L l-cz.zj 


l-< 2 r.z _ 


//(l-p)(l-«r)(l-z) 


//{I - p)(l - az)(l - z) 


fi{\-z){\-a.z)-X.: 

. Mi-p)(i-«2) - 


1 


fO-Zl 

1-a.r 

(l-yO)(I-gz) 
/i(l-a.2)-2.2 \~\a-¥(}.-a)p\z 

1 


lj.{\-z)-l.. 


1 


1 - a .2 


P (z) = (1 - p)^ 


az 


\\-[a^-{\-a)p]z l-[a + (l-a)p ]2 

( CO «> 

P (z) = (1 - p) £[(« + (1 - a)p)zY - YX^ia + (1 - a)p)r z 


n 




So that, 

Pn= o-P)t^^+(i-«)pr-4a+o-«)pr’} 

Pn = 0 - p)[« + 0 - «)p]”~' 0 -a)p (n>0) 

6.4 LENGTH QUEUE : 

The line length of queue (Ls) of the system is given by 


Ls 




(63.13) 


(6.4 J) 


>r=\ 


Y «■(! - p)[« + (1 - (1 - «)P 




the waiting time of data packet of jobs in system for each path is as 


Therefore Mean waiting time of data packet of the job, is defined as 


6.5 MEAN SERVICE TIME 


Service of the system, which buUdly transmits the job fix>m one 
system to another, depends upon a batch of size K, which is 
exponentially distributed with the mean (i, j). Therefore the mean 
service time for each path, out going fi*om the transmitting systems 


given as 


6.6 MEAN RECEIVING TIME 


Similarly, for the service of the system, wlwse receives the job, it 
is assumed that the receiving time of this system with the mean, pr (i, j) 
obeys exponential law and is as. 


6.7 MEAN RESPONSE TIME : 


Mean Response Time (T^rt) is defined as the sum of the mean 
Waiting Time and Mean Service Time of processors, which transmit, 
plus the Mean Receiving Time of the processor, which receives, which 
receives the jobs i.e. 

Tmrt (i) = [Tw (i) + Ts (i) + Tr (i)] (6.7. 1 ) 

Where Tw (i), Ts (i) and Tr (i) are given in the equation (6.4.4.), (6.4.1), 
(6.5.1), respectively 

6.8 Computation of Numerical Results : 

For the computation of numerical results. Mean Response Time 
has been obtained for different value s of x i.e. different data packet size, 
which are shown through table 1 to 4. Let us assume the service rate for 
Si, S2, S3, S4 are 1.5, 2.0, 2.5 and 2.5 respectively and receiving rate for 
Si, S2, S3, S4 are 0.5, 1.0, 1.5 and 2.0 respectively. For the various path 
sets of the network the different values of batch size x and % are 1,2,3 E 
(x) are 0.5, 0.7, 0.8, and 0.9, 0.5, 0.7 for a. The Mean Waiting Time for 
each has been obtained. It is given in the following table-1 . 


TwO) 

M^n Waiting Time 


X=1 

X = 2 

X = 3 


0.0234 

0.0771 

0.0231 


0.0182 

0.0358 

0.0086 


0.0065 

0.0061 

0.0050 


I 












The Mean Receiving Time for each have been calculated. It is 
given in the following table -3 : 

Table -3 



TrO) 


Tr(1) 


Tr(2) 


Tr(3) 


Mean Receiving Time 


2.6666 


3.6666 


3.1666 


Finally the Mean Response Time, for the different values of 
batch size x and X i.e., 1,2,3 for each path have been obtained. It is 
shown in the following table-4 


Table -4 













Tmht (i) 

Mean Response Time 

Tmrt(1) 

3.7566 

3.8103 

3.7563 

f MRT (2) 

5.2514 

5.2690 

5.2418 

Tmrt(3) 

4.6397 

4.6393 

1 , 

1 4.6382 


6.9 DISCUSSION : 

The mean response time M’^/M/l bulk arrival qu^ue model to the 
various path sets for the bridge network as shown in figure- 1, has been 
obtained for the different value of batch size and arrival rate i.e. 1,2,3> 
respectively are shown in Table-4. Graphical representation of the batch 
size/arrivals rate vs mean response time has been shown in graph- 1. It 
can be stated, that the mean response time is lesser for tire path set-1. 
Further for the path set-2, the mean response time is higter as tte 
service rate is |X)or. In the case the path set-3, it is between the path set- 
1 and path set-2, shows that the mean response tiii» decrease with the 
increase service rate. The model discuss in this chapter, would be usefiil 
to the network system designer workii^ in the field of networking. The 
decision of routing strategy in the design of networic may also be 
facilitated by the present 
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